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state manipulations are described in terms 
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ditional dynamics implicit in quantum logic 
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mental issues in the proposal for trapped- 
ion quantum computation by J. I. Cirac and 
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cussed. Several possible decoherence mech- 
anisms are examined and what may be the 
more important of these are identified. 



Some potential applications for entangled 
states of trapped-ions which lie outside the 
immediate realm of quantum computation 
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1. Introduction 

A number of recent theoretical and experimental 
papers have investigated the ability to coherently 
control or "engineer" atomic, molecular, and optical 
quantum states. This theme is manifested in topics such 
as atom interferometry, atom optics, the atom laser, 
Bose-Einstein condensation, cavity QED, electromag- 
netically induced transparency, lasing without inver- 
sion, quantum computation, quantum cryptography, 
quantum-state engineering, squeezed states, and 
wavepacket dynamics. In this paper, we investigate a 
subset of these topics which involve the coherent 
manipulation of quantum states of trapped atomic ions. 
The focus will be on a proposal to implement quantum 
logic and quantum computation using trapped ions [1]. 



However, we will also consider related work on the 
generation of nonclassical states of motion and entan- 
gled states of trapped ions [2-39]. Many of these ideas 
have been summarized in a recent review [40]. 

Coherent control of spins and internal atomic states 
has a long history in NMR and rf /laser spectroscopy. 
For example, the ability to realize coherent "77 pulses" 
or "77/2 pulses" on two-level systems has been routine 
for decades. In much of what is discussed in this paper, 
we will consider entangling operations, that is, unitary 
operations which create entangled states between two or 
more separate quantum systems. In particular, we will 
be interested in situations where the interaction between 
quantum systems can be selectively turned on and off. 
For brevity, we will limit discussion to these types of 
operations in experiments which involve trapped 
atomic ions; however, many of the discussions, in 
particular those concerning single trapped ions, will 
also apply to trapped neutral atom experiments where 
the atoms can be treated as independent. The aspect of 
entangling operations is shared by atom optics and atom 
interferometry [41, 42] and, as described below, there 
are close parallels between the ion trap experiments and 
those of cavity QED [43]. 

Earlier experiments on trapped ions, where the zero- 
point of motion was closely approached through laser 
cooling, already showed the effects of nonclassical 
motion in the absorption spectrum [44^6]. These same 
effects can be used to characterize the average energy of 
the ion. More recent experiments report the generation 
of Fock, squeezed, coherent [21], and Schrodinger cat 
[47] states. These states appear to be of fundamental 
physical interest and possibly of use for sensitive detec- 
tion of small forces [26, 48]. For comparison, experi- 
ments which detect quantized atomic motion in optical 
lattices are reviewed by lessen and Deutsch [49]. Also, 
through the mechanism of Bose-Einstein condensation, 
which has recently been observed in neutral atomic 
vapors [50-53], a macroscopic occupation of a single 
motional state (the ground state of motion) is achieved. 

Simple quantum logic experiments have been carried 
out with single trapped ions [17]; the emphasis of future 
work will be to implement quantum logic on many ions 
[1]. The attendant ability to create correlated, or entan- 
gled, states of atomic particles appears to be interesting 
from the standpoint of quantum measurement [54] and, 
for example, for improved signal-to-noise ratio in 
spectroscopy of trapped ions (Sec. 3.4). 

Therefore, we will be particularly interested in study- 
ing the practical limits of applying coherent control 
methods to trapped ions for (1) the generation and 
analysis of nonclassical states of motion, (2) the imple- 
mentation of quantum logic and computation, and (3) 
the generation of entangled states which can improve 
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signal-to-noise ratio in spectroscopy. We will briefly 
describe the experimental results in these three areas, 
but the main purpose of the paper will be to anticipate 
and characterize decohering mechanisms which limit 
the ability to produce the desired final quantum states in 
current and future experiments. This is a particularly 
important issue for quantum computation where many 
ions (thousands) and coherent operations (billions) may 
be required in order for quantum computation to be 
generally useful. Here, we generalize the meaning of 
decoherence to include any effect which limits the 
purity of the desired final states. A fundamental source 
of decoherence will be the coupling of the ion's motion 
and internal states to the environment. Also important is 
induced decoherence caused by, for example, technical 
fluctuations in the applied fields used to implement the 
operations. This division between types of decoherence 
is arbitrary since both effects can be regarded as 
coupling to the environment; however, the division will 
provide a useful framework for discussion. As a unify- 
ing theme for the paper, we will find it useful to regard, 
as much as possible, the quantum manipulations we 
discuss in terms of quantum logic. Of course, the subject 
of decoherence is much broader than the specific 
context discussed here; the reader is referred to more 
general discussions such as the papers by Zurek 
[55,56,57]. 

The paper is organized as follows. In the next section, 
we briefly discuss ion trapping. In Sec. 3, we consider 
in somewhat more detail the three areas of application 
enumerated in the previous paragraph. Since cooling of 
the ions to their ground state of motion is a prerequisite 
to the main applications discussed in the paper, we out- 
line methods to accomplish this in the beginning of Sec. 
3. Section 4 is the heart of the paper; here, we attempt 
to identify the most important sources of decoherence. 
Section 5 briefly discusses some variations on proposed 
methods for realizing quantum logic in trapped ions. 
Section 6 suggests some additional applications of the 
ideas discussed in the paper and Sec. 7 provides a brief 
summary. 

Such a treatment seems warranted in that several lab- 
oratories are investigating the use of trapped ions for 
quantum logic and related topics; the authors are aware 
of related experiments being pursued at IBM, Almaden; 
Innsbruck University; Los Alamos National Labora- 
tory; Max Planck Institute, Garching; NIST, Boulder; 
and Oxford University. This analysis in this paper 
necessarily overlaps, but is also intended to complement, 
other investigations [58-66] and will, by no means, be 
the end of the story. We hope however, that this paper 
will stimulate others to do more complete treatments 
and consider effects that we have neglected. 



2. Trapped Atomic Ions 
2.1 Ions Confined in Paul Traps 

Due to their net charge, atomic ions can be confined 
by particular arrangements of electromagnetic fields. 
For studies of ions at low energy, two types of trap are 
typically used — the Penning trap, which uses a combi- 
nation of static electric and magnetic fields, and the Paul 
or rf trap which confines ions primarily through pon- 
deromotive forces generated by inhomogeneous oscil- 
lating fields. The operation of these traps is discussed in 
various reviews (see for example, Refs. [67]-[70]), and 
in a recent book by Ghosh [71]. For brevity, we discuss 
one trap configuration, the linear Paul trap, which may 
be particularly useful in the context of this paper. This 
choice however, does not rule out the use of other types 
of ion traps for the experiments discussed here. 

In Fig. 1 we show a schematic diagram of a linear 
Paul trap. This trap is based on the one described by 
Raizen et al. [72] which is derived from the original 
design of Drees and Paul [73]. It is basically a quadru- 
pole mass filter which is plugged at the ends with static 
electric potentials. A potential Vgcos fljt + U^ is applied 
between diagonally opposite rods, which are fixed in a 
quadrupolar configuration, as indicated in Fig. 1. We 
assume that the rod segments along the z direction are 
coupled together with capacitors (not shown) so that the 
rf potential is constant as a function of z ■ Near the axis 
of the trap this creates a potential of the form 







(VoCOS/2rf + Ur) 



1 +■ 



R' 



(1) 



where R is equal to the distance from the axis to the 
surface of the electrode. (Unless the rods conform to 
equipotentials of Eq. (1), this equation must be multi- 
plied by a constant factor on the order of 1; see for 
example, Ref. [72].) This gives rise to (harmonic) 
ponderomotive potentials in the x and y directions. To 
provide confinement along the z direction, static poten- 
tials Uo are applied to the end segments of the rods as 
indicated. Near the center of the trap, this provides a 
static harmonic well in the z direction 



<t>, = kUo 



^(x^-/) 



2q 



i(x^./) 



(2) 



where k is a geometric factor, m and q are the ion mass 
and charge, and w^ = (iKqUo/my^ is the oscillation 
frequency for a single ion or the center-of-mass (COM) 
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sections, O; (or a) will represent the harmonic oscillator 
lowering operator and qi will represent the normal mode 
coordinate for the ( th mode). The solution of Eqs. (3) to 
first order in at and second order in q, is given by 



li 



Ui{t) = A\ COS(Wit + (j)i) 



+ ^ cos(2/3tO 



1 + ^ cos(/2rO 



+ A f sin(w,? + 4>,) sin(/2rO ) , (4) 



*♦*♦» *»##*«♦***•*# ***•*•*••»•* 



Fig. 1. The upper part of the figure shows a schematic diagram of 
the electrode configuration for a linear Paul-rf trap (rod spacing = 
1 mm). The lower part of the figure shows an image of a string of 
"'Hg* ions, illuminated with 194 nm radiation, taken with a uv-sensi- 
tive, photon counting imaging tube [74]. The spacing between adja- 
cent ions is approximately 10 [Lva. The "gaps" in the string are occu- 
pied by impurity ions, most likely other isotopes of Hg*, which do not 
fluoresce because the frequencies of their resonant transitions do not 
coincide with those of the 194 nm ^Si;2 — > ^Pi;2 transition of "'Hg . 



where Ui = x or y. A, depends on initial conditions, and 



«, = A^, A^ [fl, + ^?/2]"^ 



(5) 



The large amplitude oscillation at frequency Wt is typi- 
cally called the "secular" motion. When a, « qf « 1 
and [/r — 0, if we neglect the micromotion (the terms 
which oscillate at Ht and IHt), the ion behaves as if it 
were confined in a harmonic pseudopotential 0p in the 
radial direction given by 



oscillation frequency for a collection of identical ions 
along the z direction. Equations (1) and (2) represent the 
lowest order terms in the expansion of the potentials for 
the electrode configuration of Fig. 1. When the size of 
the ion sample or amplitude of ion motion is 
comparable to the spacing between electrodes or the 
spacing between rod segments, higher order terms in 
and 0s become important. However for small oscilla- 
tions of the COM mode, which is relevant here, the 
harmonic approximation will be valid. In the x and y 
directions, the action of the potentials of Eqs. (1) and (2) 
gives the (classical) equations of motion described by 
the Mathieu equation 



g0p = 2 mwXx^ +y'-) , 



(6) 



d^ 
AC 

d^ 



Ux + 2qxC05{2^) 
Uy + 2gj,cos(2^) 



x = 



3^ = 0, 



(3) 



where ^ = [htll, a, = {AqlmD^){UJR^ -KUJzi), 
ay = -{4q/mn^){Ur/R^ + k[/o/Zo), q,= -qy = 2,VJ 
{O^mR^). The Mathieu equation can be solved in 
general using Floquet solutions. Typically, we will have 
Ui < qf « 1, ( £ {x, y}. (Keeping with the usual 
notation in the ion-trap literature, in this section, the 
symbols Ui and qi are defined as above. In all other 



where w, == qVo/(2"^OTmR^) = q^n^lilVl) is the radial 
secular frequency Wj. For most of the discussions in this 
paper, we will assume Ur = 0; however it may be useful 
in some cases to make Ur i= Oio break the degeneracy 
of the X and); frequencies. Figure 1 also shows an image 
of a "string" of '"Hg* ions which are confined near the 
z axis of the trap described in Ref. [74]. This was 
achieved by making Wr » w^ , thereby forcing the ions 
to the axis of the trap. The spacings between individual 
ions in this string are governed by a balance of the force 
along the z direction due to 0^ and the mutual Coulomb 
repulsion of the ions. Example parameters are given in 
the figure caption. 

When this kind of trap is installed in a high-vacuum 
apparatus, ions can be confined for days with minimal 
perturbations to their internal structure. Collisions with 
background gas can be neglected (Sec. 4.1.9). Even 
though the ions interact strongly through their mutual 
Coulomb interaction, the fact that the ions are localized 
necessarily means that the time-averaged value of the 
electric field they experience is zero; therefore electric 
field perturbations are small (Sec. 4.2.3). Magnetic field 
perturbations to internal structure are important; how- 
ever, the coherence time for superposition states of two 
internal levels can be very long by operating at fields 
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where the energy separation between levels is at an 
extremum with respect to field. For example, in a 'Be"^ 
(Penning trap) experiment operating in a field of 0.82 T, 
a coherence time between hyperfine levels exceeding 10 
min was observed [75, 76]. As described below, we will 
be interested in coherently exciting the quantized modes 
of the ions' motion in the trap. Here, not surprisingly, 
the coupling to the environment is relatively strong be- 
cause of the ions' charge. One measure of the decoher- 
ence rate is obtained from the linewidth of observed 
motional resonances of the ions; this gives an indication 
of dephasing times. For example, the linewidths of 
cyclotron resonance excitation in high resolution mass 
spectroscopy in Penning traps [77-79] indicate that 
these coherence times can be at least as long as several 
tens of seconds. Decoherence can also occur from tran- 
sitions between the ions' quantized oscillator levels. 
Transition times out of the zero-point motional energy 
level have been measured for single ^'^Hg* ions to be 
about 0. 15 s [44] and for single 'Be"^ ions to be about 
1 ms [45]. These relatively short times are, so far, unex- 
plained; however, it might be possible to achieve much 
longer times in the future (Sec. 4.1). 

In the linear trap, the radial COM vibration frequency 
(Or must be made sufficiently higher than the axial COM 
vibrational frequency w^ in order for the ions to be 
collinear along the z axis of the trap. This configuration 
will aid in addressing individual ions with laser beams 
and will also suppress rf heating (Sec. 4.1.5). To prevent 
zig-zag and other complicated shapes of the ion crystal, 
we require 0)^/(0^ > 1 for two ions, and Wr/Wj > 1 .55 for 
three ions. For L > 3 ions, the critical ratio {(aju)^)^ for 
linear confinement has been estimated analytically [80, 
81] yielding (w,/wjc ^ 0.13L°-^^ [60]. Other estimates 
are given in Refs. [82] [{wjo},), ^ 0.63L'"'^^] and [64] 
[((Wr/wjc — 0.59L''*'*'^]. An equivalent result is obtained 
if we consider that as the potential is weakened in the 
radial direction, ions in a long string which are spaced 
by distance Sc near the center of the string, will first 
break into a zig zag configuration. At the point where 
the ions break into a zig-zag, the net outward force from 
neighboring ions is equal to the inward trapping force. 
If we equate these forces, we obtain 



«.^ = 



877^0 ^ msl ' 



(7) 



where ^ is the Riemann zeta function. As an example, 
for 'Be"^ ions, m — 9 u (atomic mass units) and 
^c = 3 (jum, we must have (ajli; > 7.8 MHz to keep the 
ions along the axis of the trap. 

The equilibrium spacing of a linear configuration of 
trapped ions is not uniform; the middle ions are spaced 
closer than the outlying ions, as is apparent in Fig. 1 . The 



separation of two ions is ^2 = 2'"^, where 
s = (g^/477SoOT<wf )"^ is a length scale of ion-ion 
spacings; the adjacent separation of three ions is 
^3 = (5/4)"^^. For L large, estimates of the minimum 
separation of the center ions are given by ^c(^) 
=- IsL-^-^'' [60], 2.01 8^L-''-''' [61], 2.29L-''-''"^ [83], and 
\.92sL-^'^ [ln(0.8L)]"' [59, 81]. For typical trapping 
parameters, the ion-ion separations are on the order of a 
few (xm and the spatial spread of the zero-point vibra- 
tional wavepackets are on the order of 10 nm. Thus there 
is negligible wavefunction overlap between ions and 
quantum statistics (Bose or Fermi) play no role in the 
spatial wavefunction of an array. 

Of the 3L normal modes of oscillation in a linear trap, 
we are primarily interested in the L modes associated 
with axial motion because we will preferentially couple 
to them with applied laser fields. A remarkable feature 
of the linear ion trap is that the axial mode frequencies 
are nearly independent of L[ 1,60,6 1,84]. For two ions, 
the axial normal mode frequencies are at u)^ and \/2)U)^; 
for three ions they are w^, vSw^, and (5.8)^'^ w^. For 
L > 3 ions, the Lth axial normal mode can be deter- 
mined numerically [60,61,84]. 

2.2 Ion Motional and Internal Quantum States 

A single ion's motion, or the COM mode of a collec- 
tion, has a simple description when the ions are trapped 
in a purely static potential, which is the case for the axial 
motion in a Penning trap or the axial motion in the trap 
of Fig. 1. We will assume that the trap potentials are 
quadratic [Eqs. (1) and (2)]. This is a valid approxima- 
tion when the amplitudes of motion are small, because 
the local potential, expanded about the equilibrium point 
of the trap, is quadratic to a good approximation. In this 
case, motion is harmonic. An ion trapped in a pondero- 
motive potential [Eq. (6)] can be described effectively as 
a simple harmonic oscillator, even though the Hamilto- 
nian is actually time-dependent, so no stationary states 
exist. For practical purposes, the system can be treated 
as if the Hamiltonian were that of an ordinary, 
time independent harmonic oscillator [34,35,85-91] 
although modifications must be made for laser cooling 
[92]. The classical micromotion (the terms which vary 
as COS. flit and cos2flit in Eq. (4)] may be viewed, in the 
quantum picture, as causing the ion's wavefunction to 
breathe at the drive frequency fli. This breathing mo- 
tion is separated spectrally from the secular motion [(at 
frequencies u):, and uiy in Eq. (4)]. Since the operations 
we will consider rely on a resonant interaction at the 
secular frequencies, we will average over the compo- 
nents of motion at the drive frequency fli. Therefore, to 
a good approximation, the pseudopotential secular 
motion behaves as an oscillator in a static potential. 
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main consequence of the quantum treatment is that 
transition rates between quantum levels (Eq. (18), 
below) are altered [34, 35]; however, these changes can 
be accounted for by experimental calibration. In any 
case, for most of the applications discussed in this paper, 
we will be considering the motion of the ions along the 
axis of a linear Paul trap where this modification is 
absent. 

Therefore, the Hamiltonian describing motion of a 
single ion (or a normal mode, such as the COM mode, 
of a collection of ions) in the (th direction is given by 



ffosc = hoiihi, H 



[x,y,z\ 



(8) 



where hi = a/ at and a, and a, are the usual harmonic 
oscillator raising and lowering operators and we have 
suppressed the zero-point energy \l2Au)i. The operator 
for the COM motion in the z direction is given by 



z = Zo{a + a^) 



(9) 



where Zo = {^l2mu)^Y^ is the spread of the zero-point 
wavefunction and m is the ion mass. That is, zo = 
((0|z^|0))"^, where \n} is the nth eigenstate ("number" 
or Fock state) of the harmonic oscillator. For 'Be"^ ions 
in a trap where o)^/2t7= 10 MHz, we have zo = 7.5 nm. 
Therefore, a general pure state of motion for one mode 
can be written, in the Schrodinger picture, as 



^ m 



= Ec„^ 



\n). 



(10) 



where C„ are complex and the | n ) are time-independent. 
For applications to quantum logic, we will be interested 
in motional states of the simple form a 1 0) + 

;Sexp(-((WiO|l)- 

We will be interested in the situation where, at any 
given time, we interact with only two internal levels of 
an ion. This will be accomplished by insuring that the 
internal states are nondegenerate and by using resonant 
excitations to couple only two levels at a time. We will 
find it convenient to represent a two-level system by its 
analogy with a spin- 1/2 magnetic moment in a static 
magnetic field [93, 94]. In this equivalent representation, 
we assume that a (fictitious) magnetic moment 
/J, = fiuS, where S is the spin operator (S = 1/2), is 
placed in a (fictitious) magnetic field B = BoZ- The 
Hamiltonian can therefore be written 



H\r 



^WqS^ 



(11) 



where S^ is the operator for the z component of the spin 
and ft)o = - ixmBoI^- Typically, the internal resonant 
frequency will be much larger than any motional mode 



frequency, Wo » w^. We label the internal eigenstates 
|Mj) = |T) and |4-) representing "spin-up" and "spin- 
down" respectively, and for convenience, will assume 
/Am < so that the energy of the |T) state is higher than 
the |4-) state. A general pure state of the two-level system 
is then given by 



^i„ 



Qe'f |i) + Cte-'T|T) 



(12) 



where | Cj,p + | Ct| ^ = 1. Of course, the two-level system 
really could be a 5 = 1/2 spin such as a trapped electron 
or the ground state of an atomic ion with a single 
unpaired outer electron and zero nuclear spin such as 

''^Mg^ 

2.2.1 Detection of Internal States 

The applications considered below will benefit from 
high detection efficiency of the ion's internal states. 
Unit detection efficiency has been achieved in experi- 
ments on "quantum jumps" [95-99] where the internal 
state of the ion is indicated by light scattering (or lack 
thereof), correlated with the ion's internal state. (More 
recently, this type of detection has been used in spec- 
troscopy so that the noise is limited by the fundamental 
quantum fluctuations in detection of the internal state 
[100]. In these experiments, detection is accomplished 
with a laser beam appropriately polarized and tuned to 
a transition that will scatter many photons if the atom is 
in one internal state (a "cycling" transition), but will 
scatter essentially no photons if the atom is in the other 
internal state. If a modest number of these photons are 
detected, the efficiency of our ability to discriminate 
between these two states approaches 100 %. We note 
that for a string of ions in a linear trap, the scattered light 
from one ion will impinge on the other ions; this can 
affect the detection efficiency since the scattered light 
will, in general, have a different polarization. 

The overall efficiency can be explained as follows. 
Suppose the atom scatters N total photons if it is 
measured to be in state |4-) and no photons if it is mea- 
sured to be in state |T). In practice N will be limited by 
optical pumping but can be 10^ or higher [101]. Here, 
we assume that N is large enough that we can neglect its 
fluctuations from experiment to experiment. We 
typically detect only a small fraction of these photons 
due to small solid angle collection and small detector 
efficiency. Therefore, on average, we detect «d = T]iN 
photons where 7]^ « 1 is the net photon detection 
efficiency. If we can neglect background, then for each 
experiment, if we detect at least one scattered photon, 
we can assume the ion is in state |4-). If we detect 
no photons, the probability of a false reading, that is, the 
probability the ion is in state |4-) but we simply did not 
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detect any photons, is given by Pn(0) = 
(1 - Vdf =" exp(-nd). For n^ = 10, Pn(0) == 4.5 X 10-^ 
for Md = 100, Pjv(O) == 4 X 10'^^ Therefore for Md > 10, 
detection can be highly efficient. 

Detection of ion motion can be accomplished directly 
by observing the currents induced in the trap electrodes 
[77-79, 102-104]. However, the sensitivity of this 
method is limited by electronic detection noise. Because 
the detection of internal states can be so efficient, mo- 
tional states can be detected by mapping their properties 
onto internal states which are then detected (Sec. 3.2). 

2.3 Interaction With Additional Applied Electro- 
magnetic Fields 

2.3.1 Single Ion, Single Applied Field, Single Mode 
of Motion 

We first consider the situation where a single, period- 
ically varying, (classical) electromagnetic field propa- 
gating along the z direction is applied to a single trapped 
ion which is constrained to move in the z direction in a 
harmonic well with frequency co^. We consider situa- 
tions where fields resonantly drive transitions between 
internal or motional states and when they drive transi- 
tions between these states simultaneously (entangle- 
ment). If we assume that the internal levels 
are coupled by electric fields, then the interaction 
Hamiltonian is 



Hi = - fii- E(z, t) , 



(13) 



where /td is the electric dipole operator for the internal 
transition and E is from a uniform wave propagating 
along the z direction and polarized in the x direction, 
E = £'iXCOs(fe -cot+ (p), where co is the frequency, k is 
the wavevector 277/A, and A is the wavelength. In the 
equivalent spin- 1/2 analog, we assume that a traveling 
wave magnetic field propagates along the z direction, is 
polarized in the x direction [B = B\X cos(fe - (at + (^)], 
and interacts with the fictitious spin (/t = (XmS). 
Therefore, for the spin analog, Eq. (13) is replaced by 



Hi = - fjL ■ B{z,t) 



= hniS^ + S.) (e 



i(kz-o>t+(}>) 



+ e 



■i(kz-o>t + (}>)• 



•), 



(14) 



where Afl = - fjumBilA (or - ixJixIA for an electric 
dipole), S+ = Sx + iSy, S- = S^- iSy and z is given by 
Eq. (9). We will assume that the lifetimes of the levels 
are long; in this case, the spectrum of the transitions 
excited by the traveling wave is well resolved if O is 
sufficiently small. 

It will be useful to transform to an interaction picture 
where we assume Hq = Hi„t„nai + ^osc and yi„t„actio„ = Hi. 



In this interaction picture, if we make the rotating-wave 
approximation (neglecting exp(±i(<w+ <Wo)0 terms), 
the wavefunction can be written 



M,=l,t n = Q 



(15) 



where \M^) and \n) are the time-independent internal 
and motional eigenstates. In general, this wavefunction 
will be entangled between the two degrees of freedom; 
that is, we will not be able to write the wavefunction as 
a product of internal and motional wavefunctions. 
We have H ^ H\ = Ul{t)H^Uo{t) where Uo{t) = 
exp(-i(Ho/^)t), resulting in 

H'l = ^nS^exp(i[T](ae-''"'' + a^e""^') - St + cj)]) 

+ h.c . (8 = o) — oio) , (16) 

where 77 = kzo is the Lamb-Dicke parameter and 
S+, S-, a, and a^ are time independent. We will be 
primarily interested in resonant transitions, that is, 
where 8= co^(n' -n) where n' and n are integers. How- 
ever, since we want to consider nonideal realizations, we 
will assume 8 = (n' - n) C0z + ^ where \A \ « co^, O. If 
we can neglect couplings to other levels (see Sec. 4.4.6), 
transitions are coherently driven between levels |4-, n) 
and |T, n') and the coefficients in Eq. (15) are given by 
Schrodinger's equation ihV^/dt = HI ^ to be 

a,„ = -,-(>-l"'-"l)e«^'-*)/3„,„CT,„., (17) 

where /3„> is given by [105, 106] 

= /3exp[-77V2](M<!/M>!)"' r]^"'-'"L„J"'-"^(r]^) , (18) 

where n^(n^) is the lesser (greater) of n' and n, and L" 
is the generalized Laguerre polynomial 



l:(x) = 2 (-ir 



n + a\ X- 



n —ml ml 



(19) 



Since we will be particularly interested in small values 
of n and a, for convenience, we list a few values of 

l:(X). 
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L|!(Z) = 1, L?{X) =l-X, LUX)l-2X + ^ , L!{X) = 1 - 3Z + | Z' - ^Z" 

2 lb 

Li{X) = 1, L,'(Z) = 2 - Z, Ll{X) = 3 - 3Z + :^Z', LJ{X) = 4 - 6Z + 2Z' - 7Z' 

2 O 



Ll{X) = 1, Lf{X) = 3 - Z, L|(Z) = 6 - 4Z + iz', L^X) = 10 - lOZ + |z' - 7Z' 

2 2 o 



(20) 



Equations (17) can be solved using Laplace transforms. 
The solution shows sinusoidal "Rabi oscillations" 
between the states |T, n') and \-l, n), so over the 
subspace of these two states we have 



^(t) 



, Z„'„ \ A . / Z„'„ 



-2/|^e-<I'-*-fMsinfer 



Z„. 



.4 



.2,-^e-<l'-^-fl"'-"l)sinfer 
Z„',„ \ 2 



.,— sm( — rl+cos( — r 



^(0), 



(21) 



where Z„.,„ = (4^ + 4/2„^.,„)"^, 4 = w - Wq - («' - «)<«,, and 
^ is given by 



^=Cu\i,n} + Cu'\'^,n'} = 



.Ci,„. 



(22) 



For the resonance condition A = 0, Eq. (21) simplifies 
to 



^(t) 



cos/2„',„i 



_;g.[^+2l"'-"lisin/2„.,„i 



^(0) 



(23) 



When the atom starts in an eigenstate, for each value of 
n' -n, the phase factor </>+ 7r\n'-n\l2 can be chosen 
arbitrarily for the first application of Hi, however once 
chosen, it must be kept track of if subsequent applica- 
tions of Hi are performed on the same ion. For conve- 
nience, we can choose it to be zero, although in most of 
what follows we will include a phase factor as a re- 
minder that we must keep track of it. In these expres- 
sions, we assume /2„',„ to be constant during a given 
application time t; this condition can be relaxed as 
discussed in Sec. 4.3.2. A special case of interest is 
when the Lamb-Dicke criterion, or Lamb-Dicke limit, is 
satisfied. Here, the amplitude of the ion's motion in the 
direction of the radiation is much less than A/2 tt which 
corresponds to the condition (''?^raotion|A:^z^|''?Cotion)"^ 
« 1. This should not be confused with the less restric- 
tive condition where the Lamb-Dicke parameter is less 
than 1 (17 « 1); if the Lamb-Dicke criterion is satisfied. 



I then 17 « 1, but the converse is not necessarily true. If 
the Lamb-Dicke criterion is satisfied, we can evaluate 
/2„',„ to lowest order in 17 to obtain 

a> = a,„' = /2t,I"'-"I (M>!/n<!)"^ {\n' - n \ !)"' . (24) 

We will be primarily interested in three types of transi- 
tions — the carrier (n' = n), the first red sideband 
(n' = n -1), and the first blue sideband (n' = n+ 1) 
whose Rabi frequencies, in the Lamb-Dicke limit, are 
given from Eq. (24) by H, r]n"^n, and r]{n+ l)"'/2 
respectively. 

In general, the Lamb-Dicke limit is not rigorously 
satisfied and higher order terms must be accounted for 
in the interaction [16,21,106]. As a simple example, 
suppose ^(0) = |J')|«) and we apply radiation at 
the carrier frequency (8 = 0). From Eq. (23), the wave- 
function evolves as 
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■^(t) = cosn„j\i,n) - i&'l' sin/2„,„f |T, n) . (25) 

For n = 0, we have i7o,o = /2exp(- rj^l). The exponential 
factor in this expression is the Debye- Waller factor 
familiar from studies of x-ray scattering in solids; for a 
discussion in the context of trapped atoms see Ref. [106] 
and Sec. 4.4.5. This factor indicates that the matrix 
element for absorption of a photon is reduced due to the 
averaging of the electromagnetic wave (averaging of the 
e* factor in Eq. (14)) over the spread of the atom's 
zero-point wavefunction. 

As a second simple example, we consider 
^(0) = \-l,n) and 8 = + co^ (first blue sideband). Equa- 
tion (23) implies 

^(t) = cos/2„ + i,„f|4, n) + e'* sin/2„ + ,,„f|t, n+ 1). (26) 

At any time t i= rmT/(2fl„+,„) (m an integer), ^is an 
entangled state between the spin and motion. If the 
excitation is left on continuously, the atom sinusoidally 
oscillates between the state |i, n) and |T, n+ 1). This 
oscillation has been observed in Ref. [21] and is repro- 
duced in Fig. 2. 

When the Lamb-Dicke confinement criterion is met 
and when the radiation is tuned to the red sideband 
(8 = - Wj), we find (choosing <^ = - 77/2) 



Hi = ^rin(S+a + S.a^) 



(27) 



This Hamiltonian is the same as the "Jaynes-Cummings 
Hamiltonian" [107] of cavity QED [43], which 
describes the coupling of a two-level atom to a single 
mode of the (quantized) radiation field. The problem we 
have described here, the coupling of a single two-level 
atom to the atom's (harmonic) motion is entirely 



analogous; the difference is that the harmonic oscillator 
associated with a single mode of the radiation field in 
cavity QED is replaced by that of the atom's motion. 
The suggestion to realize this type of Hamiltonian (in 
the context of cavity-QED) with a trapped ion was 
outlined in Refs. [2], [3], and [7]; however its use was 
already employed in the g-2 single electron experiments 
of Dehmelt [108]. 

Driving transitions between the \-l) and |T) states will 
create entangled states between the internal and 
motional states since, in general, the Rabi frequency 
will depend on the motional states (Eq. (18)). This 
"conditional dynamics," where the dynamics of one 
system is conditioned on the state of another system, 
provides the basis for quantum logic (Sec. 3.3). 

In this section, we have assumed that the atom inter- 
acts with an electromagnetic wave (Eq. (14)), which 
will usually be a laser beam. However, the essential 
physics which gives rise to entanglement is that the 
atom's internal levels are coupled to its motion through 
an inhomogeneous applied field. In the spin- 1/2 analog, 
the magnetic moment fi couples to a magnetic field 
B = B(z,t)x, yielding the Hamiltonian 



H, = -fjL,B(z,t) = 



Biz-- 



^0,0 + ^ 



1 d^B 

z=o 2 dz 



Z^ + 



(28) 



where, as above, jx^ '^ S+ + S_ and z is the position 
operator. The key term is the gradient dBldz- From 
the atom's oscillatory motion in the z direction, it expe- 
riences, in its rest frame, a modulation of B at frequency 
o)^. This oscillating component of B can then drive the 
spin-flip transition. As a simple example, suppose B is 
static (but inhomogeneous along the z direction so that 




\i0m 




40 60 

time (fis) 



BO 



100 



Fig. 2. Experimental plot of the probability Pi(t) of finding a single 'Be* ion in the | i) state after first preparing 
it in the |i) |0) state and applying the first blue side band coupling (Eq. (16), for S = + m^) for a time t. If there 
were no decoherence in the system, Pi{t) should be a perfect sinusoid as indicated in Eq. (26). Decoherence 
causes the signal to decay as discussed in Sec. 3.2.1. The solid line is a fit to an exponentially decaying sinusoid 
as indicated in Eq. (43). Each point represents an average of 4000 observations [21]. 



267 



Volume 103, Number 3, May-June 1998 

Journal of Research of the National Institute of Standards and Technology 



dBldz 4" 0) and w^ is equal to the resonance frequency 
Wo of the internal state transition. In its reference frame, 
the atom experiences an oscillating field due to the 
motion through the inhomogeneous field. Since w = (Uo, 
this field resonantly drives transitions between the inter- 
nal states. Because this term is resonant, it is the domi- 
nant term in Eq. (28), so //i — - ix^{dB/dz)z =^(5+ + 
S2)(a + a'') — 5+0 + 5-0^, where the last equality ne- 
glects nonresonant terms. If the extent of the atom's 
motion is small enough that we need only consider the 
first two terms on the right hand side of Eq. (28), Hi is 
given by the Janes-Cummings Hamiltonian (Eq. 27)). 
This Hamiltonian is also obtained if B is sinusoidally 
time varying (frequency w), we satisfy the resonance 
condition S = w - Wo = - w^ , and we make the rotating- 
wave approximation. This situation was realized in the 
classic electron g-2 experiments of Dehmelt, Van Dyck, 
and coworkers to couple the spin and cyclotron motion 
[108]. Higher-order sidebands are obtained by consider- 
ing higher order terms in the expansion of Eq. (28). 

One reason to use optical fields is that the field gradi- 
ents (for example, o'/(o'z)[e*^] = de*^ can be large 
because of the smallness of A. Stated another way, 
single-photon transitions between levels separated by rf 
or microwave transitions, which are driven by plane 
waves, may not be of interest because k is small (A large) 
and exp((fe) — 1 which implies o'/(o'z)[e*^] — 0. This 
makes interactions which couple the internal and exter- 
nal states as in Eqs. (27) and (28) negligibly small. This 
is not a fundamental restriction because electrode struc- 
tures whose dimensions are small compared to the wave 
length can be used to achieve much stronger gradients 
than are achieved with plane waves. Microwave or rf 
transitions can also be driven by using stimulated- 
Raman transitions as discussed in Sec. 2.3.3 below. A 
second reason to use laser fields is they can be focused 
so that, to a good approximation, they interact only with 
a selected ion in a collection. 

The unitary transformations of Eqs. (21) and (23) 
form the basic operations upon which most of the 
manipulations discussed in this paper are based. In this 
section, they were used to describe transitions between 
two states labeled |4-)|n) and |T)|n'). In what follows, we 
will include other internal states of the atom which will 
take on different labels; however, the transitions 
between selected individual levels can still be described 
by Eqs. (21) and (23). Sequences of these basic opera- 
tions can be used to construct more complicated opera- 
tions such as logic gates (Sec. 3.3). 

2.3.2 State Dynamics Including Multiple Modes of 
Motion 

In what follows, we will generalize the interaction 
with electromagnetic fields to consider motion in all 3L 



modes of motion for L trapped ions. Here, as was as- 
sumed by Cirac and Zoller [1], we consider that, on any 
given operation, the laser beam(s) interacts with only the 
yth ion; however, that ion will, in general, have compo- 
nents of motion from all modes. In this case Eq. (14) for 
thejth ion becomes 

Hi, = na (S+j + S.J ) [e '* ■ "J -"'^-^i^ + kc] , (29) 

where we now assume k has some arbitrary direction. 
We will write the position operator of they th ion (which 
represents the deviation from its equilibrium position) 



Xj = ujx + UL+jy + U2L+jZ, jG.{l,2, . . .L} . (30) 

We can express the Uj in terms of normal mode coordi- 
nates qic(kG[l, 2, . . . 3L }) through the matrix Df , by 
the following relations [109] 



Up= ^ Dj^ qk, qk= 2 D^Up, qk = qkoittk + al), 

k=\ p=\ 

(31) 

where qi, is the operator for the kih normal mode and a^ 
and al are the lowering and raising operators for the A:th 
mode. We have assumed that all normal modes are 
harmonic, which is a reasonably good assumption as 
long as the amplitude of normal mode motion is small 
compared to the ion spacing. (For two ions, the axial 
stretch mode's frequency is approximately equal to 
Wz( V3 - 9{a^la^^) where a^ is the (classical) amplitude 
of one ion's motion for this mode and a^ is the ion 
spacing). Following the procedure of the last section, 
we take Ho to be the Hamiltonian of the 7 th ion's inter- 
nal states and all of the motional (normal) modes 



Hoj = hbioSjj + 2] ^0^1 



'ktik 



(32) 



where hic = aluk- In the interaction picture (and making 
the rotating wave approximation), we have H{j = 
UojHij Uqj where Uqj = e\pi-iiHoj/^)t, yielding 



Hy = ^fi S+j exp 



(2 Vi("ic e-'""" + al e'"") - ;(8f - <^,) 



+ h.c. 



(33) 



where Ti=(k ■ xDl + k ■ yDt' + k ■ zDl'-*^)qm- (For 
the linear trap case, motion will be separable in the 
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X, y, and z directions and rjl will consist of one term.) 
In this interaction picture, the wavefunction is given by 



qr 






,(f)|M,), |{«,}) 



(34) 



where the coefficients are slowly varying and | {«* }) are 
the normal mode eigenstates (we have used the short- 
hand notation [rik] = rii, 112, . . . , Msl). In analogy with 
the previous section, we will be primarily interested in 
a particular resonance condition, that is, where S — 
u)k{n'k -Ilk) and il is sufficiently small that coupling to 
other internal levels and motional modes can be ne- 
glected. In this case, Eqs. (21) and (23) apply to the 
subspace of states |sL)j|nt) and |T)j|nt) if we make the 
definitions 



^=^; = Ci,„li),k) + a,„Jt),|nO 






, A = S-itik- nk)(Ok 



(35) 



and 



xi,^„^ ^ {A^ + 4{ni,,„fy'', 



(36) 

The last expression is the Rabi frequency for particular 
values of the {Wp^t }. More likely, the other mode states 
(p¥= k) will correspond to a statistical distribution; this 
is discussed in Sec. 4.4.5. For the application to 
quantum logic (Sec. 3.3) the COM mode appears to be 
a natural choice since rjl will be independent of 7. The 
dependence of rjl on j for the other modes is not a 
fundamental problem, but requires accurate bookkeep- 
ing when addressing different ions. The values of rjl can 
be obtained from the normal mode coefficients as 
described by James [61]. 

2.3.3 Stimulated-Raman Transition 

As indicated in the discussion following Eq. (28), we 
want strong field gradients to couple the internal states 
to the motion. If the internal state transition frequency 
coo is small, one way we can achieve strong field gradi- 
ents is by using two-photon stimulated-Raman transi- 
tions [45, 48] through a third, optical level as indicated 
in Fig. 3. In this case, as we outline below, the effective 
Hamiltonian corresponding to that in Eq. (14) is 
replaced by 



where ki, k2 and cou, 0)^2 are the wavevectors and 
frequencies of the two laser beams and the resonance 
condition between internal states corresponds to 
\(Oli - (OliI = (Oo- Even if Wo is small compared to optical 
frequencies, |A:i -k2\ can correspond to the wavevector 
of an optical frequency by choosing different directions 
for kj and ^2; this choice can thereby provide the desired 
strong field gradients. 

In Fig. 3 we consider that a transition is driven 
between states |-i) and |T) through state |3) by stimu- 
lated-Raman transitions using plane waves. Typically, 
we consider coupling with electric dipole transitions in 
which case 



Ei = SiEiCos{ki • X - cout + (pi), jG{1, 2} 



(38) 



C0„ 




Ii> 



//i = ;?/2(5+ + 5_)[e'K*>-*^'- 



- (WLl - WL2)' + <W 



+ h.c.] , (37) 



Fig. 3. Schematic diagram relevant to stimulated-Raman transitions 
between internal states \-l) and |T). Two plane wave radiation fields 
couple to a third state |3). The radiation fields are typically at laser 
frequencies; they are characterized by frequencies and wavevectors 
ft)L; and A:,- , / E { 1 , 2 ) . The couplings are typically described by electric 
dipole matrix elements. For simplicity, we assume field 1 only cou- 
ples states \-l) and |3); and field 2 only couples states |T) and |3). In 
this diagram, we do not show the additional energy level structure of 
the 3 L modes of motion. 



For simplicity, we assume laser beam 1 has a cou- 
pling only between intermediate state |3) and state \-l). 
Similarly, laser beam 2 has a coupling only between 
state 1 3) and state |T). Not shown in Fig. 3 are the 
energy levels corresponding to the 3L motional modes. 
Laser detunings are indicated in the figure, so 
(On - {(Ol2 + S) = coo, and we assume Ar » S, [cok] 
where {wt} are the 3L mode frequencies. We will 



269 



Volume 103, Number 3, May-June 1998 

Journal of Research of the National Institute of Standards and Technology 



assume the Raman beams are focussed so that they 
interact only with thej th ion. In the Schrodinger picture, 
the wavefunction is written 



xp... 



M, = i, T, 3 [nt] = 



X exp[- i{(OM^ + ni<wi + n2<W2 + ■ ■ ■ n^L(03L)t^ 
X |M,),|{«3l}), 



(39) 



Since A-r is large, state |3) can be adiabatically elimi- 
nated in a theoretical treatment (see, for example, 
Refs. [32], [48], [110], and Sec. 4.4.6.2). If we assume 
the difference frequency is tuned to a particular reso- 
nance 8 = w^{nl -n^), we can neglect rapidly varying 
terms and obtain 



O 



= i^C 



T,ni, ...nj nji A '^T.n i, . . . nj n 



where 
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(40) 



(41) 



where g, = E,e {-lie ■ r\3)e\pi-i(j>i)/i2A), 
g2^E2e{^\s ■ r|3)exp(- i(^)/(2;l), r)i = (Ak ■ xDi + 
Ak ■ yD^*' + Ak ■ zDl''*')qm and AA: = A:, -Aj. The 
terms |g^|/zlR and |g^|/zlR are the optical Stark shifts of 
levels 1 1 ) and 1 2) respectively. They can be eliminated 
from Eqs. (40) by including them in the definitions of 
the energies for the |4-) and |T) states or, equivalently, 
tuning the Raman beam difference frequency 8 to 
compensate for these shifts. If the Stark shifts are equal, 
both the 1 4-) and |T) states are shifted by the same 
amount, and there is no additional phase shift to be 
accounted (Sec. 4.4.3). Equations (40) for stimulated- 
Raman transition amplitudes are the same as for the 
two-level system (Sec. 2.3.1) if we make the identifica- 
tions cf)! - <p2 ^^ cf) and AA: <=^ k. Although the 
experiments can benefit from use of stimulated-Raman 
transitions, for simplicity, we will assume single photon 
transitions below except where noted. 

Another advantage of using stimulated-Raman 
transitions on low frequency transitions, as opposed to 
single-photon optical transitions, is the difference 



frequency between Raman beams can be precisely 
controlled using an acousto-optic modulator (AOM) to 
generate the two beams from a single laser beam. If the 
laser frequency fluctuations are much less than Ar, 
phase errors on the overall Raman transitions can be 
negligible [111]. Other advantages (and some disadvan- 
tages) are noted below. 



3. Quantum-State Manipulation 

3.1 Laser Cooling to the Ground State of Motion 

As a starting point for all of the quantum-state manip- 
ulations described below, we will need to initialize the 
ion(s) in known pure states. Using standard optical 
pumping techniques [112], we can prepare the ions in 
the \-l) internal state. Laser cooling in the resolved side- 
band limit [106, 113] can generate the |n = 0) motional 
state with reasonable efficiency [44, 45]. This type of 
laser cooling is usually preceded by a stage of 
"Doppler" laser cooling [106,114,115] which cools the 
ion to an equivalent temperature of about 1 mK. For 
Doppler cooling, we have {n)> 1, so an additional stage 
of cooling is required. 

Resolved sideband laser cooling for a single, harmon- 
ically-bound atom can be explained as follows: For sim- 
plicity, we assume the atom is confined by a 1-D 
harmonic well of vibration frequency co^. We use an 
optical transition whose radiative linewidth y^ is rela- 
tively narrow, y,^i « co^ (Doppler laser cooling applies 
when yrad — w^). If a laser beam (frequency co) is inci- 
dent along the direction of the atomic motion, the bound 
atom's absorption spectrum is composed of a "carrier" 
at frequency coo and resolved frequency-modulation 
sidebands that are spaced by co^, that is, at frequencies 
0)0 + (n' - n)co^ (Sec. 2.3). These sidebands in the spec- 
trum are generated from the Doppler effect (like vibra- 
tional substructure in a molecular optical spectrum). 
Laser cooling can occur if the laser is tuned to a lower 
(red) sideband, for example, at co= coo- co^- In this case, 
photons of energy ^(coo - &>z) are absorbed, and sponta- 
neously emitted photons of average energy ^coo - R 
return the atom to its initial internal state, where R = 
(Mf'/2m = hu)R is the photon recoil energy of the atom. 
Overall, for each scattering event, this reduces 
the atom's kinetic energy by hu>^ if w^ » cor, a condition 
which is satisfied for ions in strong traps. Since 
wjw^ = yf where 17 is the Lamb-Dicke parameter, this 
simple form of sideband cooling requires that the Lamb- 
Dicke parameter be small. For example, in 'Be*, if the 
recoil corresponds to spontaneous emission from the 
313 nm 2p ^Pyi -^ 2s ^S\i2 transition (typically used for 
laser cooling), cor/2tt — 230 kHz. This is to be 
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compared to trap oscillation frequencies in some laser- 
cooling experiments of around 10 MHz [45]. Cooling 
proceeds until the atom's mean vibrational quantum 
number in the harmonic well is given by 
<n)mi„ =^ {y/2(o,f « 1 [106,115,116]. 

In experiments, we find it convenient to use two- 
photon stimulated Raman transitions for sideband cool- 
ing [45, 117], but the basic idea for, and limits to, 
cooling are essentially the same as for single-photon 
transitions. The steps required for sideband laser cooling 
using stimulated-Raman transitions are illustrated in 
Fig. 4. This figure is similar to Fig. 3, except we include 
the quantum states of the harmonic oscillator for one 
mode of motion. Part (a) of this figure shows how, when 
the ion starts in the \-l) internal state, a stimulated- 
Raman transition tuned to the first red sideband 
InL)!^) -^ |T)|n-l) reduces the motional energy by 
^0)^. In part (b), the atom is reset to the \-l) internal state 
by a spontaneous-Raman transition from a third laser 
beam tuned to the |T) ^ |3) transition. We assume that 
there is a reasonable branching ratio from state |3) to 
state I i), so that even if the atom decays back to level | T) 
after being excited to level |3), after a few scattering 
events, the atom decays to state \-l). If Wr « w^, step (b) 
accomplishes the transition |T)|n - 1) — > |4-)|n - 1) with 



high efficiency. Steps (a) and (b) are repeated until the 
atom is optically pumped into the jJ-)! 0) state. When this 
condition is reached, neither step (a) or (b) is active and 
the process stops. In this simple discussion, we have 
assumed the transition |4')|«) — > |T)|« - 1) is accom- 
plished with 100% efficiency. However since, in 
general, the atom doesn't start in a given motional state 
\n), and since the Rabi frequencies (Eq. (18)) depend on 
n, this process is not 100 % efficient; nevertheless, the 
atom will still be pumped to the |4')|0) state. The only 
danger is having the stimulated-Raman intensities 
and pulse time t adjusted so that for a particular n, 
O^ ^J = rmr {m an integer), in which case the atom is 
"trapped" in the \\')\n) level. This is avoided by varying 
the laser beam intensities from pulse to pulse; one par- 
ticular strategy is described in Ref. [17]. 

So far, laser cooling to the \n = 0) state has been 
achieved only with single ions [44, 45]; therefore an 
immediate goal of future work is to laser cool a collec- 
tion of ions (or, at least one mode of the collection) to the 
zero-point state. Cooling of any of the 3L modes of 
motion of a collection of ions should, in principle, work 
the same as cooling of a single ion. To cool a particular 
mode, we tune the cooling radiation to its first lower 
sideband. If we want to cool all modes, sideband cooling 




\i) 



(a) stimulated Raman (An = -1) 



(b) spontaneous Raman (An = 0) 



Fig. 4. Schematic diagram relevant to laser cooling using stimulated-Raman transitions. In (a), we 
show that when (Dli - ft'L2 = (^o- (^z, stimulated-Raman transitions can accomplish the transition \i) 
I") ~> \')\" — !)■ In the figure, the transition for « = 2 is shown. In (b), spontaneous-Raman transtions, 
accomplished with radiation tuned to the |T) — > |3) transition, pumps the atom back to the |T) state, 
thereby realizing the transition |T) |« - 1) — > |4-) \n — 1). When atomic recoil can be neglected, one 
application of steps (a) and (b) reduces the atom's motional energy by ^o)^ unless « = 0, in which case 
the atom is in it's motional ground state. 
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must be cycled through all 3L modes more than once, or 
applied to all 3L modes at once, since recoil will heat all 
modes. For the COM mode, the cooling is essentially 
the same as cooling a single particle of mass Lm ; how- 
ever, the recoil energy upon re-emission is distributed 
over the 3L-1 other modes. Other methods to prepare 
atoms in the \n = 0) state are discussed in Refs. [5], [10], 
[12]. Morigi et al. [118] show that it is not necessary to 
satisfy the condition cor « co^ (17 « 1) to achieve 
cooling to n = 0. 

3.2 Generation of Nonclassical States of Motion of 
A Single Ion 

We begin with a discussion of the generation of non- 
classical motional states of a single trapped ion. This 
seems appropriate because the other applications 
discussed in this paper incorporate similar techniques. 
Much of the original interest in nonclassical states 
of mechanical motion grew out of the desire to make 
sensitive detectors of gravitational waves using (macro- 
scopic) mechanical resonators [119, 120]. For example, 
parametric amplification of mechanical harmonic oscil- 
lations can lead to quantum mechanical squeezing of the 
oscillation. In the meantime, nonclassical states of the 
radiation field were observed [43]. The close relation- 
ship of these two problems was pointed out above: in 
quantum optics, H^^^ of Eq. (8) represents a single mode 
of the radiation field, and Hi of Eq. (14) represents the 
coupling between the (quantized) field and atom. The 
nonclassical states of motion considered here, such as 
squeezed states, are the direct analogs of the nonclassi- 
cal photon states in quantum optics. They appear to be 
of intrinsic interest because, as in cavity QED, they 
allow the rather complicated dynamics of the simple 
quantum system [described by the Hamiltonian in 
Eq. (27)] to be studied. Before discussing some methods 
to create nonclassical states, we consider one method 
for analyzing them. 

3.2.1 Population Analysis of Motional States 

As described below in this section, from the |-1-)|0) 
state, it is possible to coherently create states of the form 
I-1-) ^motion where ^„,otio„ is given by Eq. (10). One way 
we can analyze the motional state created is as follows 
[21]: To the state |-1-) ^motion, we apply radiation on the 
first blue sideband (n ' = n + 1 in Eq. (23)) for a time t. 
We then measure the probability Pi{t) that the ion is in 
the I -1-) internal state. In the experiments of Meekhof et 
al. [21], the internal state |-1-) is the 2s ^S,/2 (F = 2, 
Mf = 2) state of 'Be*, and | T) corresponds to the 2s ^Sm 
(1, 1) state as shown in Fig. 5. The |-1-) state is detected 
by applying nearly resonant cr*-polarized laser radiation 



between the |-1-) and ^^y2{F = 3, Mp= 3) energy levels. 
Because the only decay channel of the ^'Py2iF = 3, 
Mf = 3) state is back to the |-1-) state, this is a cycling 
transition, and detection efficiency is near 1 (Sec. 2.2. 1 ). 
The experiment is repeated many times for each value of 
T, and for a range of t values. We find 



Pi{T) =^(\ + Y, Pn e-^"^cos(2/3„,i,„T) 



(42) 



where P„ = \ Ci,„ p is the probability of finding the ion in 
state |-l-)|n). The phenomenological decay constants 7„ 
are introduced to model decoherence that occurs during 
the application of the blue sideband. The measured sig- 
nal Piij) can be inverted (Fourier cosine transform), 
allowing the extraction of the probability distribution of 
vibrational state occupation P„ . 




Internal state 
energy 



(2,-2) 



(2,-1) 



(2,0) 



(2.1) 



(2,2) 



I 



aux 



Fig. 5. Hyperfme levels of the 2s ^Si/2 ground state of 'Be* in a weak 
magnetic field (not to scale). The energy levels are designated by 
horizontal lines. Above the lines, the levels are represented by atomic 
physics labels (F , Mf) where F is the total angular momentum (elec- 
tron plus nuclear angular momentum) and Mf is the projection of the 
angular momentum along the magnetic field axis. The separation of 
Zeeman substates in the different F manifolds is approximately equal 
to 0.7 X 10"'BoHz where Bo is expressed is teslas. The separation of 
the F = \ and F = 2 manifolds is approximately 1.25 GHz at Bo = 0. 
For simplicity of notation, in most of the paper we make the identifi- 
cations \F = 2,Mf = 2)= \i), |1, 1) = |T), |2, 0)= |aux). 



3.2.2 Fock States 

In Fig. 2, we show an experimental plot [21] of the 
probability Piir) of finding the ion in the |-1-) internal 
state after first preparing it in 

the |-1-)|0) state, and applying the first blue sideband 
for a time t. From Eq. (23), we would expect 
Piij) = cos^/3|,ot; however, we clearly see the effects 



272 



Volume 103, Number 3, May-June 1998 

Journal of Research of the National Institute of Standards and Technology 



of some decoherence process which we can represent 
adequately by the first term in Eq. (42) 



Plir) 



1 + e"^'"'cos2/2i oT 



(43) 



In this experiment, we think the decoherence is not 
simply caused by fundamental (radiative) decoherence 
but has contributions from fluctuations in laser power 
(which cause fluctuations in i7i,o), fluctuations in trap 
drive voltage Vo (which cause fluctuations in co^), and 
fluctuations in the ambient magnetic field (which cause 
fluctuations in coq. 

Neglecting for the moment the effects of decoher- 
ence, we see that for times r = m77/(2i7i,o) {m an in- 
teger), the ion is in a nonentangled state (|4-)|0) or 
|T)|1)). Therefore, if the ion starts in the |~l)| 0) state, we 
can prepare the atom in the |T)| 1) state (the |n = 1 ) Fock 
state) by applying the blue sideband for a time 
T = 77/2i7i,o, a so called Rabi 77 pulse. For other times, 
the ion is in an entangled state given by Eq. (26). This 
operation and the analogous operation on the first red 
sideband will form key elements of quantum logic using 
trapped ions. 

We can generate higher-n Fock states of motion by a 
sequence of similar operations. For example, to generate 
the |~l)|2) state, we start in the |~l)|0) state, apply a 77 
pulse on the first blue sideband, followed by a 77 pulse 
on the first red sideband. This leads to the sequence 
|i)|0) -> |t)|l) -> |i)|2) (neglecting overall phase fac- 
tors). In a similar fashion, Fock states up to |n = 1 6) have 
been created [21]. Other methods for creating Fock 
states have been suggested in Refs. [5], [8], [12], and 
[13]. 

3.2.3 Coherent States 

We can also create coherent states of motion; these 
states are closest in character to classical states of mo- 
tion. This can be accomplished if the atom is subjected 
to a spatially uniform classical force, or any force 
derived from a potential -f{t) • z, where / is a real 
c- number vector. For an ion which starts in the |n = 0) 
state, this force creates a displacement leading to a co- 
herent state la) defined hy a \a) = a\a) where a is a 
complex number [121]. This classical force can be real- 
ized by applying an electric field which oscillates at 
frequency co^. For example, if we apply a (classical) 
electric field E{t) = zE^ sin{cot -<f>), the corresponding 
interaction Hamiltonian (in the interaction frame for the 
motion) is given by 

Hi = - qEMae''"'' + a^e'"=') sm{cot - </>) . (44) 
If we express the motional wavefunction as in Eq. (10), 



Schrodinger's equation yields for the coefficients of the 
wavefunction 



C„ = Q{\/n Cn-\ — fl'vn + 1 C„ + 1 , 



(45) 



where fli = -qE^zo&^'^l^l/i). Equivalently, in the inter- 
action picture for the motion, the Hamiltonian of Eq. 
(44) leads to the evolution operator 



U{t) 



, [(filOa^ - (filO'o] 



D(/2,0 , 



(46) 



where D is the displacement operator with argument fl\t 
[121]. 

We can achieve the same evolution if we superim- 
pose two traveling wave fields which drive stimulated- 
Raman transitions between different \n) levels of the 
same internal state, and we make the difference fre- 
quency between the Raman beams equal to the trap 
oscillation frequency. For example, assume the ion is 
subjected to two lasers fields given by Eq. (38), where 
a)Li - a)L2 — ("z « (Uo- The dynamics can be obtained 
following the analysis in Sec. 2.3.3, except we replace 
level |i)|n) (|T)|m')) with level \g)\n){\g)\n') where \g) 
can be any ground state which has a matrix element 
with level 1 3) . For the coefficients of Eq. (10), we find 



C„ 



?.|^ + 



^21- 



C„ - 2/ (n|i^exp[(i7(ae ' 



"0 - ('((WlI - <Wl2)^] + h.c. \n')Cn' 



(47) 



where Q = -g'gjIAR, gi = qEi{g\si ■ r \3) exp(- i4>i)l 
(2^), 17 = {ki -k-i) ■ zzo, and /Ir = co^g - «li- The first 
term on the right side of Eq. (47) corresponds to a Stark 
shift of level g ; this Stark shift can be absorbed into the 
definition of the ground state energy (see, for example. 
Sec. 4.4.6.2). If this is done, the same equations for the 
C„ [Eqs. (47)] are obtained from the Hamiltonian (in the 
oscillator interaction picture) 



HI = ^Oexpl i[r}{ae "^'' + aV"^') - («li - «l2)^] 

+ h.c. , (48) 

In the Lamb-Dicke limit, if we choose the resonance 
condition where Wli - <Ul2 = <Uz and assume i7« co^, 
Eqs. (47) are the same as Eqs. (45) where Hi = rjfl. 
Therefore, in the Lamb-Dicke limit, the applied laser 
fields act like a uniform oscillating electric field which 
oscillates at frequency Wli - <Ul2- This can be under- 
stood if we consider that the two laser fields give rise 
to an optical dipole force which is modulated in such 
a way to resonantly excite the ion motion. To see 
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this, assume for simplicity that, in Eq. (38), £^=£2 = Eo, 
Si = S2 = s, and ki-k2 = z |AA:|. It is useful to write the 
total electric field as 

E = Et+E2= slE{t) 

X cos((A:i + ^2) • z/2 - wf + ^) , (49) 

where w = (<Wli + <Wl2)/2, 4> = ((pi + 02)/2, and E(t) is a 
slowly varying function 



E(t) = Eocos\—z 2 ^ + 2 ) ■ ^^^^ 



On a time scale long compared to 1/Ar but short com- 
pared to 1/(<Wli - <Wl2), the atom experiences a nonreso- 
nant electric field of amplitude E(t) which is nearly 
constant in time. If we consider coupling of this electric 
field between the ground state and state |3) (for example, 
see Sec. 4.4.6.2), this electric field leads to a spatially- 
dependent Stark shift of the ground state equal to 



4;l|g(z,0pMR where 



\8(z,tf = \ 



Ak 



«L1 



■(Ol2 <t>i 



0: 



g£o(g|g-r|3)e--^ 



(51) 



This Stark shift leads to an optical dipole force 
122-124] F^=- d{AEstM\)ldz. On a longer time scale, 
this dipole force is modulated at frequency Wli - ui-li 
which can resonantly excite the ion's motion when 
(wli - oi-Li) = u)z- This leads to Eqs. (45). When ki and 
k2 are both directed along the z axis (but in opposite 
directions), the dipole force potential can be viewed as 
a "moving standing wave" in the z direction which slips 
over the ion and whose accompanying dipole force reso- 
nantly excites the ion's motion [125]. Both methods 
have been used to excite coherent states in Refs. [21 ] and 
[47]. Other methods for generating coherent states are 
suggested in Refs. [4] and [48]. 

In the experiment of Ref. [47], a dipole force oscillat- 
ing at the ion oscillation frequency was created with 
particular polarizations of the laser fields. This led to a 
force which was dependent on the ion's internal state, 
enabling the generation of entangled "Schrodinger-cat" 
states of the form ^= (|i)|ae'*) + |t)|ae-'*))/V2. 

3.2.4 Other Nonclassical States 

When (wli - ui-lt) = 2co,a similar analysis shows that 
the moving standing wave potential discussed in the last 



section has a component which acts like a parametric 
excitation of the ion's harmonic well at frequency 2co 
[21]. This can produce quantum mechanical squeezing 
of the ion's motion. Squeezing could also be achieved by 
amplitude modulating Uo at frequency 2co^ by a nonadi- 
abatic change in the trap spring constant [48], or 
through a combination of standing and traveling wave 
laser fields [4]. A quantum mechanical treatment of the 
motion in an rf trap shows the effects of squeezing from 
the applied rf trapping fields [34, 35, 91]. More general 
nonlinear effects in the interaction can lead to higher 
"nonlinear coherent states" as discussed by de Matos 
Filho and Vogel [126]. 

Other methods for generation of Schrodinger-cat like 
states in ions are suggested in Refs. [8], [19], [29], [31], 
[33], [36], and [38]. Additional nonclassical states are 
investigated theoretically by Gou and Knight [23], Gou 
et al. [37], and Gerry et al. [37]. Schemes which can 
generate arbitrary states of the single-mode photon field 
[127,128] can also be applied directly to generate 
arbitrary motional states of a trapped atom and perform 
quantum measurements of an arbitrary motional observ- 
able [129]. A scheme which can generate arbitrary 
entanglement between the internal and motional levels 
of a trapped ion is discussed by Kneer and Law [130]. 

The procedure for analyzing motional states outlined 
in Sec. 3.2.1 yields only the populations of the various 
motional states and not the coherences. Coherences 
must be verified separately [21, 47]. The most complete 
characterization is achieved with a complete state 
reconstruction or tomographic technique; a description 
of how this has been implemented to measure the 
density matrix or Wigner function for trapped atoms is 
given in Leibfried et al. [131, 132]. These experiments 
represent the first measurement of negative values of the 
Wigner function in position-momentum space. Wigner 
functions for free atoms have also been recently deter- 
mined experimentally by Kurtsiefer et al. [133]. Other 
methods for trapped atoms have been suggested in Refs. 
[15], [22], [24], [27], [34], and [35]. These techniques 
can also be extended to characterize entangled motional 
states [39] and states which are entangled between the 
motional and internal states. 

3.3 Quantum Logic 

Significant attention has been given recently to the 
possibility of quantum computation. Although this field 
is about 15 years old [134-138], interest has intensified 
because of the discovery of algorithms, notably for 
prime factorization [139-142], which could provide 
dramatic speedup over conventional computers. 
Quantum computation may also find other applications 
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[142-152]. Schemes for implementing quantum compu- 
tation have been proposed by Teich [153], Lloyd [143, 
144], Berman et al. [154], DiVincenzo [155, 156], Cirac 
and Zoller [1], Barenco et al. [157], Sleator and 
Weinfurter [158], Pellizzari et al. [159], Domokos et al. 
[160, 161] Turchette et al. [162], Lange et al. [163], 
lorma and Stenholm [164], Gershenfeld and Chuang 
[165], Cory et al. [166], Privman et al. [167], Loss and 
DiVincenzo [168], and Bocko et al. [169]. In this paper, 
we focus on a scheme suggested by Cirac and Zoller [1], 
which uses trapped ions. Since, in general, any quantum 
computation can be composed of a series of single-bit 
rotations and two-bit controlled-not operations 
[140,155,156,170,171], we will focus our attention on 
these operations. 

In the parlance of quantum computation, we say that 
two internal states of an ion can form a quantum bit or 
"qubit" whose levels are labeled |0) and |l) or, equiva- 
lently, \-l) and |T). Single-bit rotations on ion j can be 
characterized by the transformation (Eq. (23) for n' = n) 



R{e,ci>){Ci,\i)j + Cy\r)j) 



cos(0/2) 



(e*"^sin(0/2) 



-(e-''^sin(0/2) cos(0/2) 



X 



Cij. 



(52) 



In the spin- 1/2 model, this transformation is realized by 
application of a magnetic field BJl which rotates at 
frequency coo and in the same sense as /u, which is 
applied slong the direction xcos<f> -ysincf) in the rotat- 
ing frame. This is equivalent to application of the field 
BiXCOs{kz-cot+<f>) in Eq. (14). In this expression, 
9 = IflnJ is the angle of rotation about the axis of 
this field. For 9 = tt and 4> = 0, R{9, <f>) is a logical 
"not" operation (within an overall phase factor). 
R{tt/2, - 77/2) (plus a rotation about z) is essentially a 
Hadamard transform. 

A fundamental two-bit gate is a controlled-not (CN) 
gate [142,155,156,157]. This provides the transforma- 
tion 



\e,) \e2) -^ \et) \e,®e2) , 



(53) 



where Si, ^iSfO, 1} and @ is addition modulo 2. 
Although Eq. (53) is written in terms of eigenstates, the 
transformation is assumed to apply to arbitrary superpo- 
sitions of states |£i)|£2)- In this expression, Si is the 
called the control bit and S2 is the target bit. If Si = 0, the 



target bit remains unchanged; if ^i = 1, the target bit 
flips. 

A spectroscopy experiment on any four-level quan- 
tum system, where the level spacings are unequal, shows 
this type of logic structure if we make the appropriate 
labeling of the levels. For example, we could label these 
four levels as in Eq. (53). If we tune radiation to the 
|l, 0) — > |l, 1) resonance frequency and adjust its dura- 
tion to make a 77 pulse, we realize the logic of Eq. (53). 
Similarly, an eight-level quantum system with unequal 
level spacings realizes a Toffoli gate [142], where the 
flip of a third bit is conditioned upon the first two bits 
being I's — and so on (see Sec. 5.2). This basic idea can 
be applied to molecules composed of many interacting 
spins such as in the proposals of Gershenfeld and 
Chuang [165] and Cory et al. [166]. For quantum 
computation to be most useful, however, we need to 
perform a series of logic operations between an arbi- 
trary number of qubits in a system which can be scaled 
to large numbers, such as the scheme of Cirac and Zoller 

[1]. 

Another type of fundamental two-bit gate is a phase 
gate, which could take the form 



\e,)\e2) -^ c'^'''''\e,)\e2) 



(54) 



This type of conditional dynamics has been demon- 
strated in the context of cavity QED [162,172] and for 
a trapped ion [17] (step (lb) below). 

The Cirac/Zoller scheme assumes that an array of 
ions are confined in a common ion trap. The ions 
are held apart from one another by mutual Coulomb 
repulsion as shown, for example, in Fig. 1. They can be 
individually addressed by focusing laser beams on the 
selected ion. Ion motion can be described in terms of 
normal modes of oscillation which are shared by all of 
the ions; a particularly useful mode might be the COM 
axial mode. When quantized, this mode can form the 
"bus qubit" through which all gate operations are per- 
formed. We first describe how logic is accomplished 
between this COM mode qubit and the internal-state 
qubit of a single trapped ion. In particular, the transfor- 
mation in Eq. (53) has been realized for a single trapped 
ion [17]. In that experiment, performed on a trapped 
'Be* ion, the control bit was the quantized state of one 
mode of the ion's motion (labeled the x mode). If 
the motional state was \n = 0), this was taken to be a 
1^1 =0) state; if the motional state was \n = 1), this was 
taken to be a l^i = 1) state. The target states were two 
ground-hyperfine states of the ion, the \F = 2, Mp = 2) 
and |F = 1, Mf = 1) states, labeled \i) and |T) 
(Fig. 5), with the identification here \-l) <=^ l^i = 0) 
and |T) <=^ \e2 = 1). Transitions between levels were 
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produced using two laser beams to realize stimulated- 
Raman transitions. The wavevector difference ki - k2 
was chosen to be aligned along the x direction. The CN 
operation between these states was realized by applying 
three pulses in succession: 

(la) A 77/2 pulse (fit = tt/4 in Eq. (25), where we 
assume i7o,o = /2i,i = O) is applied on the carrier 
transition. For a certain choice of initial phase, 
this corresponds to the operator V"^ ('n"/2) of 
Cirac and Zoller [1]. 

(lb) A 277 pulse is applied on the first blue sideband 
transition between levels |T) and an auxiliary 
level |aux) in the ion (the \F = 2, Mp = 0) level in 
'Be"^; see Fig. 5). This operator is analogous to the 
operator Un'' of Cirac and Zoller [1]. This opera- 
tion provides the "conditional dynamics" for the 
CN operation. It changes the sign of the |T) 
|m = 1) component of the wavefunction but leaves 
the sign of the |T) |« = 0) component of the wave- 
function unchanged; that is, the sign change is 
conditioned on whether or not the ion is in the 
\n =0) or \n = 1) motional state. Therefore, this 
step is the phase gate of Eq. (54) with = 77, 
where we make the identifications (ei = 0,1) <=4> 
(n = 0,1) and (gj = 0,1) o (internal state = -t,t). 

(Ic) A 77/2 pulse is applied to the spin carrier transi- 
tion with a 180° phase shift relative to step (a). 
This corresponds to the operator y"^(-77/2) of 
Cirac and Zoller [1]. 

Steps (la) and (Ic) can be regarded as two resonant 
pulses (of opposite phase) in the Ramsey separated-field 
method of spectroscopy [173]. If step (lb) is active 
(thereby changing the sign of the |T) \n= 1) component 
of the wavefunction), then a state change (spin flip) is 
induced by the Ramsey fields. If step (lb) is inactive, 
step (Ic) reverses the effect of step (la). 

Instead of the three pulses (la - Ic above), a simpler 
CN gate scheme between an ion's internal and motional 
states can be achieved with a single laser pulse, while 
eliminating the requirement of the auxiliary internal 
electronic level [174], as described below. These simpli- 
fications can be important for several reasons: 

(1) Several popular ion candidates, including ^''Mg"^, 
"''Ca% ''St\ "'Ba\ ^'^Yb% ''^Yb% and "^Hg% do 
not have a third electronic ground state available 
for the auxiliary level. These ions have zero nu- 
clear spin with only two Zeeman ground states 
(M^ = -i,T). Although excited optical metastable 
states may be suitable for auxiliary levels in some 
of these ion species, use of such states places 
stringent requirements on the frequency stability 
of the exciting optical field to preserve coherence 
(see Sec. 4.4.3). 



(2) The elimination of an auxiliary ground state level 
removes "spectator" internal atomic levels, which 
can act as potential "leaks" from the two levels 
spanned by the quantum bits (assuming negligible 
population in excited electronic metastable 
states). This feature may be important to the 
success of quantum error-correction schemes 
[142, 175-187] which can be degraded when 
leaks to spectator states are present [188]. 
(Specific error-correction schemes for ions are 
suggested in Refs. [182] and [187].) 

(3) The elimination of the need for an auxiliary level 
can dramatically reduce the sensitivity of a CN 
quantum logic gate to external magnetic fields 
fluctuations. It is generally impossible to find 
three atomic ground states whose splittings are all 
magnetic field insensitive to first order. However, 
for ions possessing hyperfine structure, the tran- 
sition frequency between two levels can be made 
magnetic field independent to first order at partic- 
ular values of an applied magnetic field (see Sec. 
4.2.2). 

(4) Finally, a reduction of laser pulses simplifies the 
tuning procedure and may increase the speed of 
the gate. For example, the gate realized in Ref. 
[17] required the accurate setting of the phase and 
frequency of three laser pulses, and the duration 
of the gate was dominated by the transit time 
through the auxiliary level. 

The CN quantum logic gate can be realized with a single 
pulse tuned to the carrier transition which couples the 
states \n)\-l) and \n) |T) with Rabi frequency /2„,„ [see 
Eqs. (18), (36), and (41)]. Considering only a single 
mode of motion. 



/2„,„ = /2 1 (m I e 



irj(a + a ) 



= fh-^''^ £e„iv') , (55) 



where 2!„(r]^) = L°„(r]^) (Eq. (18)). Specializing to the 
|« = 0) and |« = 1) vibrational levels relevant to quantum 
logic, we have 






(56) 



The CN gate can be achieved in a single pulse by 
setting 7] so that fli,i/flo,o = (2k + l)/2m, with k and m 
positive integers satisfying m > A: > 0. Setting iliJila.o 
= 2m/(2k+ 1) will also work, with the roles of the |0) 
and |l) motional states switched in Eq. (53). By driving 
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the carrier transition for a duration t such that 
/2i_iT = (k + \l2)Tr, or a 'V-pulse" (mod 2Tr) on the 
\n) = |l) component, this forces flo,oT = rmr. Thus the 
states |-l)|l) and |T)|l) are swapped, while the states 
|4-)|0) and |T)|0) remain unaffected. The net unitary 
transformation, in the {4-0, TO, 4-1, Tl } basis is 



collection of L ions. Here, we choose the particular ion 
oscillator mode to be a COM mode of the collection. 
Specifically, to realize a controlled-not Cc,t between two 
ions c = control bit, t = target bit), we first assume the 
COM mode is prepared in the zero-point state. The 
initial state of the system is therefore given by 



1 














1 

















/e'*(- 1) 








!e-'*(-l)*-" 






\[r — 



(57) 



This transformation is equivalent to the reduced CN of 
Eq. (53), apart from phase factors which can be elimi- 
nated by the appropriate settings of the phase of subse- 
quent logic operations [157]. 

The "magic" values of the Lamb-Dicke parameter 
which allow the above transformation satisfy 
idiirf') = \ - 7]^ = {2k + l)/2m, and are tabulated in 
Table I of Ref. [174] for the first few values. [For rf 
(Paul) trap confinement along the COM motional mode, 
the Rabi frequencies of Eqs. (55) and (56) must be 
altered to include effects from the micromotion at the rf 
drive frequency /2r. In the pseudopotential approxima- 
tion (o) « /2r), this correction amounts to replacing the 
Lamb-Dicke parameter rj in this paper by rjll - co/ 
(2V2/2r)], as pointed out by Bardroff et al. [34, 35]. 
However, there is no correction if the COM motional 
mode is confined by static fields (such as the axial COM 
mode of a linear trap.)] It may be desirable for the 
reduced CN gate to employ the |n) = |2)or|n) = |3) state 
instead of the \n) = |l) state for error-correction of 
motional state decoherence [182]. In these cases, the 
"magic" Lamb-Dicke parameters satisfy ^liv^) = 
1 -277^ + 77 ^2 = {2k + l)/2m for quantum logic with 
\n) = |0) and |2), or ^^{r]^) = 1 - 3 77' + 3r]V2 - 7]% = 
{2k + l)/2m for quantum logic with \n) = |0) and |3). 

This scheme places a more stringent requirement on 
the accuracy of O and rj , roughly by a factor of m . In 
the two-photon Raman configuration (Sec. 2.3.3), the 
Lamb-Dicke parameter 77 = I/IA: |zo can be controlled by 
both the frequency of the trap (appearing in zo) and by 
the geometrical wavevector difference Ak of the 
two Raman beams. Accurate setting of the Lamb-Dicke 
parameter should therefore not be difficult. Both 
CN-gate schemes are sensitive to excitation in other 
modes as discussed in Sec. 4.4.5. 

The CN operations between a motional and internal 
state qubit described above can be incorporated to 
provide an overall CN operation between two ions in a 



2 2 

Ml =i,T M2 = 4.,T 



2 c, 

Ml = -I, t 



M], M2, . . .Ml 



X |M,),|M2)2 . . . iMOi |0) . (58) 



Cc,t can be accomplished with the following steps: 



(2a) 



(2b) 



(2c) 



Apply a •;7-pulse on the red sideband of ion c. 
This accomplishes the mapping (a|4-)c + 
13 1 t)c) 1 0) ^ I l)c(a 1 0) + ;8 1 1 ), and corresponds to 
the operator U'm of Cirac and Zoller [1]. We note 
that in the NIST experiments [17], we prepare the 
state (a|4-) + (3|t))|0) from the |4-)I0) state using 
the carrier transition. We can then implement the 
mapping (a |i) + (3|t))|0)^|i)(a|0) + ;8|l)) by 
applying a •77 -pulse on the red sideband. This is 
the "keyboard" operation for preparation of arbi- 
trary motional input states for the CN gate of 
steps 1 a - 1 c above. Analogous mapping of inter- 
nal state superpositions to motional state super- 
positions are reported in Refs. [47], [131], and 
[132]. 

Apply the CN operation (steps la - Ic or, the 
single carrier pulse described above) between the 
COM motion and ion t. 
Apply the inverse of step (2a). 



Overall, Cc t provides the unitary transformation (in the 

{|i),|iX, |4|tX, lUl^X, |t),|TX} basis) 



^c.t — 



10 
10 
1 
10 



(59) 



which is the desired logic of Eq. (53). Effectively, Cc,t 
works by mapping the internal state of ion c onto the 
COM motion, performing a CN between the motion and 
ion t, and then mapping the COM state back onto ion c. 
The resulting CN between ions c and t is the same as the 
CN described by Cirac and Zoller [1], because the oper- 
ations V"\e) and U'^ commute. 
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A third possibility, which also uses only one internal 
state transition on each ion, is the following. We employ 
two nondegenerate motional modes, which we label 
here as 1 and 2. These might be the COM modes in two 
different directions. We first map the internal state infor- 
mation from two qubits j and k onto the separate 
motional modes (which are both initially in the \n = 0)i 
\n = 0)2 zero-point state). This can be accomplished as 
described in step (2a) above. We then apply a condi- 
tional phase gate (Eq. (54) with <^ = 77) to the two 
motional modes. This could be accomplished by driving 
a 277 transition on a second order red sideband, at 
frequency Wo - <wi - <W2, on a particular (extra) ion "g" 
which is initially in the |4-)g state. This ion is not used 
to store information; it is only used for this one particu- 
lar purpose. This would be followed by operations 
which map the motional states back onto the internal 
states of ions 7 and k (like step (2c) above). Overall, this 
provides a phase gate (Eq. (54) with cf) = 77) between 
ions j and k. To make a CN gate between ions j and k, 
we need to precede the above operations with a 77/2 
pulse on the internal state of iony (or k) and follow the 
above operations with a 77/2 pulse on the internal state 
of iony (or k). 

In this section, we have assumed that each ion can be 
addressed independently. Also, since very many such 
operations will be desired for a quantum computer, the 
accuracy or fidelity of these operations is of crucial 
importance. These issues are confronted in Sec. 4. As 
noted in Sec. 2.3, in each separate operation involved in 
a quantum computation, such as application of the red 
sideband in step 2(a) to ion j, a definite phase of the 
applied fields is assumed. This phase for each ion can be 
chosen arbitrarily for the first operation, but upon 
successive applications of the same operation to the 
same ion, it must be held fixed, or at least be known, 
relative to the initial phase. An exception to this is appli- 
cation of 277 -pulses as in step 1(b) where the phase of 
the fields does not enter into the final result of the 
operation. 

3.4 Entangled States for Spectroscopy 

A collection of atoms, whose internal states are entan- 
gled through the use of quantum logic, can improve the 
quantum-limited signal-to-noise ratio in spectroscopy. 
Compared to the factorization problem, this application 
has the advantage of being useful with a relatively small 
number of ions and logic operations. For example, for 
high-accuracy, ion-based frequency standards [74,189, 
190,191], use of a relatively small number of trapped 
ions (L < 100) appears optimum. As outlined here, the 
states involving L ions which are useful for spectroscopy 
and frequency standards can be generated with L logic 



gates. For L — 100, a significant improvement in perfor- 
mance in atomic clocks could be expected. 

In spectroscopy experiments on L atoms, in which the 
signal relies on detecting changes in atomic populations, 
we can view the problem in the following way using the 
spin- 1/2 analog for two-level atoms. The total angular 

momentum of the system is given by J =^ Si where Si 

i=l 

is the spin of the ith atom. The basic task is to measure 
Wo, the frequency of transitions between the |4-) and |T) 
states Eq. (11). We first prepare an initial state for the 
spins. We assume spectroscopy is performed by apply- 
ing (classical) fields of frequency cor for a time Tr 
according to the method of separated fields by Ramsey 
[173]. The same field is applied to all atoms. After 
applying these fields, we measure the final state popula- 
tions; that is, we detect, for example, the number of 
atoms Ni in the |4-) state. In the spin- 1/2 analog, this 
is equivalent to measuring the operator J^, since 
A^j, = Jl-Jz where /is the identity operator. We assume 
the internal states can be detected with 100 % efficiency 
(Sec. 2.2.1). If all sources of technical noise are elimi- 
nated, the signal-to-noise ratio (for repeated measure- 
ments) is fundamentally limited by the quantum fluctua- 
tions in the number of atoms which are observed to be 
in the |4-) state. These fluctuations can be called quan- 
tum projection noise [100]. Spectroscopy is typically 
performed on L initially nonentangled atoms (for exam- 
ple, ^{t = 0) = nil |i)i). With the application of the 
Ramsey fields, the atoms remain nonentangled. For this 
case, the imprecision in a determination of the fre- 
quency of the transition is limited by projection noise to 
(A(w)nieas = l/(i7RT)"^ where t » Tr is the total averag- 
ing time [100]. If the atoms can be initially prepared in 
particular entangled states, it is possible to achieve 
(A(w)„,eas < \/{LTRTy'^. Initial theoretical investigations 
for ions [3, 9] examined the use of correlated states 
which could achieve (A(w)meas < i/iLTRT)"^ when the 
population (J^) was measured. These states are 
analogous to those previously considered for interferom- 
eters [192, 193]. More recent theoretical investigations 
[194] consider the initial state to be one where, after the 
first Ramsey pulse, the internal state is the maximally 
entangled state 



t// : 



V2 



{\i\\i)2... |i). + e'«''|T),|T)2...|t),) 



(60) 



where (j){t) = (po- Lwot. After applying the Ramsey 
fields, we measure the operator O = Ilf^i S^i instead of 
J^. This gives (A(w)nieas= l/(i^7RT)"^, which is the max- 
imum signal-to-noise ratio possible and corresponds to 
the Heisenberg limit [194]. In the language of quantum 



278 



Volume 103, Number 3, May-June 1998 

Journal of Research of the National Institute of Standards and Technology 



error correction, if we express ^ in terms of the basis 
states |0)'= (|0) + |l))/2'° and |l)'= (|0) - |l))/2"', we 
find that cj)(t) is determined from a parity check of the 
total state in this second basis [142]. For an atomic clock 
where the interrogation time Tr is fixed by other con- 
straints, this means that the time t required to reach a 
certain measurement precision is reduced by a factor of 
L relative to the nonentangled atom case. This improve- 
ment is of significant practical importance since, to 
achieve high measurement precision, atomic clocks are 
run for averaging times t of weeks, months, and even 
longer. 

Cirac and Zoller [1] have outlined a scheme for pro- 
ducing the state in Eq. (60) using quantum logic gates. 
Starting with the state '^(t = 0) = nf=i|4-)i|0), we first 
apply a 77/2 rotation (flo,ot = ttIA, cj)i = +ttI2 in Eq. (23)) 
to ion 1 to create the state "^ = 2'"^(|4-)j + |T)i)|4-)2|4-)3 
. . . |4-)l|0). We then apply the CN gate Ci,,, sequentially 
between ion 1 and ions i = 2 through L to achieve the 
state of Eq. (60). An alternative method for generating 
this state is described in Ref. [194]. As a final example, 
we consider a method for generating the maximally 
entangled state which requires a fixed number of steps, 
independent of the number of ions. For simplicity, we 
illustrate the method for three ions. Starting with the 
state ■^(f = 0) = |4-)|4-)|4-)|0), we first selectively drive 
one of the ions (say ion 3) with a carrier 77/2 pulse 
followed by a red sideband 77 pulse to give the sequence 
|i)|i)|i)|0) ^ 2-- |i)|i)(|i) + |T)) |0) ^ 2-- |i)|i) 
|-i-)(|0)) + 1 1)). We now use the Lamb-Dicke dependence 
of the carrier transition (Eqs. (56)) to make an odd 
(even) integer number of 77 flips correlated with the 
|n = 0) state and an even (odd) integer number of 77 flips 
correlated with the |n = 1) state (the laser beam intensity 
is assumed to be the same on all ions). We can now 
employ a transition to an auxiliary level. For example, if 
the state after the last step, is the state 2'"\ | t)|t)|T)| 1) 
+ |-l)|-i-)|-l)|0)), we could apply a blue sideband 77 pulse 
between states |T) and |aux) of one of the ions (say the 
third) followed by a carrier 77 pulse on this transition to 
carry out the steps 2-"^(|T)|T)|T)|l) + |i)|i)|i)|0)) -> 
2-"^(|T)|T)|aux) + |i)|i)|i))|0)) ^ 2-"^(|T)|T)|t) + 
1 4-) 1 4-) 1 4-)) 1 0). (If the state after the previous step was 
2-"\|t)|T)|T)|0) + 1 4.) 1 4.) 1 4.) 1 1)), we would sandwich 
the last operation between 77 pulses on the |4-) — > |T) 
transition of the selected ion.) 

Other correlated states can also be useful for spec- 
troscopy. A strategy which essentially measures the 
variance of A^j, is discussed by Holland and Burnett 
[195] and Kim et al. [196]. This method has also been 
incorporated into a proposed technique for spectroscopy 
of internal states of Bose-Einstein condensates [197]. 

In comparing the case for entangled vs non-entangled 
states in spectroscopy, the above discussion has assumed 



that Tr is fixed. This constraint would be valid if the ions 
were subject to a constant heating rate and we desired to 
maintain the second-order-Doppler (time dilation) shift 
below a certain value, for example. However, the use of 
entangled states may not be advantageous, given other 
conditions. For example, Huelga et al. [198] assume that 
the ions are subject to a certain dephasing decoherence 
rate (decoherence time less than the total observation 
time). In this case, there is little advantage of using 
maximally entangled states over non-entangled states. 
The basic reason is that the maximally entangled state 
decoheres L times faster than the states of individual 
atoms. Therefore, when using the maximally entangled 
state, Tr must be reduced by a factor of L for optimum 
performance. Because of this, the gain from use of the 
maximally entangled state is offset by the reduced value 
of Tr. (Huelga et al. [198] actually show that a modest 
improvement can be obtained under these conditions by 
use of partially entangled states.) In appendix A, we 
compare entangled vs nonentangled states in the context 
of a practical atomic clock application where a refer- 
ence oscillator is locked to the atomic resonance. 

4. Decoherence 

The atomic motional and internal states, and the 
(logic) operations, were described above in an idealized 
fashion. In this section we consider some of the practical 
limitations to these idealizations. These limitations can 
generally be grouped under the heading of decoherence 
if, by decoherence, we mean any effect which limits the 
fidelity of these operations (see Sec. 4.3). This is a more 
general use of the term decoherence; in some treat- 
ments, decoherence refers only to dephasing of qubit 
states and does not include state changes. Although 
somewhat arbitrary, we also find it convenient to break 
decoherence into categories: (1) decoherence of the ion 
motion, (2) decoherence of the ion internal levels, and 
(3) decoherence caused by nonideal applied fields 
which are responsible for the logic operations. 

4.1 Motional Decoherence 

For the trapped ion system discussed in this paper, 
decoherence may be dominated by that of the motional 
state. Scaling of decoherence will depend on the physi- 
cal system being treated and the mechanism of decoher- 
ence [199, 200]. For quantum computation with ions, 
motional decoherence is somewhat easier to character- 
ize than for a general motional state since we are 
primarily interested in relaxation of the |n = 0) and 
|n = 1) motional states for a particular mode (for exam- 
ple, the center-of-mass (COM) mode along the axis of a 
linear trap). 
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4.1.1 Phase Decoherence Caused by Unstable Trap 
Parameters 

A simple form of motional decoherence is caused by 
fluctuations in trap parameters. Most likely, by employ- 
ing electronic filtering, these parameters fluctuate 
slowly on the time scale of the basic operations 
(— l/i7„>), therefore, the motion is subject to dephasing 
due to the corresponding adiabatic changes in motional 
frequency. In the linear trap, if we assume a, « qf , 
iG{x, y} (Sec. 2.1), then for small fluctuations 8Vo, 
SUo, SOy, Sk , and SR , we have 



where h is the mean number of motional quanta when 
the ion is in equilibrium with the environment 
(n = [exp(^coJkBT) - 1]"'), and yis the relaxation rate of 
the energy to thermal equilibrium. Since ion trap exper- 
iments will typically operate in the situation where 
^Wj « ksT, then h — k^TlAw^. We will assume the 
ion(s) start in the subspace of density matrix elements 
Poo, Poi, Pio, and pii where py = {i\p\i). Equation (62) 
implies 

p^,{t = 0) = Vlyifi + \)pn - r(2n + l/2)po, 



Ox.y/^x,y ■ 



f>VJVo - hnr/nr-iM/R, 



-2nypQi , Poo = - r[MPoo - (n + l)pii]. 



hcojco, = l/2(8/<//< + Wo/Uo) 



(61) 



2yin + l)p22 - r(3n + 1) pn + ^wpoo • (63) 



The relationship between these frequency fluctuations 
and phase fluctuations in a series of logic operations is 
discussed in Sec. 4.3.2. The effects of modulation of 
these parameters by high-frequency noise is considered 
in Sec. 4.1.3 and in Ref. [66]. Although an experimen- 
tally open question, it is expected that all of these 
parameters could be controlled sufficiently well that 
they should not be the primary cause of decoherence. 

4.1.2 Radiative Decoherence 

Decoherence has received considerable attention in 
connection with quantum measurement [201] and has 
been put forth as a practical solution to the quantum 
measurement problem [55,56,57]. In quantum optics, a 
paradigm for decoherence has been to consider relax- 
ation of the harmonic oscillator associated with a single 
mode of the radiation field by coupling to the environ- 
ment [202-206] . This kind of fundamental decoherence 
has recently been observed in the context of cavity QED 
by Brune et al. [207]. An important result from these 
studies is that relaxation of superposition states occurs at 
a rate which increases with the separation of the states 
in Hilbert space and almost always precludes the exis- 
tence of "large" Schrodinger-cat-like states except on 
extremely short time scales. 

A fundamental source of decoherence for the COM 
mode of ion motion is understood by considering that 
the dipole associated with the oscillating charged ion(s) 
is radiatively coupled to the thermal fields of the envi- 
ronment, at temperature T. The master equation describ- 
ing the evolution of the density operator p for the motion 
(in the interaction picture) can be written [208,209,210] 

'V 

p = ^ (n + 1) (2ap^a^ — a^ap — pa^a) 
+ -i n(2a pa — aa p — paa ) , (62) 



General expressions for dpidt and d{h)ldt are given in 
Appendix B. Based on these expressions and in the limit 
that h» 1 , we will characterize the motional decoher- 
ence by the time i* = \l{hy), which is approximately the 
time for the ion to make a transition from the ground 
state. This agrees with a classical estimate [103]. 

4.1.3 Radiative Damping/Heating 

The electric dipole associated with the an ion's COM 
oscillatory motion will couple to thermal (black body) 
or ambient radiation in the environment. However, since 
the wavelength corresponding to ion oscillation frequen- 
cies will typically be much larger than the trap electrode 
spacings, this coupling can be described by lumped-cir- 
cuit models [102, 103]. In these models, we assume the 
ion's motion induces currents in the neighboring elec- 
trodes; these currents, in turn, couple to the resistance of 
the electrodes or circuit elements attached to the elec- 
trodes. In experiments where these resistances are pur- 
posely made high to maximize damping of the ion COM 
motion, observed time constants agree with the model 
[77,78,79,102,103,104]. In the two (single-ion) experi- 
ments which have been able to achieve cooling to the 
n = motional state, the measured value oft* was about 
0.15 s for a "^Hg"^ ion [44] and about 1 ms for a single 
'Be^ ion [17,21,45,47,131,132,211]. In these experi- 
ments, t* was intended to be made as long as possible; 
however, the observed values of t* were considerably 
shorter than what we would calculate from the model, 
as shown below. 

The model considers the electric-circuit equivalent 
shown in Fig. 6(a). Effectively, the electric-dipole 
oscillator formed by the ion COM motion can be consid- 
ered to be confined in a cavity formed by the trap 
electrodes. A useful representation of this situation is to 
model the COM motion (in one direction) as a series 
inductive-capacitive (/c ) circuit which is shunted by the 
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Fig. 6. Schematic diagrams of the lumped circuit equivalents for the trap electrodes and trapped ion(s). (a) The 
left part of the diagram shows, schematically, the electrodes for a Paul trap with hyperbolic electrodes and a 
collection of trapped ions. On the right is shown the corresponding lumped circuit equivalent; Cj represents the 
inter-electrode capacitance (the combined effects of the capacitances shown in the left part of the figure); r 
represents the resistive losses in the electrodes and connecting wires, /l and cl represent the equivalent inductance 
and capacitance for the COM mode of oscillation in the z (vertical) direction [103]. (b) A schematic diagram of 
the endcaps electrode for the trap of Jefferts et al. [211] which was used in the NIST experiments (the ring is 
not shown). Induced currents in the z (vertical) direction are assumed to follow a path indicated by shading; the 
resistance in this path represents r in part (a) of the figure, (c) The rf potential between ring and endcaps electrode 
(or between pairs of rods as indicated in Fig. 1) is typically generated with a resonant rf step-up transformer. The 
resistance in this transformer can, in principle, couple to the ion motion as discussed in the text. 



capacitance of the trap electrodes as shown [9, 103]. 
The resistance r is due to losses in the electrodes and 
conductors which connect the electrodes. The Johnson 
noise associated with this resistance can heat the ions. 
The equivalent inductance of the ion COM motion is 
given by II — m(^IL{aqY where d is the characteristic 
internal dimension of the ion trap electrodes, L the num- 
ber of ions, and a is a geometrical factor on the order of 



1 which can be computed [212]. For traps with hyper- 
bolic electrodes, if we consider motion in the z direc- 
tion, d = 2zg (the separation of the endcap electrodes) 
[67-70] and a =- 0.8 [213,214,215]. For the trap used in 
the NIST experiments, where zo= 130 (xm, /j — 60 000 
henries! The resistance r yields a time constant 
It/r = My. This implies [9] 
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where Q = wji/r is the quahty factor of the ion oscilla- 
tor. The last expression in Eq. (64) shows f * in terms of 
the spectral density of electric field fluctuations at the 
site of the ion which can be written Se = Ak^Tr{ald)^ 
where Ak-^Tr is the Johnson noise voltage associated 
with the effective resistor r. The trap reported in the 
NIST experiments had the endcaps made of a single 
piece of molybdenum as shown schematically in Fig. 
6(b) (ring electrode not shown). We assume the induced 
currents flow in the endcaps electrode as indicated in the 
shaded portion of Fig. 6(b), where Ss is the skin depth. 
This seems to be a conservative estimate since currents 
will also flow in the sides of the endcap electrodes and 
will not be confined to the skin depth [216], thereby 
reducing the effective value of r. Taking the resistivity 
of molybdenum to be p(Mo) = 5.7 X 10"* fi cm, Wt = 
1 25 (xm, and Xt = 1 mm, we find r — Ipx-Y/S^w-r = 0. 04 1 5 
Cl. If we assume T — 300 K, t* — 4.6 s, considerably 
longer than the observed value of f *(— 1 ms). An alter- 
native model for dissipation of charges moving parallel 
to a nearby surface [200] predicts a much larger value 
off*. Lamoreaux [217] has derived an expression which 
agrees with Eq. (64), however he chooses a value of r 
higher than what we calculate. 

Faster heating will occur if T » 300 K. This can be 
expected at the relatively high powers delivered to the 
step-up transformer used to generate Vo (— 1 W in the 
NIST experiments), but this alone cannot explain the 
difference between what the model predicts and the 
observed heating rate. Conversely, if the trap can be 
operated at cryogenic temperatures, this kind of heating 
should be substantially reduced. 

Heating can occur in the axial and radial directions 
due to the interplay of a stray static field (e.g., from 
patch potentials on the electrodes) and noise on Uo or Vo 
near one of the secular frequencies. Here, we explain 
what appears to be one important case, a fluctuation of 
Uo in the presence of a stray static field along the z 
direction of a linear trap. This and other cases are dis- 
cussed more fully in Ref. [66]. 

In equilibrium, the force on an ion from a stray static 
field E = E^z is balanced by the field from the trap given 
by Eq. (2). We have E^ = iKqUoZ^uih where Zequii is the 
equilibrium position of the ion (here, we assume 
Zequii = in the ideal case). A fluctuation in Uo therefore 
causes a fluctuation in the electric field seen by the ion. 
We can characterize the spectral density of these field 
fluctuations as Se{o^) = {EJUo)^Suo{o^) where Suoioa) is 
the spectral density of potential fluctuations. From 
Eq. (64), we have 



For a very small linear trap where k — (0.3 mm) ^, and 
for m = 9 u (e.g., 'Be*) and wJItt = 10 MHz, we have 
t/o == 17 V. For £, == 100 V/m and Suo = {l nV)'/Hz (the 
Johnson noise voltage from a 60 fi resistor at room 
temperature), we have t* — 430 s. Since t* «: w^d^, we 
see there is a premium on having a relatively large trap 
with large values of Uo to keep w^ as large as possible. 

In the above, we assumed that the ions couple to the 
surroundings through the oscillating electric dipole due 
to their COM motion. In situations where the extent of 
the ion sample is small compared to the distance to the 
electrodes, the induced currents result dominantly from 
the COM mode; therefore radiative decoherence from 
modes other than the COM mode can be substantially 
suppressed [1, 103]. For example, for two trapped ions 
aligned along the z axis, we would expect electric fields 
from stray (fluctuating) potentials on one of the end 
electrodes to cause an excitation force on the z stretch 
mode which is suppressed by a factor equal to the ratio 
of the ion spacing to trap dimensions compared to the 
force on the COM mode. 

Fluctuations in Vo and Uo for E^ = can also cause 
heating of the ions. These sources are discussed in Ref. 
[66]. Heating might be caused by parametric processes. 
For example, heating could be induced if the trap pseu- 
dopotential is modulated (coherently or by noise) at 
twice the secular frequency. This problem has been 
treated by Savard et al. [218] in the context of optical 
dipole traps for neutral atoms (a kind of Paul trap for an 
electron to which the atomic core is attached). For the 
conditions of the NIST 'Be* experiments, this kind of 
heating was estimated to be too small to account for the 
observed value of f * [66]. 

4.1.4 Injected Noise 

Noise from various ancillary electronic devices might 
be injected onto the electrodes; this additional electronic 
noise could then heat the ions. Added electronic noise 
can be modeled as a resistor r in Fig. 6 that has a 
temperature much higher than the ambient temperature. 
These sources of noise can be tested by injecting noise 
at a level equal to or above the ambient noise level and 
looking for a shortening of t*. For this test to be valid, 
we must have a reliable means of sensing the noise at the 
trap electrodes. This may be difficult to achieve in prac- 
tice, since, in the experiments, it is usually desirable to 
filter the electrodes from the rest of the environment at 
the motional frequencies. This was the case in the NIST 
experiments, where electronic filtering at the motional 
frequencies precluded the direct observation of voltage 
noise on the electrodes. 
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4.1.5 Motional Excitation From Trap RF Fields 

The rf fields used for trapping in a Paul trap can lead 
to excitation of the ion motion. We will consider four 
types of effects in which the rf micromotion can, indi- 
rectly, cause heating. For the first type of effect consid- 
ered, we will analyze heating of the axial motion of a 
single ion in a conventional spherical-quadrupole Paul 
trap; the results can be generalized to other cases such 
as the heating of radial modes in a linear trap. For a 
spherical quadrupole trap, motion in the z (axial) 
direction has the same form as Eq. (4). If we assume a 
potential yocos(/2r0 is applied between the ring and 
endcap electrodes, we have 



z{t) == A^cos(w^r+ <f),) 



1 + ^ cos(/2tO 



(66) 



where A^ and 0^ are set by initial conditions, 
q^ = ^qVol{mfl^{ri + 2zo)), w^ — q^nj/H"^, ro is the inner 
radius of the ring electrode, 2zo is the distance between 
endcaps, and we assume q^« 1 . In the radial direction, 
the motion will be similar with radial secular frequency 
(Wr = <Wz/2. From this equation, we see that the ion's 
motion in the z direction has components at frequencies 
ft) = /^T ± (o^ Since the rf voltage Vo is typically applied 
through a resonant step-up transformer [shown sche- 
matically in Fig. 6 (c)] the ion's motion at these frequen- 
cies might be expected to couple to the resistance Rt 
between the ring and endcaps associated with this step- 
up transformer. At a frequency w near /2r, the 
impedance between the ring and endcaps electrode can 
be represented by a parallel tuned circuit as shown in 
Fig. 6(c). This impedance is given by 



Z{co) = 



1 + 2iQ 



n. 



/2t 



r,(«) + iXiw) , (67) 



where Q is the quality factor for the circuit 
(Q — Rilifl-iLi) = RtOtCt). Coupling to the effective 
series resistance rs((w) should not occur if the endcap 
electrodes are placed symmetrically around the ring 
electrode (as intended). However, since the relative elec- 
trode positions are difficult to control in small ion traps, 
a displacement of the ring electrode toward one of the 
endcaps will cause a net induced current from ion 
motion to flow between ring and endcaps at frequencies 
o) = Or ± o)^. We characterize this current (in the z 
direction) by / = fSqz (2zo), where z is the ion velocity 
and /3 is a geometrical parameter which expresses the 
coupling to the electrodes (/3 = when the ring is placed 



symmetrically between endcaps). The effective induc- 
tance of the ions for this type of coupling is given by 
IL = IJP^- Associated with r^{(o) in a small bandwidth 
Aj' around w is a series Johnson noise (Vn) = A-kBTAiTs- 
The electric field associated with this noise at frequen- 
cies ft) = i7T ± <Wj can heat the ion motion in a way 
similar to the way in which the motion can be excited by 
a coherent excitation at these frequencies [68]. From 
Eq. (3.7) of Ref [68], we see that an electric field Eiz 
applied at a frequency O^ ± <Wz is equivalent to an 
electric field [coJ{fiT\/2)]E ,z applied at frequency co^. 
Therefore, the Johnson noise from the series resistance 
r^{(o = /^T ± <Wz) is equivalent to that from a series 
resistance rl = [w^ /{2flj)]rs (w = /^t ± w^) at frequency 
(Wj. The heating from this source is characterized by the 
heating time /*' = lll{r'sn). For the NIST single 'Be"^ ion 
experiments, this source of heating was estimated to be 
negligible. 

A second type of rf heating can occur due to the 
Coulomb interaction between ions. In a collection of 
ions, such a string of ions in an ideal linear trap, the 
Coulomb coupling between ions makes all of the 
motional modes, except the COM modes, nonlinear. 
This can lead to excitation of these modes in a Paul trap 
by the driving fields at frequency ili. This excitation 
and the resulting chaotic motion have been studied 
extensively for two ions trapped in a conventional Paul 
trap. Experiments have been performed at Munich and 
IBM; these studies are discussed in Ref. [219]. More- 
over, even for a single (harmonically bound) ion, nonlin- 
ear subharmonic excitation can occur if the exciting 
field is inhomogeneous [220]. Both types of heating can 
be made negligible when the mode frequencies are not 
submultiples of /2r, when a,, qf « 1, and when all 
modes are sufficiently cooled and therefore very linear. 

Another type of rf heating occurs in some experi- 
ments when the condition a,, qf « 1 is not rigorously 
satisfied and the trapped ions are fairly energetic. The 
motion of single ions (or multiple ions when the mutual 
Coulomb interaction or "space charge," can be 
neglected) will be unstable when the condition 
pu)^ + mu)r = fi-Y ^^ satisfied in a spherical quadrupole 
trap or when /jw^: + mw^, = /^t in a linear trap) ip and m 
are integers). This type of heating has been observed in 
some beautiful experiments [221] and has been ex- 
plained theoretically [222]. These "heating resonances" 
arise from terms in the trap potential which are higher 
order than quadratic. We briefly explain their origin for 
a single ion. For simplicity, we neglect the contribution 
of the static potentials; their inclusion will not change 
the analysis significantly. 
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In general, the potential of the trap can be expanded 
about the equilibrium position of the ion and written in 
spherical coordinates (r, 6, cj)) as 



<P=VoCosnTt^ 2 Ci.' 



d 



n„(0,0), (68) 



where F; ,„ are the spherical harmonics and c? is a charac- 
teristic dimension of the trap. We take d = R for the 
linear trap (and d=(ro + Izaf^ for the spherical quadru- 
pole trap) [67,68,69,70]. For the ideal linear trap 
(Eq. (1) with t/r = 0, and neglecting a constant term), 
only two terms in the expansion in Eq. (68) contribute, 
C2,2 = C2,-2 = -(2i7/15)"l (For the ideal spherical 
quadrupole trap, neglecting a constant term, only one 
term contributes, C2,o = 4(i7/5)"^.) For a nonideal linear 
trap, the resonance heating can be explained as due to 
terms in Eq. (68) which give rise to a Hamiltonian of the 
form 



Hi = q<l>=qVp,, 



cos HtI 



(69) 



In the interaction picture for the motion, this becomes 



m 



qVp,, 



(flj)" ia;r + h.i 



(70) 



where this last expression is the leading term which 
satisfies the resonance condition poi^ + mWy = O^. This 
interaction will be suppressed because of the inherent 
smallness of high-order anharmonic terms Vp,i„ (for sim- 
ple trap electrode shapes) and the smallness of the terms 
(xo/Rf and (yo/R)"'. Furthermore, if the x y modes are 
cooled to near the zero-point energy, matrix elements of 
the motional operators will be near 1 . For large ampli- 
tudes of motion, the mode frequencies are not well de- 
fined because of anharmonic terms and heating from 
this coupling would be expected to occur. In any case, 
it is easy to check for a resonant heating of this type by 
varying the resonant frequencies w^ and W3, relative to 
/2r. It can also be checked by varying the initial ampli- 
tudes of motion in the modes. These tests were used for 
the NIST single 'Be* ion experiments; no change in the 
heating was observed. 

Similarly, in a linear trap with L ions, we would 

expect resonances to occur when 2 m-kCOk = fir where m^ 

are integers and w^ are normal mode frequencies (see 
Sec. 2.3.2). If all modes are cooled to the point where 
only a few motional states are excited, then the mode 
absorption spectrum will consist of sharp features 
around the mode frequencies, and the resonances can be 
avoided by changing the trap parameters. Moreover, the 



coupling parameters (Sec. 4.1.8) will, in general, be 
very small. 

A fourth type of heating due to rf trap fields is ex- 
plained as follows. A common problem in ion trap ex- 
periments is the presence of stray static electric fields. 
These fields can give rise to coherent motion at fre- 
quency /2r and potential to heating, which must be ac- 
counted for. Stray static electric fields can arise from 
potential variations on the electrode surfaces ("patch" 
fields) due, for example, to the finite crystalline grain 
size of the electrode material [223], or charge buildup 
on the trap electrodes. Charge buildup can occur be- 
cause, typically, ions are created by electron impact 
ionization of neutral atoms which pass through the trap. 
Often, the ionizing electrons are also collected by the 
electrode surfaces. Electrode charging is particularly 
important at low temperatures where, apparently, ad- 
sorbed gases on the electrodes can provide an insulating 
surface upon which stray charge resides for long periods 
of time (hours). 

If stray static electric fields are present, the equi- 
librium position of an ion is shifted to a place where the 
force from the stray field is counterbalanced by the force 
from the pseudopotential. We will analyze the effects of 
such stray fields using a classical treatment of the 
motion of a single trapped ion. In general, £s„ay = E^^ 
X + Esyy + EsyZ- Stray fields along the z direction in a 
linear trap merely shift the origin along this direction 
and can therefore be neglected. For balance in the x and 
y directions, we have F^ = qE^ - d(q<l>^)/dx = and 
Fy = qEsy - diq(l>p)dy = where 0p is given by Eq. (6). 
This leads to offsets Ax and Ay and resultant motions 
which are, to first order in q^ and qy given by 



x(f) — (Ax + A^ cos(w^f + (p^)) 
Ax = ^, 

yit) — (Ay + Ay cos(w/ + (jty)) 



^y = €5 



1 + ^cosin-rt) 



1 + ^cosiOjt) 



(71) 



where, as in Eq. (4), A^ and Ay are the amplitudes of 
secular motion. The presence of the offsets Ax and Ay 
means that the ion motion has an additional component 
at frequency /2r. This motion will effectively give rise 
to sidebands on the applied radiation as seen by the ion, 
thereby reducing the size of matrix elements between 
states. For example, for single-photon transitions driven 
by a traveling wave with wavevector k = k^x + kyy , the 
electric field from this traveling wave at the site of the 
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ion is proportional to exp(( [^^(0 + kyy(t) - cot + (p]) + 
c.c. Terms in the exponential like k^Ax are (constant) 
phase shifts which can be neglected. Terms like A^. 
cos(wi + cj)^) [1 + (^;c/2)cos(/2rO] are just the motion of 
the ions in the ponderomotive potential in the ideal case. 
The factor from the remaining term, 



exp(( ^ [kjiAx - kyAy] cosCilxO) 



= cos[0fl cos(i7T0] + ism[cj)Q cos(i7 jt)] , (72) 



where cj)o = (q^ /2)[t Ax -kyAy] gives rise to frequency 
modulation sidebands on the spectrum which are spaced 
by /2r. If we consider the carrier or central part of this 
spectrum (the first term in the expansion of 
cos[<^flCos(/2r^], we find that the matrix elements are 
reduced by the factor Jo(0fl) compared to the case 
where the static fields are absent. For cj)o « 1, Jo((t>n) 
— 1 - (cj)Q/2y — exp(- (cJ)q/2)^). Therefore the effect of 
the micromotion looks like an additional Debye- Waller 
factor due to the smearing out of the atom's position 
over the exciting wave (see the discussion following Eq. 
(25) and Sec. 4.3.5). 

To the extent that offset fields are constant, they 
should not cause heating unless the Or sidebands give 
rise to unwanted spectral components that are close to 
transition frequencies of interest. However, the offsets 
Ax and Ay can lead to a problem if the trapping field has 
noise V„ at frequencies /2r ± co„ that is, VoCosi^T^ -^ 
VoCosi^T^ + V„cos(i7T + (Oi)t + y„cos(i7T - (Oi)t. In this 
case (assuming Ay = for simplicity), the ion experi- 
ences noise fields at frequencies Oj ± w, equal to 
E^ = - d(t>ldx = - V„Ax/R^ [Eq. (1)]. From the first part 
of this section, this is equivalent to noise fields of ampli- 
tude V„Axqx/(4R^) applied at frequency w^ [66]. For the 
NIST single 'Be"^ ion experiments, this effect was esti- 
mated to be negligible. Moreover, this source of heating 
was tested for by purposely applying a static field offset 
and seeing if the observed heating rate increased; a null 
result was obtained. Experimentally, it has been possible 
to reduce stray static electric fields by heating the elec- 
trodes [191, 224] or cleaning the electrodes with elec- 
tron bombardment [225]. Alternatively, they can be 
compensated for with the use of correction electrodes 
[226-228]. 



dent or very slowly varying (< 500 Hz) components 
[223]. However, fluctuating patch fields caused by fluc- 
tuating adsorbate coverage has been studied in some 
cases [229, 230]. These studies differ somewhat in the 
low frequency behavior at time scales comparable to 
diffusion times, but at frequencies v > v„ where l/v^ is 
a time constant characteristic of surface diffusion, they 
predict S((l>n, v) °^ V" where S{<i>n,v) is the spectral 
density of rms potential fluctuations <J>„ (in units of 
V^ - Hz"') and a =^ 111 [230]. To estimate the effects 
of time- varying patch potentials on a single trapped ion, 
we assume the ion is sensitive only to the potential on a 
portion of a nearby electrode. We take the area of this 
portion equal to TrOp where a^ is the distance between 
the ion and the nearest part of the electrode surface. The 
effects of these potentials on the ion motion in one 
direction is then estimated by assuming the ion is cen- 
tered between two capacitor plates of area tio^ sepa- 
rated by a distance 2a^. The fluctuating potentials on 
these plates give rise to a fluctuating field at the site of 
the ion which can then excite its motion. 

Patch potential fluctuations can be caused by the 
fluctuations in the surface coverage from adsorbed 
background gas molecules (or atoms). High frequency 
fluctuations appear to be dominated by surface diffu- 
sion rather than adsorption and desorption [230]. The 
number of adsorbed molecules in an area TrOp can be 
approximated by Af(&) — 9(TTa^)/(TTr^) where 9 is the 
fractional coverage and r^ is the radius of the adsorbed 
molecule. For low coverages (9 « 1), the number of 
molecules will fluctuate randomly, AN — N^'^, which 
leads to fluctuations in coverage A9 — 9^'^rja^. A sim- 
ple model for changes in the surface potential due to 
adsorbed molecules is that the molecules are polarized 
by the surface and effectively screen the surface poten- 
tial. We can relate the change in potential of the plate to 
the change in surface coverage by A<J> — kA9 = K9"^rJ 
flp, where k is a proportionality constant. 

If we take S{<i>n, v) = So (constant) for v < v^ and 
5(0„, v) = Soivjvf'^ for v> v„ we have {(p^) = JS((P„, 
v)dv = 3SoVo- Here (<J'„^)'° is taken to be equal to the 
value of Afp estimated in the previous paragraph. The 
cutoff frequency v^ is given by v^ — 1/idiff where tuff can 
be approximated by idiff = ld/(4D) where D is the diffu- 
sion constant and l^ is the diffusion length [231]. Here, 
we take l^ to be the radius of the effective patch (l^ — Op). 
We then find for v > v„ 



4.1.6 Fluctuating Patch Fields 

Electrode patch fields might also vary in time; if the 
spectrum of these variations overlaps the mode frequen- 
cies, this could lead to ion heating. Investigations into 
patch fields have primarily been done for time indepen- 



5(0„,i') = 4 @D' 



2,al 



(73) 



where an extra factor of 2 has been included to account 
for the two capacitor plates which are placed on either 
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side of the ion. To calculate the heating rate from these 
potential fluctuations, we first note that they will pri- 
marily act over a narrow bandwidth associated with the 
ion's motional frequency. In this case, we can represent 
the fluctuations as coming from the Johnson noise of a 
resistor r at temperature T connected between the capac- 
itor plates, that is, S{(l>n, v) = ^-ksTr (assuming the 
capacitive impedance is much greater than r). There- 
fore, we can rewrite Eq. (64) as 



t* 






(74) 



where w^ is the ion oscillation frequency. 

Since the polarizability of molecules and atoms does 
not change dramatically for different species, we will 
estimate k from a measurement of change in surface 
potential for potassium atoms on tungsten. From Fig. 2 
of Schmidt and Gomer [232], we find k — 3 V. For an 
estimate of 6, we extrapolate the data presented in 
Fig. 6.6 of Tompkins [233] for H2 on tungsten and find 
— 0.13 at a partial pressure of 10^* Pa. The diffusion 
constant for H2 on Mo is approximately equal to 10^" 
cm^ • s"' [233]. To make a comparison with the heating 
observed on a single 'Be"^ ion in the NIST experiments, 
we take h — a^ — 130 |jum, ui^IItt =11 MHz, L = \, 
Z, = 6.2X lO'* H, r, = 1 nm (Sec. 4.1.3), and we find v, 
^ 2.4 X 10^' Hz and t* ^ 3000 s. 

This model is very sensitive to the high frequency 
dependence of 5'(<J'„, v) on v, and because of the very 
low value of v^ estimated here, the model should be 
refined. However, we note that if the value of t*{— 1 
ms) for 'Be* observed in the NIST experiments is caused 
by fluctuating potentials on the surfaces of the elec- 
trodes, this would correspond to 5'(<J'„, 11 MHz) 
— (1.3 nV)^/Hz. This should be detectable with a sensi- 
tive amplifier. Therefore, independent of the model, this 
type of noise may be detectable in a straightforward way. 
Conversely, we note that a single trapped ion in the 
experiments considered here will be an extremely 
sensitive detector of potential fluctuations on electrodes 
in vacuum. 

4.1.7 Field Emission 

Field emission from sharp protrusions on the elec- 
trode surfaces can cause ion heating, either from the 
direct electron-ion Coulomb coupling or from associ- 
ated electronic noise on the electrodes. Field emission 
caused by the trap potentials is not unexpected, and field 
emission points have been observed to grow in a number 
of ion trap experiments. For typical values of Vo and trap 
dimensions, the electron transit times from one 



electrode surface to another are much less than \lili so 
that field emission occurs as if the fields are quasi-static. 
Since the onset of field emission varies exponentially 
with the applied voltage between electrodes, it is 
possible to check for field emission by varying the trap 
potentials by small factors, and monitoring the ion heat- 
ing rate. This technique appears to have ruled out field- 
emission heating in the NIST experiments since the 
change in ion heating was much less than a factor of two 
when Vo was reduced by a factor of two. This argument 
assumes the exponential variation of field emission with 
applied voltage; if the emission or some leakage current 
is less sensitive to voltage changes, this test may not be 
valid. If field emission points are formed, it is usually 
straightforward to remove them by momentarily apply- 
ing a large negative potential to the electrode in 
question. The resulting high current is usually sufficient 
to "burn out" the field emission tip. 

4.1.8 Mode Cross-Coupling From Static Electric 
Field Imperfections 

According to the scheme of Cirac and Zoller [1], the 
operations which provide quantum entanglement of the 
internal states of L trapped ions involve the coherent 
manipulation of a single mode k of collective motion. In 
the quantum logic scheme discussed in Sec. 3.3, this 
mode is typically taken to be the COM mode along the 
axis of a linear trap. A potential source of motional 
decoherence is caused by the coupling of this kth mode 
to one or more of the 3L-\ other spectator modes of 
vibration in the trap. If the 3L-\ other modes of oscilla- 
tion are not all laser-cooled to their zero-point energy, 
then energy can be transferred to the A:th mode of inter- 
est. Even when the spectator modes are cooled to the 
zero-point state, they can act as a reservoir for energy 
from the COM mode. This does not lead to heating but 
can cause decoherence. Ideally, the ions are subjected to 
quadratic potentials as in Sec. 2.1. In practice, higher- 
order static potential terms are present; these terms can 
induce a coupling between the modes. Similar couplings 
are induced by the intentionally-applied time-varying 
fields necessary for providing entanglement; these are 
discussed in Sec. 4.4.7 below. 

We will assume that the higher order field gradients 
act as a perturbation to the (harmonic) normal mode 
solution. Following the convention of Eq. (30), these 
fields will be specified by £, for Z £ { 1,2, . . . 3L} where 
the index i specifies both the ion and direction of E . We 
write the electric field at thej'th ion as 



Ej = Ejx + El+j y + E2L+J f , 7 e { 1 ,2, , 



..L}. (75) 
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From Eq. (31), we can write the equation of the A:th 
normal mode as [109, 234] 



dt m j^j 



In general, we can write 

£, = £,({Mp}) =£,({?;}) 



(76) 



1 



Ei{{qj 



3L 3L 



3L 



dEj 
dq-n. 



(9j) = 



+ 9 E E ?/?™ 

^ /=1 m=l 



_d%_ 

dqidq„, 



\ijV 



(77) 



where the derivatives are evaluated at the equilibrium 
positions. The first term on the right side of this equa- 
tion just gives rise to a shift of the equilibrium positions, 
and the second term can be absorbed into new defini- 
tions of the normal mode frequencies w, . The second- 
order term (last term shown in this equation) can reso- 
nantly couple two modes of oscillation (Z and m) to the 
normal mode of interest k. We find a possible resonant 
term: 



d^qk 



3^2" + COkqk - — 2j ^k 



dt' 



■'k qtqm 



dqidq„, 



(9j) = 



(78) 



where the I and m mode frequencies satisfy |&)/ ± u>m\ = 
u)k- This type of coupling can either add to or extract 
energy from mode k, depending on the relative phases 
of motion in the three modes. By substituting the free 
solution to modes I and m [qj(t) = Qjexp(±i(wjt + cfij)] 
into the last equation, we find that if qic(t = 0) = 
(dqic/dt)t=o = 0, the driven solution to the amplitude of 
mode k initially grows linearly with time: 



For a single ion, the three normal modes are just the 
oscillation modes along the x, y, and z directions 
(qi = X , q2 = y , q^ = z; Dj, = 5,jt). The mode frequencies 
were ( co^, w,, co,)/2t7 == (11.2, 18.2, 29.8) MHz, thus 
approximately satisfying the condition Wj^ + w^, = w^ . For 
sake of argument, we assume this resonance condition to 
be exactly satisfied. We consider heating of the x motion 
assuming both the y and z modes are excited. From Eq. 
(79), we find |x(0| = \qtAyAAd^EJdy dz],y = z = o) /i2mco^)\ 
where Ay and A^ are the amplitudes of motion in the y 
and z directions. For simplicity, we neglect the fact that 
energy is exchanged between the three modes; that is, 
we assume the amplitudes of the y and z motion remain 
fixed. In this approximation, if Ay = A^= ^, the time it 
takes to excite the x motion to the same amplitude (^) 
is given by t = \2mQ)Jiq^[d^EJdydz](y = , = o))\- If ^ == 
10 nm (corresponding to {Hy) — {n^) — 1 for the condi- 
tions of the single 'Be"^ ion NIST experiments, the field 
gradient required to drive the x motion to an amplitude 
of 10 nm {{fix) — 1) in the observed time of 1 ms is 
approximately d^E^/dydz = 1000 V/mm^. It is highly 
unlikely the gradient was this large for the NIST exper- 
iments, and, furthermore, the resonance condition was 
only approximately satisfied. Moreover, this source of 
heating was easily tested by varying the initial values of 
Ay and A^ (by varying the Doppler-cooling minimum 
temperature through laser detuning) and studying the 
heating rate of the x motion which had previously been 
cooled to the zero point of motion. No dependence on 
the initial values of A3, and A^ was found. In any case, if 
all modes of motion are initially cooled to the zero-point 
state this source of heating vanishes because the as- 
sumed coupling only provides an exchange of energy 
between modes. This places a premium on cooling all 
modes to as low an energy as possible. Finally, it appears 
that this single-ion example gives a worst case scenario 
since, for large numbers of ions, the force on the gener- 
alized coordinates [right hand side of Eq. (76)] requires 
a high-order field gradient to be nonzero. These gradi- 
ents are highly suppressed in the typical case where 
ion-ion separation is much smaller than the distance 
between the ions and the trap electrodes. 
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We illustrate with an approximate numerical example 
which might have been expected to play a role in the 
heating that was observed in the NIST experiments. 
In those experiments, performed on single 'Be"^ ions, the 
heating that was observed was such that the ion made a 
transition from the n = to m = 1 level in about 1 ms. 



4.1.9 Collisions With Background Gas 

Although trapped ion/quantum logic experiments will 
typically be carried out in a high-vacuum environment 
(P < 10"^ Pa), residual background gas collisions can be 
important. The effects of collisions can be broken up 
into two classes: (i) inelastic collisions, which alter the 
internal state of the trapped ion or even change the 
species of the ion, and (ii) elastic collisions, which only 
add kinetic energy to the ion. Both types of collisions 
will cause decoherence, although heating from elastic 
collisions is expected to be to be the chief concern. 
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The most troublesome inelastic processes are chemi- 
cal reactions and charge exchange. A background gas 
atom or molecule can collide and chemically react with 
the trapped ion, creating a different ionic species which 
is no longer useful. For the reactions to occur, they must 
be energetically favorable (exothermic), and, if the back- 
ground neutral is a molecule, the reactions almost always 
proceed since the internal degrees of freedom of the 
molecule can help satisfy energy and momentum con- 
servation in the reaction. In the ion trap experiments 
discussed here, the ion can spend an appreciable amount 
of time in the (optically) excited state due to laser exci- 
tation; in this case the extra energy due to laser excita- 
tion can make an otherwise endothermic reaction be- 
come exothermic. For example, in experiments on 
laser-cooled Hg* ions [235], when ions were excited to 
the metastable 5d'^6s^ ^Dsq level (approximately 4.4 eV 
above the ground state) they reacted with neutral Hg 
atoms in the background gas to cause loss of the Hg"^ 
ions (presumably due to radiative association causing 
Hg2 dimer formation). As a second example, in experi- 
ments on 'Be* ions [47,76,21 1], the ions were converted 
to BeH* upon collision with an H2 molecule when reso- 
nant light was applied to the ^Sig — > ^P 1/2,3/2 transitions. 
As a final example, the formation of YbH* by a similar 
process has been carefully studied by Sugiyama and 
Yoda [236]. The second form of inelastic collision is 
charge exchange, where a neutral background atom 
gives up an electron and neutralizes the trapped ion. 
Both types of inelastic collision depend critically on the 
particular constituents involved. 

Chemical reactions and charge exchange can occur 
only if the interparticle spacing of the two colliding 
partners approaches atomic dimensions. An upper limit 
on these rates is given by the Langevin rate, for which 
background neutrals penetrate the angular momentum 
barrier and undergo a spiraling-type collision into the 
ion [237]. In these collisions, the electric field from the 
trapped ion polarizes the background neutral (polariz- 
ability a), resulting in an attractive interaction potential 
U{r) = — aq^/i^TTEor'^). Impact parameters less than a 
critical value focnt = (a^V'n'eo/AV^)"'* will result in spiral- 
ing collisions, where fi and v are the reduced mass and 
relative velocity of the pair. (Since the ions are assumed 
to be nearly at rest from laser cooling, v is simply given 
by the velocity of the background constituent.) The ve- 
locity-independent Langevin rate constant 
A:Langevin = (TV = TTfocritV Icads to au overall reaction rate 
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used in the NIST 'Be"- experiments [47,76,211], the 
dominant background gas constituent is usually H2. 
ForH2, we obtain A:Langevin = 1.64 X 10"' cm^ — s~'. At a 
pressure of 10~* Pa and a temperature of 300 K, we 
have yLangevin — 0.004 s~'. Other candidate background 
molecules and atoms have similar polarizabilities and 
values of A:Langevin- Experimentally, we observe lifetimes 
of several hours for trapped 'Be"^ and ^"^Mg"^ ions at 
pressures of around = 10"^ Pa, indicating that at least for 
ground state 'Be"^ and ^"^Mg"^ ions the probability of 
chemical reactions with the background gas constituents 
is small. In a cryogenic ion trap for "'Hg"^, the lifetime 
is many days [74]. 

Background gas can heat the trapped ions by transfer- 
ring energy during an elastic collision. The Langevin 
rate above gives too low an estimate of the rate of colli- 
sions which transfer energy to the ion, since a "heating" 
collision need not penetrate the angular momentum bar- 
rier. A conservative estimate for the heating rate can be 
given from the more general expression for the total 
collision cross section (r^]^stic (in the quasi-classical limit) 
in a CA/r"^ potential [238]. We take C4 = agVSTreo and 
find 



0-elastic = 'ITr(l/3) 
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where n is the density of the background gas. In a 
metal/glass room-temperature apparatus such as was 



If we average over a thermal distribution of background 
H2 velocities, this results in the rate constant feastic = 
< oiiasdcv > = 1.23 X 10'a'"(v)"' where v = {2k^TliJLy'\ 
For H2 at 300 K and a pressure of 10~* Pa, we find 
faasdc — 1-24 X 10"** cm' - s "' and yeiasdc — 0.03 s"'. 
Although each collision on average transfers a large 
amount of energy to the trapped ion, we conclude that at 
typical UHV pressures, such collisions will also be rare. 
Collisional heating can be tested by raising the back- 
ground gas pressure. A simple way to do this is simply 
turn off the vacuum pump; collisions as a source of 
heating in the NIST experiments was eliminated in this 
way. However, doing so does not insure that all partial 
pressures increase by the same factor. Preferably, the 
partial pressure of selected gases should be increased by 
leaking them into the vacuum system and looking for an 
increase in ion heating. 

When ions are first loaded into a trap, their kinetic 
energy is in general comparable to the depth of the trap 
(typically greater than 1 eV). In this regime, elastic 
collisions with the background gas are actually benefi- 
cial, as the background gas can provide a viscous damp- 
ing medium and bring the temperature of the trapped 
ions into thermal equilibrium with the surrounding gas 
[239]. This allows initial laser-cooling to proceed much 
faster. 
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4.1.10 Experimental Studies of Heating 

Some experimental diagnostics for heating have been 
discussed above in Sees. 4.1.4^.1.9. As discussed in 
Sec. 4.1.2, we are primarily interested in determining 
the heating from the |n = 0) or |n= 1) states. In the 
NIST single 'Be* ion experiments, heating could be esti- 
mated by first preparing the ion in the I 'I) 1 0) state, 
waiting a certain delay time, and then measuring the 
strength and ratio of the first red and blue sidebands 
[44,45]. This method gives a simple indication of the 
heating rate, but a more complete method is described 
here. 

More recently, we have determined that depletion of 
the I n = 0) state was dominated by a nearly continuous 
and smooth heating which initially took the atom from 
M = to M = 1 . This was established in experiments 
where the distribution of n level populations was mea- 
sured after a delay time using the technique of Ref. [21] 
and the tomographic technique of Leibfried et al. 
[131,132]. For example, in the experiments of Ref. [21], 
many measurements were repeated for each value of the 
time delay in order to extract the population distribu- 
tions from Pi{t) [Eq. (42)]. In this way the time t* [Eq. 
(64)] was determined to be about 1 ms. In this same 
apparatus, if we make the assumption that the heating 
was caused by a coupling to the environment at 300 K, 
the overall time constant to reach equilibrium is 
Tequil — t'ikBTIAo)^) — 570 s. 

Interestingly, under typical operating conditions, 
these tests of heating are highly insensitive to heating 
caused by collisions. In the previous section, we saw that 
at the typical background pressure of 10 ^ Pa for the 
experiments of the NIST single 'Be* ion experiments, 
collisions would be expected to occur at a maximum rate 
of about 0.03 s~\ as estimated from total elastic scatter- 
ing. Strong heating collisions, given by the Langevin 
rate, occur at a rate of about 0.004 s~'. These rates 
indicate a time constant to reach equilibrium with the 
300 K background gas of between 30 s and 250 s, 
shorter than the 570 s time estimated from the observed 
M = — > n = 1 heating rates (previous paragraph). How- 
ever, since each heating-rate experiment takes only 
about 1 ms, at most, only one experiment in about 
30 000 indicates a background gas collision; these events 
are simply lost in the experimental noise. In effect, the 
technique of measuring t' is able to detect the continu- 
ous, smooth heating occurring in between collisions, 
even though collisions are expected to be a stronger 
overall heating source. This effect was even more pro- 
nounced in the experiments of Ref. [44] . 

Although the source of heating in the NIST experi- 
ments is not understood at the present time, it should be 
possible to tell if the heating is radiative or from some 



other cause by comparing the heating rates of the axial- 
COM and axial-stretch modes of two ions in the trap. 
The heating of the COM mode should be nearly the 
same as for a single ion since it can be excited by a 
spatially uniform (oscillating) field. However, since the 
stretch mode will be excited only by a field gradient, 
radiative heating of this mode by fields emanating from 
the trap electrodes should be significantly less. 

4.1.11 Experimental Studies of Motional 
Decoherence 

It is desirable to have some methods to test for deco- 
herence. Full characterization of decoherence in a mo- 
tional state could be accomplished by reconstruction of 
the density matrix [131,132] coupled with a time delay 
between creation and measurement. This complete char- 
acterization may not be necessary and other methods 
have been used. One possibility is to create an interfer- 
ence signal between two states that depends on decoher- 
ence mechanisms and monitor the contrast of that signal 
in time. This type of measurement was used to charac- 
terize the purity of Schrodinger-cat states for trapped 
atoms [47]. Decoherence of Schrodinger-cat states of 
the electromagnetic field, caused by radiative damping, 
has been monitored by studying a correlation between 
observed states of two successive atoms which probe the 
field [207]. One kind of Schrodinger-cat for motional 
states has the form ( I a) + I — a)) where I a) denotes a 
coherent state. A measurement of the value of the 
Wigner function at the origin W{a = 0) may be suffi- 
cient to characterize decoherence [131,240]. If a is suf- 
ficiently big ((0|q!) = 0), we would expect an initial 
nonzero value for W{0) that is damped toward zero in 
the course of decoherence [240,241]. For quantum logic 
with ions, it appears that decoherence of the COM mo- 
tional states in the submanifold of states 1 0) and 1 1 ) is 
of primary importance. It may therefore be sufficient to 
reconstruct the 2X2 matrix spanned by the 1 0) and I 1 ) 
number states of the COM mode and characterize its 
evolution due to decoherence. One scheme is outlined 
here. 

Suppose we initially create the state ^{t = 0) = 
U)(a,o (0)l0) + a,i(0)|l)). in a particular experi- 
ment, after a time t^ the state becomes ^(^d) = 
I 4-)(2Jo Ci,, (/d)l 0) due to motional decoherence. We 
can then apply two analysis pulses of radiation to this 
state. We first apply a it pulse on the I ~l) 1 1) — > I T) 1 0) 
red sideband transition. We assume this pulse takes time 
t\ {flo,\t\ = 17/2). We follow this with a t7/2 pulse on the 
I 'I) — > IT) carrier transition. Here, we simplify the dis- 
cussion by assuming the time for both of these pulses is 
short enough that decoherence during the pulses can be 
neglected. The wavefunction is now given by 
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where (pi and 02 are the phases (Eq. (23)) of the analysis 
pulses relative to the pulses used to create the state 
^(0). We now measure the probability Pj, of finding the 
atom in the | -l) state and obtain 



Pi{A<P) = 72 



ICi,o(fd)P+ \CiM)\ 



2Re[iCUk)CiMW^^] + E I Ci,,(fd)l 



(83) 



where A<t> = (pi — (j)2. In this expression, we have as- 
sumed all cross terms of the form Ci,,Cj,,,+i for / > 1 
average to zero over many measurements; this will be 
true if decoherence is caused by some process which is 
uncorrected with the creation and analysis pulses. 
We find {Re(Ci,oCi,i)) = '/2{Pi('n)- Pi(0)) and 
{Im (Cl.oCi,,)) = V2(Pi(TT/2) - Pii - tt/2)). Therefore, 
from measurements of Pi for four values of A<p,we 
measure the coherence (C|,oCi,i) (or po,i(fd), the off-di- 
agonal matrix element of the motion after time fd)- If 
necessary, we can find the amplitudes I Cj,,o(fd) I and 
I Cj,,i(fd)l by applying the blue sideband to ^(fd) for a 
time T and reducing Pi{t) as described in Sec. 3.2.1. 

4.2 Internal State Decoherence 

For many years, one of the principal applications of 
the stored-ion technique has been for high-resolution, 
high-accuracy studies of internal state structure of 
atomic ions. This capability can be applied for use in 
atomic clocks (for reviews of efforts from a number of 
laboratories, see Refs. [242] and [243]). High resolution 
and accuracy are obtained because ions can typically be 
stored for long times (many hours or days) with minimal 
perturbations to their internal structure from electric 
and magnetic fields. 

Energy level shifts caused by electric fields (Stark 
shifts) are usually small and, in many cases, magnetic 
field level shifts can be controlled well enough; these 
properties lead to very weak decoherence between inter- 



nal state superpositions. As an example, in Fig. 7, we 
show a spectrum taken of a particular hyperfine transi- 
tion in 'Be* ions [76]. This resonance, obtained with the 
Ramsey method of separated fields [173], has a 
linewidth of less than 0.001 Hz (at a frequency Wq/ 
2tt — 303 MHz) and corresponds to a coherence time 
between the two internal levels of the transition of more 
than 10 minutes. (We could independently establish that 
the noise apparent in Fig. 7 was primarily caused by the 
fluctuations in the oscillator driving the transitions and 
not due to decoherence of the 'Be"^ internal levels.) More 
recently, a slightly narrower resonance was reported by 
Fisk et al. [244]. This resonance was observed on a 
much higher frequency hyperfine transition in '^'Yb"^ 
(=-12.6 GHz) which resulted in a g factor of 1 .5 X 10" 
(Q = transition frequency/linewidth). For levels sepa- 
rated by optical energies, very long coherence times are 
also possible because of the very long radiative lifetimes 
of particular optical levels. So far however, observed 
coherence times have been limited by the linewidth 
of the probing lasers to a few tens of hertz 
[235,245,246,247]. 
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Fig. 7. Ramsey signal of the 2s ^Si,2(Mi = - '/2, M, = + '/2) -> 
2s ^Si/2(— 3/2, + V2) hyperfine transition in 'Be* at a magnetic field of 
0.8194 T). This resonance was obtained using a free precession time 
of 550 s. The data are the result of one sweep (one measurement per 
frequency point). The fluctuations in the data were due to the instabil- 
ity of the reference oscillator used to take the spectrum. These kinds 
of measurements indicated that the coherence time for superpositions 
between the two hyperfine states (which could be used as qubit levels) 
was longer than 10 min [76]. 

These relatively small rates of decoherence indicate 
that quantum states considered for clock transitions are 
also attractive as qubit levels for a quantum computer. 
However, the long coherence times obtained in the spec- 
troscopy experiments were obtained under special con- 
ditions which may not always be compatible with their 
use in quantum logic. Therefore, in this section, we 
consider various sources of internal state decoherence 
and how they might be controlled. 
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4.2.1 Radiative Decay 

A fundamental limit to internal state coherence is 
given by radiative decay. In general, we must consider 
decay from both levels of a two-level system; for sim- 
plicity, and as is often the case, we will assume the lower 
state is stable. 

For electric dipole radiation, the decay rate from up- 
per state 1 2) to lower state 1 1 ) is given by yrad = 



considered quasi-static during the time of a single oper- 
ation Top. A common source of low frequency fluctua- 
tions which would typically fit into this category are 
sinusoidal field fluctuations due to unbalanced currents 
in ac power lines. Therefore, in the spin-V2 analog, we 
assume 
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■ ^WoSAi + P(t)l 
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WI(2|/*|l)|V(377eo^c^)where«o-(£2-£i)Mand/* ^j^^^g ajo0«{T„,r'. For this Hamiltonian, 



denotes the atomic dipole operator. For magnetic dipole 
radiation, we have yrad = <Uol (2| /u, 1 1) I V(3'7T€o.^c^). For 
hyperfine transitions, which decay by magnetic dipole 
radiation, we make the approximation | (2 1 /u, 1 1 ) I — /Ab, 
the Bohr magneton. If we assume (Uo/2'7t =10 GHz, then 
Trad — 10"'^ s~\ clearly small enough for the discussion 
here. For optical electric dipole transitions, we make the 
approximation I (2 1 /u, 1 1 ) | — qao where — ^ is the elec- 
tron charge and Oo is the Bohr radius. If we assume 
Wo/2t7 = 10'^ Hz (A = 300 nm), we find %ad == 7.5 X 10' 
s~\ which is much too fast for our purposes. However, 
as is well known, a number of trapped-ion species have 
first excited optical levels which are metastable. A well- 
studied case is Ba"^ for which the first excited ^Dsg state 
radiatively decays by electric quadrupole radiation. It 
has lifetime measured to be about 35 s [95,248,249, 
250]. Other ions also have metastable optical levels 
which are interesting for the purposes here. Some of 
these are considered by Hughes et al. [63,64] and Steane 
[60]. Lifetime measurements of various ions have been 
compiled and reviewed by Church [251]. 



4.2.2 Magnetic Field Fluctuations 

In the absence of purposely applied electromagnetic 
fields which provide the logic operations, uncontrolled 
fluctuating external magnetic fields are expected to give 
the primary contribution to internal state decoherence. 
Decoherence results from the fact that the energy sepa- 
ration between two levels of interest depends on the 
external magnetic field B. In most cases, it will be 
possible to express the transition frequency between two 
levels to sufficient accuracy by 



COo + 8(U = C0o + 
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dB^ 
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where Bq is the average magnetic field and Wo is the 
transition frequency for B = Bq. 

The effects of a fluctuating magnetic field depend on 
the spectrum of the fluctuations. We first assume the 
typical case that the fluctuations are slow enough to be 



Schrodinger's equation yields 

Ci(t - to) = Ci(to)d^'\ 
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In a sequence of logic operations, (^(0 must be small 
enough or be taken into account (see also Sec. 4.3.2). 
As a second case, we assume the magnetic field varies 
rapidly; specifically, we assume B varies sinusoidally: 
B — Bo = ABcoscoBt where cob » 1/Top. We find 



8(u = 



dw 
~dB 



ABcoscoBt 
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dB' 
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(87) 



The term [d^(jj/dB^]ABV4 on the right side of this ex- 
pression is a frequency shift which can be absorbed into 
the definition of coo if AB remains constant. For simplic- 
ity, assume that one of the cosine terms in Eq. (87) 
dominates so that 8(u sinusoidally oscillates at &)„! = <Wb 
or 2cob so that I3{t) = ^ocoswn,? in Eq. (85). The effects 
of this fast modulation can be seen if we consider apply- 
ing external radiation near the carrier frequency coo of 
the internal state transition. Schrodinger's equation 
leads to expressions similar to those of Eqs. (17) which 
take the form (for A = S=0) 
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where i]^ = Wo^lw^. We have 
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where Ji(X) is the (th Bessel function with argument X. 
For (Wm » 1/Top — /2„,„, the sinusoidally varying terms 
on the right average to zero (as in the rotating wave 
approximation). Therefore the wave function evolves just 
as in Eq. (25) except we must replace I /2„,„ | by 
I ^n,ii \Joif]m)- Since the Jg factor can be absorbed into 
the definition of O, this steady modulation should not 
cause a problem. An important source of high-fre- 
quency magnetic fields is from the currents in the elec- 
trodes which oscillate at the trap drive frequency fir. 

The gate demonstrated in Ref. [17] had a strong sen- 
sitivity to magnetic field fluctuations since the qubit 
frequency and the auxiliary transition frequency had a 
dependence dw/dB — ixb/^ — 10'" Hz/T. If this scheme 
is used in future experiments, sufficient magnetic 
shielding must be provided. Alternatively, the magnetic 
field dependence can be minimized by operating at a 
magnetic field where the transition frequency becomes 
independent of magnetic field to first order (making the 
d(o/dB terms in Eq. (84) and (87) vanish). For example, 
most atomic clocks are based on hyperfine transitions in 
ground electronic ^Sig states, where the nuclear spin /is 
half integral (Vz, 3/2, 5/2, ...). In these cases, where the 
total angular momentum (F = 1 + S) is integer, the 
(F,mF = 0) ■'-^ (F',mF = 0) transition is first-order field 
independent for \b\ — 0. These particular transitions 
may not be useful in the applications discussed here 
since other transition frequencies (for example, 
(F,mF = 0) ■'-^ (F',mF= ± 1) transitions) will be very 
close to that of the (F, m/? = 0) •<->^ (F',mF = 0) transition 
frequency as \b\ — > 0. This may cause unwanted cou- 
plings to these other levels. At optical frequencies, sim- 
ilar first-order field independent transitions in ions have 
been used in frequency standards [235]. However, other 
hyperfine and optical transitions become first-order 
field independent at nonzero magnetic field, and their 
frequencies are separated from competing transitions 
(see, for example, Refs. [76] and [189]). 

For quantum logic, storing qubit information in two 
states whose energy separation is first-order field inde- 
pendent therefore appears to be attractive. However, the 
conditional dynamics which is central to quantum logic 
may necessitate a transition to an auxiliary state. In the 
notation of Sec. 3.3, a transition between the I T) and 
I aux) levels will have a transition frequency which will, 
in general, not be first-order field independent if the 
I 4-) — > IT) transition frequency is field-independent. 
Therefore, if the external field fluctuates, the fidelity of 
this operation will be compromised. As an example, 
suppose quantum logic is performed on a collection of 
'Be* ions using the 2s ^Sjc \mi= — 3/2,mj= 1/ 
2) = "\ iy and I - 1/2,1/2) = " I T)" levels as qubits 
(strong-field state representation). With these designa- 
tions, level I -l) has a higher energy than level I T). These 



are the same levels as those used for the clock reported 
in Refs. [75] and [76]; the transition frequency Wq/ 
2'iT — 303 MHz is first-order field independent at a field 
of Bo — 0.8194 T. At this magnetic field, the transition 
frequency has a second-order dependence of Swo/ 
2'IT == - 5.2(8B/Bo)' MHz. Therefore, a field fluctua- 
tion of SB = lO""* T from Bo leads to a frequency offset 
of the transition of only about 80 mHz. For the first CN 
gate described in Sec. 3.3, we could use 2'iT transitions 
between the |T) and "|aux)"= 1 1/2,1/2) level. How- 
ever, this transition has a first-order dependence on 
magnetic field at 0.8194 T given by (c)(w/2'it)/ 
as == - 22 MHz/T. Therefore if the field shifts by 
SB = 10""* T, this will cause a shift of this transition 
frequency of A/2'n = — 2200 Hz. Similar consider- 
ations must be applied to quantum logic using other ions 
such as ^"'Hg* ions where the splittings between hyper- 
fine transitions are much higher [252]. Since it will be 
difficult to find two field-independent transitions at the 
same field, it may be advantageous to use logic gates 
which require only one internal state transition as in the 
second and third gates described in Sec. 3.3. 

4.2.3 Electric Field Fluctuations 

In Sec. 4.4, we treat the effects of electric fields from 
the applied laser beams; in this section we treat the 
effects of other electric fields. Electric field shifts are 
likely to be less important as a source of internal state 
decoherence than shifts due to external magnetic fields. 
First, electric-field shifts are second-order in the field, 
except for exceptional cases involving nearly-degenerate 
states of opposite parity. Second, the fields are relatively 
small at the site of the ion. In analogy to Eq. (84), we 
have 



COo + ?>(Oe = (Oo + 



dho 



(90) 



where E is the magnitude of the electric field. Unlike 
Eq. (84), the expansion is made around E = 0, rather 
than a nonzero value Eo, since, as noted previously, the 
mean value of E is zero for an ion in a trap. If it were not 
zero, then the ion would have a mean acceleration and 
would not be spatially localized. Of course, the ion must 
experience an average electric field to counteract the 
force of gravity, but this field is extremely small. For an 
ideal rf trap, such as one in which the rf potential is 
described by Eq. (1), and in which there are no stray 
static electric fields, the ions are attracted to the line 
along which the rf electric field is zero. In this case, the 
mean-squared electric fields due to confinement experi- 
enced by cooled ions are very small. If stray fields are 
present, the mean-squared electric field seen by the ions 
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will be nonzero and can be significantly larger [228]. 
The stray fields can be reduced by applying potentials to 
compensation electrodes; to some extent this is neces- 
sary anyway for the efficient laser cooling. It should be 
possible to reduce the rms electric fields from the elec- 
trodes and neighboring ions to less than 1 V/cm. We 
note that rms electric fields of about 8 V/cm are present 
at room temperature, due to the blackbody (thermal) 
electromagnetic field. Since the spectrum of this field is 
dominated by infrared frequencies and because it is 
steady, this should lead only to a small ac Stark shift 
[253]. 

We briefly consider two cases that have been dis- 
cussed for ion trap quantum computation: narrow opti- 
cal transitions between the ground and a metastable state 
and transitions between ground-state hyperfine sublev- 
els. For an optical transition between two energy levels, 
the energy shift, derived from perturbation theory, is 
roughly (qaofE^/AE, where qao is the product of the 
proton charge and the Bohr radius, that is, a typical 
electric-dipole matrix element, AE is a typical energy 
difference between one of the two levels and another 
which is connected by the electric-dipole operator, and 
E is the electric field. For the case of the Ba"^ 6s ^Si^-to- 
5d ^DsG transition which has a wavelength of 1 .76 [jum, 
this expression predicts a shift 8&)£/2'n- of approximately 
3 mHz for £ = 1 V/cm, if AE is taken to be the 6s 
^S\/2-to-6p ^Pi/2 energy difference. This is in rough 
agreement with the experiment of Yu et al. [254]. Simi- 
lar values of 8&)£ are estimated for the 5d'^%s ^Sic -^ 
5d%s^ 'Dsc transition in '"Hg"' [228]. Transition fre- 
quencies between hyperfine levels, such as those used in 
the NIST single 'Be* ion experiments, are even less 
sensitive to electric fields, since, to a good approxima- 
tion, an electric field shifts all sublevels of the hyperfine 
multiplet by the same amount. A differential shift be- 
tween the sublevels arises in third-order perturbation 
theory [253]. For 'Be"^, for example, the fractional shift 
between hyperfine sublevels is approximately 
- 4 X 10"" for £ = 1 V/cm [255]. 

So far, we have considered only electric dipole shifts, 
that is, the shifts that are due to a uniform electric field. 
An electric quadrupole shift, which is proportional to 
the product of the electric quadrupole moment of the 
state and the applied electric field gradient, may also be 
present. The shift Swg is given approximately by 



„ qQ^ dEi 

fi drj 



(91) 



where Q^ is the atomic quadrupole moment and dEildVj 
is a typical component of the electric field gradient 
tensor, such as BE^Idz. The precise value depends on the 
whole electric field gradient tensor and the details of the 
atomic state. There is no quadrupole shift for S-states or 



for other states with electronic angular momentum 
/< 1, like the ^Po states in the Group IIIA ions, such as 
In* [256]. The trap rf electric fields have gradients, but, 
since they oscillate sinusoidally, they average to zero. 
Static electric field gradients due, for example, to patch 
effects, or the static potential to confine ions along the 
axis of a linear trap [Eq. (2)] would be more of a prob- 
lem. For a typical value of Qa = 10 ai, and an electric 
field gradient of 10 V/cm', '^WqIItt = 0.7 Hz. Fluctua- 
tions in this (small) shift will probably be small com- 
pared to those of magnetic field shifts. Electric field 
gradients due to neighboring trapped ions will also be 
present [257]. Consider an ion at position x = y = z = 0. 
The electric field gradient component BEJdz due to ions 
having charge q, located atx = 3? = 0, z=±fifis — ql 
(TTeod^). For d= 10 [jum, this is 575 V/cm'. However, this 
should be relatively constant, since the spacing between 
ions is determined by the trap fields, which must be kept 
nearly constant if motional decoherence is to be avoided. 
Thus, they may lead to energy shifts, but not to decoher- 
ence. 

Garg [58,59] has considered the decoherence of inter- 
nal states of ions due to fluctuating electric fields in- 
duced by the reservoir of vibrational modes of a string 
of trapped ions. Such decoherence is present even if the 
modes are in the ground state, because of zero-point 
motion. This form of intrinsic decoherence was found to 
be negligible, at least for the specific case in which the 
6s 'S 1/2 and the 5d 'D5/2 states of Ba* are used to store the 
qubits. For a string of 1000 ions, the decoherence time 
was found to be around lO'' — 10^ times longer than the 
natural lifetime of the upper state, which is 35 s. 

4.3 Logic Operation Fidelity and Rotation Angle 
Errors 

The third and perhaps most important class of deco- 
herence involves imperfections in the logic operations. 
Ideally, a quantum computation transforms an initial- 
ized state vector ^(0) to t/ideai^(0), where 
t/ideai = exp[— (T/ideai^/^] (wlth //ideal tlmc independent) 
is the unitary evolution operator without errors. In prin- 
ciple, both [/ideal and the corresponding interaction 
Hamiltonian //ideal operate on the Hilbert space formed 
by the I -l) and I T) states of the L qubits. In practice, for 
the case of trapped ions, the physical implementation of 
logic gates requires that we must include at least the 1 0) 
and 1 1 ) states of the chosen motional mode and perhaps 
an internal auxiliary state (Sec. 3.3). If decoherence 
mechanisms cause other states to be populated, the 
Hilbert space must be expanded. Although more 
streamlined algorithms may be available, the operator 
[/ideal can always be broken up into fundamental single 
and two-bit quantum gates [140,155,156,157]; for sim- 
plicity, this will be assumed in what follows. 



293 



Volume 103, Number 3, May-June 1998 

Journal of Research of the National Institute of Standards and Technology 



In the previous sections, we have concentrated on the 
decoherence of the state vector ^ due to environmental 
coupling, in the absence of purposely applied additional 
fields. Here, we treat decoherence due to imperfections 
in the evolution operator [/jdeai which is implemented by 
these additional applied fields. We break this decoher- 
ence into two categories: (i) errors or noise in the (clas- 
sical) gate parameters (that is, the rotation angles (p and 
i7„',„f in Eq. (23)) which result in undesired evolution 
and, (ii) coupling to the environment induced by the 
gates (for example, spontaneous emission induced by the 
applied fields). In some cases, we can model noisy logic 
gates by considering an ensemble of gates whose 
parameters are statistically distributed. This allows a 
simple characterization of several candidate noise 
sources. 

The final step in an ideal quantum computation is a 
probabilistic measurement of all the quantum bits, 
which yields a given L -bit number A:, or equivalently, the 
state k, with probability Pic= \{k\Uii^^^(0))\^ 
(A:= 0,1,2,...,2^ — 1). A particular quantum computation 
is thus completely characterized by the vector of proba- 
bility amplitudes Ak = {k\ t/ideai^(0)). The usefulness of 
a quantum computer algorithm relies on the coherent 
interference of the probability amplitudes, resulting in 
only a small fraction of the 2^ possible numbers with 
appreciable measurement probability. If the actual evo- 
lution operator U contains imperfections, the probabil- 
ity amplitude vector is instead Bk = {k\ [/^(O)) 7^ At. 
We characterize the fidelity of a computation by the 
expression 

F=(^\ (t/id.ai^(0) I UnO)) I ') = ( I E A;b, I ' y (92) 

where the average is performed over any random vari- 
ables affecting the operator U. Roughly speaking, F is 
the probability that errors in the operator U do not affect 
the result of a computation, with F= 1 representing 
error-free computation. This characterization of the fi- 
delity follows the approach of quantum trajectories 
[258,259]. A quantum computation can be interpreted as 
the evolution of a complicated path through Hilbert 
space. If the operations comprising the computation are 
imperfect, the path of computation bifurcates into a 
superpositions of the correct path and incorrect path 
[260]. In this model, the fidelity F is just the probability 
that the correct path has been followed. Schemes for 
complete characterization of a quantum process in terms 
of a input-output matrix have been proposed in Refs. 
[261] and [262]. This approach would be able to charac- 
terize the fidelity of any of the operations discussed 
here. It can also be used to test for decoherence of the 
internal and motional states, by use of time delays be- 



tween state preparation and application of the diagnos- 
tic. 

An adequate characterization of the effects of noise in 
a quantum computation requires knowing the type of 
noise processes present, and application of this noise to 
a specific computational algorithm on a specific input 
state. This is most straightforwardly accomplished using 
numerical simulations [1,263-266]. Here, we do not 
treat any specific algorithms, but will attempt to make 
some general observations on particular sources of logic 
operation errors and how these errors accumulate after 
many operations. 

4.3.1 Accumulated Errors 

Each successive operation in a quantum computation 
based on the Cirac/Zoller trapped-ion scheme can be 
regarded as a rotation between two states | A )j | « ) and 
|A')j|«') in a multidimensional Hilbert space where 
\A)j and \A')j are from | l-)j, \ 1)j, or I aux)^ and \n) and 
In') are from 1 0) or 1 1) of the selected motional mode. 
Before the A:th operation, assumed to be applied to the 
jth ion, we can, to a good approximation, write the 
wavefunction as 

^t-l = at-i I l/'t-l,„)j I ijjic-J.a) + Pk-\ I ^k-\.fi)j I lAt-L/s) 



+ yt-il ^k-\,j)j\ ^k-\,j). 



(93) 



where I i/'t-i,a)j and I ^k-\s)i ^^ the basis states for the 
rotation (for example, the I T)^ 1 1) and I aux)^ 1 0) states 
required to carry out step lb of Sec. 3.3). The state 
I ijjt-\,y)j is the part of the wavefunction for the jth ion 
which does not include states involved in the rotation 
(for example, the | i); state). The states I i/'t-i,a), 
I 4'k-i.^), and I i/'t-i,y) are states which include other ions 
and, because of gate imperfections, may involve states 
outside the Hilbert space spanned by the 1 4-), I T), 
and I aux) internal states, and |0) and 1 1) motional 
states . We choose ( i/'t- 1 ,„ I i/'t- 1 ,„ ) = ( i/'t- 1 ,,? I i/'t- 1 ,/? ) = 
(lAt-i.y I lAt-Ly) = 1 so that I a^-, I ^ + I /3t-i h -H I yt-i I ^ = 



1. From Eq. (23), the "actual'' 
written as 

Rka(Ok+^k,(t>k) = 



cos(ft+^t) -;e*sin(0t+^t)" 
-;e-'*'sin(ft+^t) cos(ft+^t) 



A:th operation can be 





"t-i 


, <A = 


>-,_ 



(94) 



where, to simplify notation, we have chosen 
6k = {fin\n t)k and ^t, 0t « tt are the errors in the ampli- 
tude and phase of the rotation. We can write 
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Rka = Rk+ ^Rk, 



(95) 



where Ri, = RtoCos^t^lRto = Rk{9k, 0)] is the ideal A:th op- 
erator and 



^Rk 



sinftsin^i 
(e"''^' - l)sinftcos^i + e"'*cosftsin^i 



(■ (e'*- l)sinftcos^i + e'*cosftsin^i 
sinftsin^i 



(96) 



A/^j. has terms which are first order (and real) in ^<. and 
(^fc. Nevertheless, we might expect that the overall fi- 
delity has errors which depend only quadratically on ^t 
and (f)k since a general computation may approximate a 
sequence of nearly random rotations on a multidimen- 
sional Hilbert space. To see how this might come about 
we express a general computation consisting of M oper- 
ations as 



O = RhlaRl 



Mai^iM- l)a ■ 



.Rka...R2.Rla^Oi+02+0,, (97) 



where 



Ol =RmRm-I-R2Ri, 

O2 ^ [Rm-.RiARi] + [«M...A«2«,]... 
+ ■■■[Rm^Rm-i---Ri] + [^RmRm-i---Ri], 

03 = Rm-Rs^Ri^Ri + Rm-AR^Ri^Ri- 



(98) 



+ ...ARmRm-AR 



Ml^M-l'-^I^M-2- 



.R, + ARmARm-iR 



M'-il\M-ll^M-2---l^l 



R, 



and where we have neglected terms which are higher 
than second order in ^^ and <f>k. In this notation, the ideal 
computation is given by the operator 



Oideal — R 



R2aR\i 



M0---"-20"-[0 



(99) 



For the fidelity, we need to calculate 

(^(0) I O^.ealO I ^(0)) = (^(0) I 0^,eal(0, + O2 + 

O3) I ^(0)) where we will neglect terms of third order 
and higher in { ARj } . For the first term, we find 

M 2 

<^(0) I OiLiO, I ^(0)) - 1 - E T ■ (100) 

For the second term, we have a sum of M terms; the A:th 
of which is given by 

(^(0)1 OiLl«M«M i...A«,...«2«, I ^(0)) 

- {nO)\R}oRL-RLAR,R,-,...R,\nO)) (101) 

= {%-JrI>ar,\%-:,o} 



where ^^-lo corresponds to the wavefunction repre- 
sented in Eq. (93) for the ideal case. We find 

{,!,,.. JrLAR, I (A,- 1,0) - - ^ sin^ft[| «,_, h+ I ft. , I '] 



Pk-i I ^)(^tCosft + sinft)(^iSinft 

(102) 



+ 2Re[ak-,l3k-,](^ sinftcosft - ^k) 



- 2/m[Q;i-i/3t-i](sinft + ^iCosft)</)iCosft]. 

The third term (^(0)1 O^ideaiOsI ^(0)) is straightfor- 
ward to evaluate and is clearly second order in ^^ and 
4>k- Even though the imaginary term in Eq. (102) has 
terms linear in ^j. and <p^ , when the fidelity F is calcu- 
lated, these terms add only in quadrature so the A:th 
operation gives only a contribution of second order in ^^ 
and <f>i,. Therefore, in general, we have 



1-E 



Ck[^k + Ck2^k(Pk + Ck3(Pk 



(103) 



where the coefficients Qi, Q2, C^ — 1- For general 
computations involving many operations M in a large 
Hilbert space, it might be expected that the errors in 
rotations are distributed in an approximately random 
fashion. In this case, we would expect the fidelity to be 
approximately given by F — 1 — M{F^i^^ + F^^<p^ + 
Fi^^<f>) where ^and <p are characteristic of the errors for 
all operations and F^, F^^, F^^ ^ 1. In the case ampli- 
tude errors (^) dominate, then the maximum number of 
operations before the fidelity drops appreciably below 1 
is M^ax — 1/^^. For example, laser intensity fluctuations 
on the order of 1 part in 10"* would give Mmax — 10^ At 
the other extreme, the errors might add coherently in 
some cases. These operations could conspire to give 
contributions to the fidelity that goes as (X^^tf. We 
can illustrate this with a simple example. Assume that 
an intended rotation ft has an error ^t « it. The state of 
the system before application of this rotation can be 
expressed by Eq. (93). From Eq. (92) we find F = 
cos'(^,/2)( I a,- , h+ 1 ft-i I '■)+A[Re { a;-,ft-,exp(;(^,)}]' 
sin^(^<./2). The existence of the second term only serves 
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to improve the fidelity, so at worst (when O/ = 1 or 
jS, = 1), the fidelity of this single operation is given by 
cos^(^j:/2) — 1 — (^t)V8. Now suppose that a second ro- 
tation 6m about the same axis has error ^„. A similar 
analysis shows that when these two operations are ap- 
plied in succession, a worst case for the fidelity is given 
by F=cos{^t/2 + ^„/2) =- 1 - (^, + ^„)78. Similarly, 
for rotations on M bits, the worst case fidelity is given 
hy F = cosC^f ^,/2) so that the errors could accumulate 
linearly. Most likely, this case will correspond to rather 
uninteresting computations and one can generally expect 
that the errors will accumulate in a way which is be- 
tween random walk and linear. 

The effects of operation errors (^,, 0,) have been 
observed in simulations of quantum computations to 
factorize small numbers [263,266]. If the errors fluctu- 
ate randomly about the correct value, the fidelity is 
given by a simple expression which is essentially the 
same as Eq. (103) for f — 1 [266]. However, for con- 
stant errors the fidelity indicated by the simulations is 
somewhat worse indicating that the computational al- 
gorithms cannot be assumed to be rotations about axes 
chosen pseudo randomly. Therefore, the simulations in- 
dicate that requirements on systematic offsets are more 
stringent than the requirements on random fluctuations 
from operation to operation. 

4.3.2 Pulse Area and Phase Fluctuations 

As in the last section, we will consider two sources of 
noise: (1) those which result in offsets from, and fluctu- 
ations in, the ideal Rabi frequencies [for example, off- 
sets and fluctuations of i7„',„f in Eq. (21)] and (2) those 
caused by phase fluctuations between applications of 
separate operations to a particular ion [fluctuations from 
some constant value of (p in Eq. (21)]. 

We first make some remarks regarding the coherent 
evolution between two states of the system which com- 
prise the basis states for a single operation. First, there 
are a number of methods of coherently transferring pop- 
ulation between a two level system with radiation, in- 
cluding constant amplitude Rabi pulses as represented 
by Eqs. (21) and (23), adiabatic rapid passage [267], and 
adiabatic transfer via dark states [268-270] . The advan- 
tage of the last two techniques is that for it pulses, their 
efficiency does not depend critically upon the pulse area 
(Eq. (104) below) of the applied field. However, for 
rotations other than it, this advantage is lost. Moreover, 
for these techniques to work with high efficiency, par- 
ticularly simple atomic level structures are required and 
coupling to any off-resonant levels can potentially con- 
taminate the transfer. For example, adiabatic transfer 
can be extremely efficient between ground state hyper- 
fine levels if the transfer occurs only through a single 



excited (electronic) state. However, in practice, other 
excited states are driven by the coupling (laser fields) 
thereby affecting the fidelity of the transfer [271]. Fur- 
thermore, the high powers required for adiabatic rapid 
passage and adiabatic transfer aggravate decoherence 
from spontaneous Raman scattering (Sec. 4.4.6.4). 
Therefore, we will only consider the case of Rabi 
pulses, where population is transferred by applying a 
resonant field with a well-controlled amplitude envelope 
and duration. 

When a resonant radiation field is applied to a two- 
level system, the resultant unitary transformation is 
given by Eq. (23). We can further generalize this expres- 
sion replacing fl„;„t by 0„_„ (f)/2 where 



t 
6)„.„(f) = 2 /2„.,„(f')df' ("pulse area"). 



(104) 



The factor of two in this equation is introduced so that 
the condition ffl„;„ (f ')df ' = 17/2 corresponds to a it pulse 
of the effective Bloch vector. Here, as opposed to the 
situation assumed in Sec. 2.3.1 (where the coupling 
radiation was assumed to be turned on instantaneously, 
left at constant amplitude, and then turned off instanta- 
neously), we now allow the coupling radiation to be 
turned on and off gradually. This more closely approxi- 
mates a real experimental situation and can be advanta- 
geous since the spectrum of i7„_„(f) can be tailored to 
suppress off-resonant transitions (Sec. 4.4.6). Here, we 
call 6)„',„ the pulse area of the applied radiation and 
2i7„_„(f) the Rabi frequency envelope. If the applied 
pulse area is noisy or set incorrectly, then the output 
state will fluctuate from it's ideal value as discussed in 
the previous section. 

We must also consider detunings between the applied 
frequencies and resonant frequencies of the ions [in Eq. 
(21)] and phase fluctuations between successive opera- 
tions on the same ion [fluctuations in (p in Eq. (21) from 
operation to operation]. This could occur from fre- 
quency or phase fluctuations in the applied fields or 
fluctuations in the ion internal or motional frequencies. 
For frequency fluctuations, we consider that for a given 
operation on ion 7, Aj{t) is slowly varying over the time 
of a single operation (Sec. 4.2.2). In the limit that detun- 
ings are small, the dominant effect can be characterized 
by assuming that A = for each operation, but that 
different phases 



(105) 



0,= AjiOdt' 



are present for successive applications of the same oper- 
ation. In this expression, to is the time that the first 
operation is applied at ion j and 4 is the time of the A:th 
operation (applied to thejth ion). 
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4.4 Sources of Induced Decoherence 

Below, we consider various sources of fluctuations 
and decoherence caused by the logic operations and how 
these might be evaluated and controlled in experiments. 
First, we are interested in controlling the accuracy and 
stability of the Rabi pulse are (9„',„ given by Eq. (104) 
and the accumulated phase as expressed in Eq. (105). 
Since gates will most likely be implemented with laser 
transitions we will examine laser intensity and timing 
fluctuations. 



4.4.1 Applied Field Amplitude and Timing 
Fluctuations 

Fluctuations in the laser intensity at the site of a given 
ion can arise from both fluctuations in the relative posi- 
tion of the beam with respect to the ion and fluctuations 
in laser power. Laser/ion position stability is particularly 
important since the Cirac/Zoller scheme of quantum 
logic assumes that ions in an array be selectively ad- 
dressed, thereby requiring a high degree of control of 
the laser beam spatial profile (Sec. 4.4.4). Of course, the 
simplest method for minimizing position fluctuations is 
to employ mechanically stiff mounts for the optics and 
ion trap electrodes, and have the laser source as close as 
possible to the ions. A quadrant detector indexed to the 
trap electrodes and placed near the ion may also be used 
to actively stabilize the beam position by feeding back to 
a galvanometer or acousto-optic modulator. If optical 
fibers are used to deliver laser beams to the ions, posi- 
tion fluctuations between the fiber and the ions could be 
made small; however, we must also consider position 
fluctuations between the laser source and the input to 
the fiber. If the position tolerances can be adequately 
controlled, the dominant source of intensity noise at the 
ion would likely be given by fluctuations in optical 
power and laser mode. Here, we estimate limits on laser 
amplitude noise. 

If we assume the laser fields responsible for quantum 
logic operations are coherent states, the fundamental 
noise floor is photon shot noise. For a laser beam of 
average power Po, the fractional level of shot noise is 



8P 
Po' 



fid 



PoTox 



(106) 



where w is the (optical) photon frequency. Top is the time 
the radiation is applied, and, for simplicity, we assume 
square pulse envelopes as in Eq. (23). Of course, lasers 
seldom produce amplitude noise at the photon shot- 
noise limit. Almost all laser sources have significant 
amplitude noise in the 10 Hz to 10 kHz range due to 
acoustic vibrations which, for example, affect the laser 
cavity resonators. Much of this noise can be removed by 



active power stabilization, where a beamsplitter directs 
a portion of the laser power to a photodetector, and an 
error signal is derived and fed back to an upstream 
modulator or, in the case of a diode laser, directly to the 
amplitude of the laser source [272,273]. The limiting 
noise of this stabilization scheme is degraded slightly by 
the imperfect quantum efficiency of the photodetector 
as well as the beamsplitter. If the beamsplitter directs a 
fraction e of the input optical power to the stabilizer 
(which then gives an optical power P^ — {\ — e)Po di- 
rected to the ion), and the quantum efficiency of the 
photo detector is T7det, the limit of fractional power noise 
in the logic pulse is (assuming no added electronic noise 
in the feedback loop) 



Pn 



fid 



/'uTop77dete(l - e) ' 



(107) 



For a laser wavelength of 3 13 nm, and assuming e = 0.5 
and 77de, = 0.5, we have ^PJP^ > 2.3 X lO"' (PoTop)"'". 
For 1 W of usable laser power and Top = 1 (xs, this 
corresponds to a fractional power fluctuation of 
> 2.3 X 10"". 

This estimate applies only to the laser power fluctua- 
tions at the beamsplitter and assumes no additional noise 
is introduced between the beamsplitter and the photo 
detector or the beam splitter and the ions. Typically, the 
usable part of the laser beam must be directed further 
through optics, the air, and a window to the vacuum 
envelope enclosing the ion trap. Fluctuating etalon ef- 
fects in the optics and air currents may therefore seri- 
ously increase the power fluctuations beyond Eq. (107). 

The effects of (Gaussian) noise on laser intensity have 
been treated by Schneider and Milburn [274]. These 
effects show up in a well-characterized way for transi- 
tions involving Fock states. 

Fluctuations in timing errors may also be important. 
Clearly, fractional fluctuations in the duration of laser 
pulses will correspond directly to the same fractional 
fluctuations in the desired value of 6)„_„ . If we require 
fractional fluctuations of 10 " on 6)„„, then we require 
timing precision of 1 ps on a 1 (xs pulse. Similar consid- 
erations apply to the stability of pulse envelope shapes. 

For both amplitude and timing fluctuations, it may be 
possible to sample a portion of the beam used for logic 
and apply it to a "check bit" ion. The response of this ion 
could then be used to monitor and control the amplitude 
and timing of the pulses. 

4.4.2 Characterization of Amplitude and Timing 
Fluctuations 

In practice, it may be easiest to characterize the am- 
plitude and timing fluctuations with a power detector 
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placed as close as possible to the position of the ions. 
However, it would be useful to characterize these effects 
using the ions themselves. One possibility is to observe 
the sinusoidal Rabi oscillations predicted by Eq. (21). 
For example, for a single ion initially prepared in the | -l) 
state, we can record the probability Pi(t) of detecting 
the ion in the I -l) state after the laser is applied to the ion 
for time t. From Eq. (23), if we first prepare the ion in 
the \-l,n) state and if the power is constant, 
Pi(t) = I Cj,,„ h = cos^ I fl„;j I = 72[1 + cos2/2„.,„/]. In 
Fig. 2, we show this type of curve for n = and n' =\. 
This plot includes a decay due to decoherence. Here, we 
consider decoherence caused by specific types of noise. 
We will cast the results in terms of fluctuations in /2„,„ , 
however, they may easily be converted to corresponding 
time fluctuations since we are really interested in fluctu- 
ations in the net rotation angle given by Eq. (104). For 
simplicity of notation, in this section, we write /2„,„ = il. 
First assume that il fluctuates slowly so that it is 
constant over the time t for an individual measurement 
of Pi{T) but fluctuates over the time taken to make an 
average measurement of {Piir)). We characterize power 
fluctuations by a spectrum D(fl — (Jo) of £2 values 
around the desired value (Jo- We have 



'I 



{Pi(t)) = 72 [1 + cos2/2t]D(/2 - /2o)d(/2 - Oo). 



(108) 



We illustrate with two possibilities for D , finding 



(n(T)) = 2 



D(n- fk) = 



l+cos2/2Te-'(^"'"^>^ 
1 



for 






2ttA/2„ 



(109) 



{Piir)) : 



1 + cos2/2t 



1 



D{n- fk) = 



1 + 2(Aan,sT)^ 

V2lfl-A)l 



for 



AAr. 



where (A/2rms)^ = {{H ~ ^Y)- In both cases the contrast 
of the signal drops about a factor of two for t— 1/ 
A/2rms- Although the second choice of D(fl — (Jo) 
seems less physical, it indicates the exact spectrum of 
the noise may not be important for a qualitative test. If 
we need to detect fluctuations of order AO.^Jfl, then 
we must measure Pi after approximately /2t/2tt — 
/2/(2TTA/2rms) Rabi cycles. For example, if we require 
sensitivity ^.Q^mJH — 10"'', we need to make measure- 
ments of Pi{T) after about 1500 Rabi cycles. 

If, instead, the intensity fluctuates slowly compared 
to Top, but fast compared to the time t, {Pi(t)) is a less 
sensitive test of the fluctuations since their effects tend 
to average out. In order to show the effects of high-fre- 
quency fluctuations of £2 for a simple case, we consider 
a sinusoidal time variation: 



(2(1) = iTo + A/2sin(a)an,p; + q}). 



(110) 



a)amp is the modulation frequency of the Rabi amplitude 
(2(1). This analysis should, at least, yield the dependence 
on the amplitude and frequency of the fluctuations. The 
measured quantity Pi(t) is given by 



Pl(r) = hl+cose(r)], 



(111) 



where the rotation angle 9(t) is given by Eq. (104) 
which, here, takes the form 



eiT) _ r 

2 Jo 



nit)dt: 



flor+ AO 

Jo 



sin(«amp^ + (p)dt 



(112) 



An 



iloT+ [cos^(l — cosWampT) + simpsmw^^^T]. 



This yields 



n(T) = 2 



1 + cosi 2iloT+ [cos^(l — coscOampT) + sin^siiicOampT] 

\ ^amp 



1 + cos(2/2ot)cos( [cos^(l — coscOampT) + sin^sincOampT] 

\ ^amp 



sin(2iloT)sin( [cos^(l — coscOampT) + sin^sincOampT] 

\ ^amp 



(113) 
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Since Pi(t) is obtained by repeating the measurements, 
and since the phase ^ will be random from one measure- 
ment to the next, we average over ^: 



, . 1 cos(2/Jo 
(n(T))^ = 2 + 4^ 



^/>( 



2A/2 



[cos^(l — cosWampT) + sin^sinWampT] Id^ 



sin(2A,T) f" . (2An ^ .. X ^ ■ ■ n u 

sm[-;-| — [cos^(l — cosWampT) + sm^smWampT] )d^. 



(114) 



4t7 



\<War 



If we assume that I AH/w^^p I « 1 , which is valid for a 
small, high-frequency, modulation, then the integrands 
can be expanded to second order in (Afl/co^^p): 



, . 1 , cos(2/Jo 
(n(T))^-2 + 4^ 



sin(2/2oT) 

4t7 



•'0 

f'^f2An\ 



1 (2AnY , 

1 - ;t I I [C0S^(1 - COSWanipT) + sm^smWampT] 

^ \ ^amp / 



d^ 



[cos^(l — cosWampT) + sin^sinWampT] 



d^ 



(115) 



1 cos(2/2ot) 

'2'^ 2 



/An Y 

1 - 2i I (1 - COSWampT) 

V^amp / 



Aside from the fast oscillations due to the cosWampT 
term, which might be difficult to observe, the effect of 
the high-frequency fluctuations of 12 is to reduce the 
visibility of the signal by 2(Ai7/<yamp)^- This result can 
be averaged over a distribution Ai7(<yamp)- 



4.4.3 Applied Field Frequency and Phase 
Fluctuations 

A simple source of fluctuations is due to frequency or 
phase fluctuations in the radiation that is used to carry 
out the basic operations. Schneider and Milbum [274] 
have considered a specific model of phase fluctuations 
for ion experiments and calculate the corresponding 
decoherence for operations of the type used in quantum 
logic. If qubit energy levels are separated by optical 
energies, the lasers that drive qubit transactions must 
have the required frequency and phase stability. Given 
the performance of current stable lasers [235,245,246, 
247], this may be a problem for long computations. As 
discussed in Sec. 2.3.3, optical stimulated-Raman trans- 
actions provide the required strong field gradients 
whereas the overall frequency sensitivity in the transi- 
tions is dominated by the difference frequency of the 
lasers, not the frequency of each laser. Since the two 
Raman beams can be derived from one laser beam with 
the use of frequency modulators, the frequency fluctua- 
tions are dominated by those of the oscillator which 
drives the modulator. The phase stability of these 
sources can be high and does not appear to pose prob- 
lems. For example, the data shown in Fig. 7 can be 



regarded as one test of the phase sensitivity of the 
oscillator used to drive the indicated internal state tran- 
sition. Stimulated-Raman transitions have the disadvan- 
tage that they give rise to ac Stark frequency shifts as 
indicated by the Igi^rMR terms on the right side of 
Eqs. (40). If these shifts are equal for both qubit levels, 
the net shift is zero. If they are different, the effective 
qubit frequency is shifted during the operation. This 
must be measured and accounted for. The shift can be 
tuned to zero by appropriately adjusting Ar and the 
relative intensities and polarizations of the Raman 
beams. A problem still exists if the Stark shifts fluctu- 
ate; we examine the consequences of those fluctuations 
here. 

From Eqs. (40), the ac Stark frequency shift of the 
qubit transition due to stimulated-Raman transition is 



COS = C0S2 - (Osi = - (\ gll 



lgih)//iR, (116) 



so that the shifted qubit frequency is given by 
coo= (Oo+ cOs- We first consider the effects of frequency 
fluctuations of the Raman beams. If the two beams are 
derived from the same beam with the use of a modula- 
tor, the frequency shift of the two laser beams will be 
the same. Therefore, to a good approximation, we find 
8<yi — — (8Ar/Ar)cOs where 8z\r represents the fre- 
quency shift of both laser beams. Typically, Ar > 10 
GHz, and frequency fluctuations of lasers can be con- 
trolled to less than 1 kHz, so this source of error should 
be small. Moreover, it is essentially absent if we tune co^ 
to be zero. 
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More serious are fluctuations in laser intensity. We 
will characterize these fluctuations by the correspond- 
ing fluctuations in g, (oc square root of laser intensity). 
From Eq. (116), we have 



8(ys = 2ws2 — - - 2wsi — ^ ■ 

g2 g\ 



(117) 



As a worst case, we will assume fluctuations in g\ and 
g2 are random and uncorrected and, for simplicity, we 
assume {{^gilgif) == {{^gjgif) = ^^ where () denotes 
an average over the spectrum of the fluctuations. There- 
fore 



((8«s)') = 4^\«|, + wh). 



(118) 



As discussed in Sec. 4.3.2, we assume the fluctuations 
are slow compared to T„p, where T„p is the time of the 
operation. In a particular operation, a frequency offset 
Sws gives rise to a phase fluctuation 8(^s = SwsTop. It is 
useful to characterize the effects of these phase fluctua- 
tions by comparing their size relative to fluctuations in 
rotation angles caused by the corresponding fluctuations 
in the Rabi rates /2„',„ . With the same assumptions re- 
garding 8gi and 8g2, we have ((8/2„,„//2„.,„)^) = 2^^. 
Therefore, fluctuations in the overall rotation angle 
0„.,„ = 2/2„..„Top [Eq. (104)] are given by ((80„.,„/ 
&n;nf) = 2^^. From these expressions we find 



((8(^s)') Wl + oi2) 
((80„..„)^) 2/2\..„ 



(119) 



For quantum logic, a worst case appears to be for side- 
band excitation n' = n±\. If we assume the Lamb- 
Dicke limit, we find 



((80„.,„)') ■ 



ist+gt) 



(120) 



From this expression, it is desirable to have gi — g2 to 
minimize the effects of 8(^s relative to 80„.„. Even in 
this case, rms phase fluctuations caused by Stark shifts 
are worse than those caused by Rabi frequency fluctua- 
tions by a factor of — I/tj. Most likely, however, fluctu- 
ations in intensity will be dominated by fluctuations 
from the primary laser from which both beams are 
generated. Therefore we expect fluctuations to be corre- 
lated 8gi/gi = ^gi/gi- When, in addition, wsi — cosi, the 
phase fluctuations caused by Stark shifts will be less 
than those caused by Rabi frequency fluctuations. If 
these conditions hold true, it appears that stimulated 
Raman transitions between two qubit levels separated 
by fairly low (for example, hyperfine) frequencies are 
superior to single-photon transitions between qubit lev- 



els separated by optical frequencies. This situation 
might change as laser oscillators become more stable. 

In the above discussion, we have neglected the alter- 
ation of the ion motional frequencies caused by the 
superimposed dipole force potentials of the focussed 
laser beams. However, if we assume the Raman beams 
have waists Wo which are approximately equal and that 
the magnitudes of the projections of their k vectors 
along the z axis are the same, it is straightforward to 
show that the relative shift of (single trapped ion) secu- 
lar frequencies for states I -l) and I T) is approximately 
equal to hcos(zo/wof- Therefore, since zo « Wo, the fre- 
quency shifting effects will be dominated by the shifts 
of the internal states. 

Phase fluctuations of the laser beams themselves will 
also directly affect the fidelity of the operations. In the 
stimulated-Raman case, where both beams are derived 
from the same beam, phase fluctuations between Ra- 
man beams will be very nearly cancelled. However, an- 
other source of laser phase fluctuations will come from 
path length fluctuations between the (laser) source and 
the trapped ions. Path length fluctuations are expected 
to be dominated by mechanical vibrations; these vibra- 
tions are typically restricted to low frequencies (< 1 
kHz). They could result from a number of causes such 
as fluctuating mirror mounts or trap mounting hard- 
ware. For single photon laser transitions, the overall 
path length between the laser and ions is important; for 
stimulated-Raman transitions, the primary problem will 
be caused by path length differences between the two 
Raman beams after the frequency modulator. Therefore, 
for stimulated-Raman transitions using overlapping, co- 
propagating beams (which can be used to drive carrier 
n' = n transitions) the paths are the same and there 
should not be a problem. However, for stimulated Ra- 
man sideband transitions, we require ki and k2 along 
different paths, so that the path length problem is 
analogous to the problem for single photon transitions. 
For brevity, we will treat the problem of single-photon 
transitions; other cases (including the stimulated-Ra- 
man case) are easily generalized from this. 

We assume the laser electric field at the exit of the 
laser oscillates as cos(coi^t). At the position of the ion, 
the field is given by the same expression with the time 
replaced by the retarded time t — djc where di, is the 
distance between the laser and ion. Therefore, the overall 
phase difference between the field at the laser and ion is 
co^diJc . We will assume the fluctuations in d^ are slow 
enough that over the time of a single operation the phase 
can be considered constant. If d^ fluctuates by an 
amount Adj^, the phase of the field fluctuates by Adi^oiJ 
c = 2'nA.diJX where A is the laser wavelength. In the 
NIST single 'Be"^ ion experiments, all of the operations 
have typically taken less than 1 ms; therefore vibrations 
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have not been a problem. However, in longer computa- 
tions, the requirements on Ad^ are stringent. Some form 
of active stabilization, such as the method described by 
Bergquist et al. [235], will probably have to be used. 
That technique can be viewed as a Doppler shift cancel- 
lation scheme and was based on a Doppler shift cancel- 
lation scheme used in spacecraft tracking [275]. To il- 
lustrate what is, in principle, possible, we note that 
optical cavities can be made to track the frequency of 
lasers to precision much smaller than 1 Hz 
[235,245,246,247]. For a laser frequency of 5 X 10'" 
Hz (A = 600 nm), and a cavity composed of two mirrors 
separated by a distance of 50 cm, this corresponds to 
holding relative positions of the mirror to within 10~" 
cm. Therefore, although they add additional complica- 
tions to the experiments, such schemes for length stabi- 
lization can be used to effectively null the effects of path 
length variations. 

4.4.4 Individual Ion Addressing and Applied Field 
Position Sensitivity 

The scheme of Cirac and Zoller for trapped ion quan- 
tum logic requires that ions (along the axis of a linear 
trap) be addressable individually with laser beams for 
logic operations. This may be difficult, because the high 
vibrational frequencies desired for efficient laser cool- 
ing and suppression of decoherence also results in 
closely spaced ions. As discussed in Sec. 2.1, the mini- 
mum separation of adjacent ions in a linear trap is be- 
tween the center ions and is approximately 
i'min — 2sL~°'^'', with s = {q^lA'n€omwff'^ where L is the 
number of ions. For 'Be"'^ ions with an axial COM fre- 
quency of wJ2'n = 1 MHz, this separation is about 4 [xm 
for 10 ions. 

The most straightforward method for individual opti- 
cal addressing is to tightly focus laser beams on the 
selected ion [1]. The transverse intensity distribution of 
a Gaussian optical beam of power P is 



/(r) = -2exp 2 

TTWo \ Wo 



(121) 



where Wo — \l{ii-NA) is the beam waist, A the radiation 
wavelength, and NA = tanfl is the numerical aperture of 
the beam with cone half-angle B (the formula for Wo in 
the paraxial ray approximation is valid only for NA < 1) 
[276]. For large numerical apertures {NA ~ 0.5), beams 
can apparently be focus sed to Wo — A, but this is diffi- 
cult to realize in the laboratory. If Wo = 5 [xm in a Gaus- 
sian beam, at a distance 10 [xm from the center of the 
beam, this would imply a relative intensity of about 
3 X 10"'' (relevant for stimulated-Raman transitions) or 
an electric field amplitude (relevant to single photon 



transitions) of 1 .8 % relative to the center of the beam — 
clearly a problem. If Wo = 2 [xm could be obtained, the 
intensity (electric field) would be down by a factor of 
1.3 X 10"'" (1.1 X 10"'). However, these estimates are 
too optimistic because they assume the laser beam is 
normal to the axis of the linear trap; to address the axial 
modes by the methods described in Sees. 2 and 3, we 
need a component of the laser beam k vector along the 
axis. Therefore, if the angle of the laser beam relative to 
the trap axis is ^k, we must replace r in Eq. (121) by 
i'minSint^j:. In addition, imperfections in the surfaces of 
the intervening vacuum port window, multiple reflec- 
tions from these windows, and diffraction typically dis- 
tribute laser intensity outside of the theoretical waist of 
the beam. The degree to which this occurs depends on 
the details of window surfaces, etc. and must be resolved 
experimentally. 

Alternatively, we could accomplish all of the opera- 
tions described in Sees. 2 and 3 by using the transverse 
gradient of the field associated with a focused beam 
[125]. Therefore to accomplish sideband transitions, we 
would displace the laser beam laterally, with respect to 
its k vector from the position of the ion. This would give 
rise to a coupling to the ion's motion which is in the 
direction of the transverse intensity gradient of the laser 
beam. In this case, we can make the direction of the 
beam normal to the trap axis. This method has the 
disadvantage that the field gradient would be reduced 
relative to the case treated above thereby leading to 
smaller Lamb-Dicke parameters and correspondingly 
reduced values of /2„.„ for n' i= n. 

The transverse intensity gradients of focused laser 
beams can also cause significant intensity fluctuations 
at the selected ion if the relative position of the beam 
with respect to the ion is not stable on the time scale of 
the computation. An alternative to using tightly focused 
Gaussian laser beams is to first feed the (expanded) 
laser beam through a sharply defined aperture, and use 
a lens to image the aperture at the position of the ions. 
With this technique, the beam intensity can be dis- 
tributed more smoothly around the selected ion and have 
very steep intensity edges (on the order of the original 
aperture sharpness) away from the ion, thus suppressing 
beam vibration problems and confining the radiation to 
a single ion. This technique has been used to make 
relatively "hard" walls for an optical dipole trap [277]. 
For this technique to work well, the imaging lens must 
collect a large fraction of the light transmitted through 
the aperture or else diffraction effects will result in light 
intensity outside the image of the aperture. To address 
individual ions, we require very small aperture images, 
which gives rise to a design tradeoff. If a one-to-one 
relay lens is used to image a small object aperture, 
effects of diffraction are enhanced. If a demagnifying 
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lens is used to reduce a large object aperture, then the 
aperture must be placed a large distance from the lens, 
requiring a relatively large lens. For two ions, imaging 
a sharp edge such as a razor blade at the space between 
the ions may be sufficient. We might also consider hav- 
ing every other ion in a string be a "garbage" ion which 
is not used in the computation, thereby increasing the 
spacing between qubit ions by a factor of two (or more, 
if more garbage ions are used between each qubit ion). 
This has the disadvantage that total number of ions (and 
spectator modes) increases, aggravating the problems 
associated with large quantum registers. If sufficiently 
good addressing on one ion in a string can be accom- 
plished, it may be simpler to adiabatically shift the posi- 
tion of the ions, rather than the laser beams, in order to 
address different ions. This could be accomplished by 
applying different static potentials t/o(0 and t/o(0 to the 
end segments of the rods in Fig. 1 . However, changes in 
t/o(0 and Uo{t) would have to be coordinated to keep the 
COM axial frequency constant or else additional phase 
shifts would be introduced. Stimulated-Raman transi- 
tions have the advantage that Afc can be made parallel to 
the axis of the trap even though each beam is at an angle 
with respect to the trap axis to facilitate selection of a 
particular ion. This is important, since coupling to trans- 
verse modes is eliminated. 

Another method of optically addressing individual 
ions is to cause a destructive optical interference at the 
position of a specific ion, with a net coupling at the 
other ion(s). For instance, if ion j is positioned at the 
node of a resonant standing wave laser field, the cou- 
pling between states I 4-); | n^) and I \)j \ rik) is propor- 
tional to {nk I sm['r]i{ai, + al)\ \ rik). In this case, the cou- 
pling of the standing wave to ion j vanishes when the 
laser frequency is tuned to an even order sideband such 
as the carrier («<.' = m<.) (also see Sec. 4.4.6.2). If, in- 
stead, the ion is positioned at an antinode, the coupling 
is proportional to (n^' I cos[i74(a<: + a^t)] I M/t); thus, the 
coupling vanishes when the laser frequency is tuned to 
an odd order sideband, such as the first blue or red 
sideband (n^' = n^ ± 1). By appropriately choosing the 
angles of focused laser beams relative to the trap axis or 
the spacing between ions, it should be possible to posi- 
tion an antinode (node) at ion 7 while positioning nodes 
(antinodes) at the ions adjacent to ion j (for equally 
spaced ions). In the case of two-photon stimulated Ra- 
man transitions, we desire to place ion j at common 
nodes or antinodes of two standing waves. Although this 
interference technique should allow individual access to 
each of two trapped ions, it appears technically difficult 
to extend this technique to more than three ions. 

Finally, we consider the application of external field 
gradients which shift the internal energy levels of ions 
depending on their position. For a magnetic field gradi- 



ent to give this selectivity, we require the Zeeman split- 
ting between adjacent ions to be much larger than the 
Rabi frequency, or Afjb(d I B \ ldz)sl/i » fl, where Afjuis 
the difference in {ijl-B)/\B \ between the two levels of 
interest, and s is the ion separation along the z direction. 
For Ajx = jXb, s = 10 (xm, and a Rabi frequency of O/ 
2tt = 1 MHz, this requires d\B\/dz»0.l T/cm. Field 
gradients of this magnitude can be achieved; however, 
they would introduce large, and not-easily-controlled 
phase shifts for the other ions in a quantum register. 

The laser beam itself can provide ion selectivity by 
employing the transverse gradient in the optical field 
intensity. For instance, if we desire to perform a 0-pulse 
on ion J without affecting neighboring ion k, the inten- 
sity profile of the laser beam can be set so that the ratio 
of field strengths (intensities for the case of two-photon 
Stimulated-Raman transitions) at ion j vs. ion k is 6/ 
iTTm , where m is an integer. Now if the pulse duration 
is set so that ion j is rotated by 6, ion k receives a rotation 
of iTTin and hence returns back to its initial state (with 
an extra phase factor of (— 1)"). 

For the case of two-photon stimulated-Raman transi- 
tions, the laser beam can provide ion frequency selectiv- 
ity by employing the Stark shift and the transverse gra- 
dient of the optical field. Here, for example, we could 
assume that two counterpropagating Raman beams of 
equal intensities and spatial profiles are offset so that 
beam 1 is centered on ion j , and beam 2 is centered on 
adjacent ion k as depicted in Fig. 8. We will assume the 
coupling scheme of Fig. 3. Let e be the fraction of peak 
intensity seen by the offset ions (that is, the intensity of 
beam 2 at ion j and beam 1 at ion k). Assume that when 
both beams are positioned on either ion, gi =g2 = g- 
When the beams are offset, the two-photon resonant 
Rabi frequency at each ion is il= e"^{g^/AR), where 
sVAr is the Rabi frequency expected if both beams were 
centered on a given ion. The Stark shifts of the two ions 
are in opposite directions: S, = + So, 81,= — So, where 
So = n{\ - €)/€"\ If we make So » H (e « 1), then 
by appropriately tuning the difference frequency of the 
laser beams, we can selectively drive transitions on ei- 
ther ion 7 or k. Alternatively, if, for example, we desire 
to perform a 0-pulse to ion j without affecting ion k in 
an "unrepairable" way, e can be adjusted to a particular 
value which results in ion k making the transformation 
.A, = « U), + /3 I T), ^ ,A; = e^ia I i\ + e'"f/3 | T),). For 
example, if we use Rabi pulses which have a square 
pulse envelope, Eq. (21) applies. For ion j, the laser 
must be tuned to the shifted frequency cojo = cojo + So; 
Eq. (21) applies with A = and | il„;„ 1 1 = 6/2 where t is 
the pulse time. For ion k, w^ = co^ — So; Eq. (21) ap- 
plies with A = 28o. To achieve the desired form of i//^, 
we want 6/2(1 + S"//2")"" = nnr, or e' - [1 + (2mtt/ 
6y]e+ 1=0, where m is an integer. For m=l and 
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Fig. 8. Scheme for differential ac Stark shifting of neighboring ions. Equal intensity counter propagating beams 
1 and 2 are centered on ionsy and k, respectively. A fraction £ of the peak intensity la of each beam is applied 
to the other ion. This results in a differential ac Stark shift of ionsy and k, allowing the possibility of individually 
accessing the ions by tuning the frequency of the laser beams. 



0=17 (a TT-pulse on ion j), we find €=0.208, and 
^= 1.7417. The phase shift ^ on ion k and the corre- 
sponding phase shift on ion j must be kept track of in 
subsequent operations on these ions. Generalizing this to 
more than two ions becomes difficult if the laser beams 
also overlap other qubit ions. This scheme places an 
additional premium on laser power stability, since the 
light shifts are bigger than the Rabi frequencies by l/e"^ 
for €« 1. In addition, in both of the above schemes 
employing the laser beams to differentially address 
neighboring ions, one major drawback is that the posi- 
tions and profiles of the laser beams must be accurately 
controlled. 

Many of the above differential addressing schemes 
appear useful when dealing with only two ions instead 
of a string of many. This leads us to seriously consider 
systems where quantum logic operations are performed 
on accumulators consisting of only two ions, with the 
other ions located somewhere else (Sec. 5.1). A pro- 
posal has also been made for transferring quantum in- 
formation from one register to another by optical means 
[151,278]. 



4.4.5 Effects of Ion Motion (Debye- Waller 
Factors) 

The Rabi frequency fln\n describes the transition rate 
between states 1 4-) I n ) and | T) I n'). To realize the condi- 
tional dynamics desired for quantum logic, we want /3„,„ 
to depend on the motional quantum numbers n and n' of 
a particular vibrational mode and be independent of the 
state of motion of other modes. In addition, for simple 
rotations on internal states, we want /3„',„ to be indepen- 
dent of the motional state for all modes. It is not possi- 
ble, in practice, to rigorously satisfy both of these 
requirements. For instance, the conventional controlled- 
not gate employs two carrier pulses (steps (la) and (Ic) 
in Sec. 3.3) which are intended to not depend on the 
state of motion; this requires the Lamb-Dicke parameter 
7] to be small [see Eq. (56)]. In the Raman configura- 
tion, 7] is proportional to the difference in two wavevec- 
tors and can be made negligible by using co-propagating 
beams (AA: — 0). On the other hand, with single-photon 
optical transitions, the Rabi frequencies depend on the 
motion of all modes which have a component of motion 
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along the direction of A:. As discussed in Sec. 2.3, we 
can take advantage of the motional dependence of the 
carrier to construct a logic gate, but in this case also, the 
Rabi frequency will depend on the motion in the other 
modes along the direction of k ot Ak. Similarly, for 
sideband operations, such as step (lb) in Sec. 3.3, it will, 
in general, be impossible to have i7„,„ , depend on only 
one mode of motion. In this section, we examine the 
influence of extraneous modes on the Rabi frequencies 

Transitions will generally be affected by motion in the 
3L-\ spectator modes of vibration. The motion in these 
other modes reduces the Rabi frequency in much the 
same way as lattice vibrations affect a single emitter or 
scatterer embedded in a crystal, as described by the 
Debye- Waller effect [279,280]. From Eq. (36), we have 
the Rabi frequency for thejth ion (assuming all modes 
are in specific Fock states) 






{{n,^k}\Y\^"'''"''"''''\{Hp*k}) 



pi=k 



(122) 



pi'k 



" e i-z 
2z"i?„(x) = -j 

n=0 V — Z 



/!=o I — Z 



2xV~z 
1 - 



(125) 



In these expressions, I z I < 1 and /o(e) is the zeroth 
modified Bessel function with argument e. 

The first expression in Eq. (124) shows the exponen- 
tial reduction of the Rabi frequency due to motion in 
modes p (Debye-Waller factor). For each mode p, we 
can write (rfp)\hp + 'h) = Qlik^^xlY) where k^a is the 
effective k vector of the radiation and xj, is the amplitude 
of the component of motion parallel to k^ff for the j\h 
ion. For L large, the product term on the right side of the 
first expression is approximately equal to exp(— 7 
2 ^eff(x/ms)^) where x4,s is the total rms amplitude of mo- 
tion parallel to fceff. This is the Debye-Waller reduction 
factor due to the thermal energy of the jth ion. 

The second line in Eq. (124) allows us to determine 
the fractional fluctuations in the Rabi frequency from 
experiment to experiment 



where flUt^m = fl\{nl\ exp(ii7i(aj. + al)) \n,,)\ is the Rabi 
frequency of the ^th mode (selected for logic), in the 
absence of the 3L-\ spectator modes labeled by index/? 
and vibrational number rip. We assume the vibrational 
states of spectator modes are independently thermally 
distributed with mean vibrational number hp : 



"" (i + «p)"' 



(123) 



If we average over this distribution, the mean and mean- 
squared Rabi frequencies are given by 



p'/=k'-np=l \ 



(124) 



pi=k'-np=l \ 



: {il^l^ n e-^<'^'^<V>'^> /o( 2{j^ifyhp{hp+l) 

pi=k 



where we have used the Laguerre polynomial sum iden- 
tities 






'^■^"i>t 



JAi'.nt 



^k' 



-pi=k 



Y\h(2{vif^np{hp+\) 



= A/E {'rfp)\inp+ 1). 



(126) 



pi=k 



In this last approximation, the Bessel function is ex- 
panded to lowest order [/o(e) = 1 + eV4 + ■■■] which is 
appropriate if the 3L-\ arguments of the Bessel function 
in Eq. (126) are all small compared to 1. This is ex- 
pected to be the case if all modes are cooled to the 
Lamb-Dicke limit [(')7^)^np« 1]. Eq. (126) describes 
the fractional rms fluctuations in the Rabi frequency 
due to thermal motion in the spectator modes of vibra- 
tion. 

Typically, the thermal motion is determined by initial 
conditions and the reservoir is "turned off once the 
experiment begins. This is expected to be the case for 
imperfect initial laser-cooling, resulting in a probability 
distribution of stable Rabi frequencies with the above 
mean and rms values. In this case, the Rabi frequencies 
for each ion maintain a constant value during a single 
run of the experiment, but deviate from the mean as 
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indicated in Eq. (126). Stating the results in another way, 
if the total number of modes is large, then the distribu- 
tion of Rabi frequencies is nearly Gaussian, and the 
probability that a given run of the experiment results in 
a Rabi frequency whi ch is fractionally within e of the 
mean Rabi frequency /2„j',„j is 



Pr\ 



■erf\ 







fl' . 




D,' . 



<€ 






(127) 



For example, if each of 100 spectator modes is laser- 
cooled to hp = 0.1 and has Lamb-Dicke parameter 
rfp = 0.01, then the probability that the Rabi frequency 
is within 10""* of the mean value is approximately 0.23. 
To see how Eq. (126) scales with the number of ions, 
assume that fceff is aligned with the axis of a linear trap, 
where the axial motion is described by L modes. If the 
frequencies and amplitudes of all modes contributing to 
the axial motion of ion j are assumed to be about the 
same, we can write rjj, — rjJL "^ and fVp — h where 171 
and h are the Lamb-Dicke parameter and mean occupa- 
tion for the axial motion of a single (thermalized) 
trapped ion. In this case, Eq. (126) becomes 



A 0*™> 



n(n + 1) 



(128) 



This expression indicates that a large number of ions is 
beneficial because it tends to average out the effects of 
motion in the (L — 1) extraneous modes. Equation (128) 
is an overestimate of the fluctuations since the L — 1 
extraneous modes will have higher frequency than the 
COM mode, leading to smaller amplitudes of motion 
than assumed in this crude estimate. 

The Rabi frequency fluctuations discussed here trans- 
late to fractional phase offsets [{^J&k of Eq. (94)] for 
each gate operation as indicated in Eq s. (124 ) and (126). 
In practice, the mean Rabi frequency fi„^\nt can be mea- 
sured by averaging over many experiments. If we as- 
sume the COM mode is used for logic, then 17^ = tji/L "^. 
Therefore, for L » 1, we cannot choose 171 to be too 
small, or else the entangling operations for quantum 
logic become too slow (Eq. (24)). To suppress the ef- 
fects of Debye- Waller factor fluctuations, it is therefore 
desirable to cool all modes (whose motion is parallel to 
A:eff) to the zero-point state. 



4.4.6 Coupling to Spectator Levels 

We have assumed that each time an external field is 
applied to form part of a logic operation, only two quan- 
tum states take part in the interaction. This assumption 
is valid when any other state is far from resonance. To 
explain this in a simple example, we refer to Fig. 9. 
Suppose we want to drive a transition between states I 4-) 
to I T) with radiation resonant with this transition. If 
level I 4-) is coupled to level I s) with this radiation, the 
wavefunction will have a small admixture of state \s) 
after the operation; the amplitude of state \s) will be- 
come larger as A becomes smaller. After a sequence of 
many operations, the amplitudes of these "spectator" 
states can build up and cause errors in the computation. 



s> 



t> 



i; 



r, 



i) 



Fig. 9. Simplified energy level diagram for characterization of cou- 
pling to, and spontaneous emission from, off resonant "spectator" 
levels. We assume coherent radiation is tuned near the transition fre- 
quency for the |4-)|n)— > |T)|n') transition (for simplicity we have 
not shown the motional substructure for the I i) and I T) states). The 
state U) is assumed to be the nearest spectator state from which 
off-resonant coupling can occur. 



4.4.6.1 Polarization Discrimination of Internal 
States 

In an array of trapped ions, two internal states of each 
ion comprise a qubit of information. In addition, a third 
auxiliary internal level | aux) may be required (tran- 
siently) for the operation of a CN gate as described in 
Sec. 3.3. This state might be a particular Zeeman sub- 
level of a hyperfine multiplet [17] or an optical 
metastable state [1]. By employing suitable polarizations 
of the driving field, particular internal state transitions 
can be selected with high discrimination. Furthermore, 
the Zeeman sublevels can be spectrally resolved (assum- 
ing the Zeeman splitting is much larger than the Rabi 
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frequency), by applying a magnetic field. This combina- 
tion will help to isolate the two internal levels of interest 
from other internal levels. Here, we give an example of 
how polarization selection can be used to suppress cou- 
pling to internal "spectator" states. 

We consider the case of Ref. [17], which uses two 2s 
^Si/2 ground state \f, Mf) hyperfine levels of 'Be"^ ions 
as qubit levels. These levels are designated I -1-) = 1 2,2) 
and I T) = 1 1,1) (Fig. 5). We also make the identification 
I aux) = 1 2,0). They are driven using stimulated-Raman 
transitions by detuning the lasers from the ^Fui state. By 
choosing the Raman beam coupling between the ^Pig 
and I -l) states to be cr polarized and the coupling be- 
tween the ^Pi/2 and I T) states to be tt polarized, the | -l) 
and I T) states form a closed family; that is, neither the 
I -1-) or I T) state can be driven to other Zeeman sublev- 
els. In this case, the nearest off-resonant transitions to 
consider are detuned by at least the hyperfine frequency 
(— 1.25 GHz in 'Be"^). Similarly, by choosing the beam 
coupling between the ^Pic and I aux) states to be cr* 
polarized (and the coupling between the ^Pic and | T) 
states to be tt polarized), the | aux) level is coupled only 
to the I T) level. 



H, = i^Wo\^)(t\ +^iwo + A)\s){s\, 

(129) 
H, = 2^coswtin[\i){1\ + |T)(i|] 

+ /2'[U)<sl + \s){i\]), 

where the expression for Hi is equivalent to that in Eq. 
(14). In the Schrodinger picture, we write ^= QI -1-) + 
CTexp(— iwot) I T) + Csexp(— ; [wo + A]t)\ s). This leads 
to equations for the amplitudes 

a=-//2e*CT-(/2'e'<«"^>'Cs, 

Ct = - ine-'"Ci, (130) 

C,= -//2'e-'<«-^"Q, 

where S= co — coo and we have neglected rapidly vary- 
ing terms proportional to exp(± i(co+ coo)t) (the usual 
rotating wave approximation). If we make the substitu- 
tion Cs = Cfi'^', the equation for C^ can be written 



Cs + iAC=-in'e-""Ci. 



(131) 



4.4.6.2 Spectral Discrimination of States 

Spectator states can include both motional and inter- 
nal states. For example, a single CN gate uses pulses 
which drive on the carrier as well as first motional 
sidebands of the COM mode. Discrimination between 
these two types of transitions cannot be made with 
polarization selection since the relative strengths of the 
matrix elements are independent of polarization. Since, 
in the Lamb-Dicke limit, the carrier operation will have 
a higher resonant Rabi frequency than the sideband op- 
erations, the largest source of contamination may be due 
to an off-resonant excitation of the carrier during the 
sideband operations. The basic problem can be illus- 
trated by considering the simple example illustrated in 
Fig. 9. Here, we assume decay from the excited states is 
negligible, ^T^ = rV = 0. We assume we want to carry out 
an operation which coherently drives the I -1-) -"-^ IT) 
transition (resonance frequency coo) and avoid the pres- 
ence of any amplitude in the "spectator" state | s) after 
the operation. We assume radiation of frequency co—Wo 
is applied which couples levels I -l) and | T) with 
strength O. We assume this radiation also couples levels 
I i) and I s) with (resonant) strength {}', but, for sim- 
plicity, does not couple levels | T) and I s). In the context 
of quantum logic, fl might correspond to the Rabi fre- 
quency for a sideband transition and O' the Rabi fre- 
quency for the carrier transition. If we take the zero of 
energy in Fig. 9 to be the energy of the I -l) level, we can 
write the Hamiltonian as H = Ho + Hi, where 



For I {}' I small enough relative to 1 4 | , we can 
"adiabatically eliminate" level s, by assuming |dCs7 
df I « I ACsl . In this case, Eqs. (130) become 



Ci 



ine'"Cr + i^^ Cu 



Ct = - iOe-'^'Ci, C =- - ^ e'-'^'Ci. (132) 



The Stark shift term S^ = (O'^/A is a downward energy 
shift of the I -l) state. It can be suppressed in the previ- 
ous equation by including it into a shift of Wq. That is, by 
making the substitution Cj, = Ciexp((5sf ) and choosing 
CO — coo= 5s, we find 



Ci=-inCt, Ct=-inCi, C=-^Ci. (133) 



In this case, Cl and Ct are given by Eqs. (17), leading 
to a the desired evolution between the I i) and I T) 
states. Moreover, if the couplings fl (and O') are turned 
on and off slowly compared to 1/A (the condition of 
adiabaticity), then the amplitude in state | s) adiabati- 
cally grows during the operation and reduces to zero 
upon completion of the operation, since it is propor- 
tional to il'. 

The general problem, however, is that we desire to 
make all transitions as fast as possible. This means that 
we cannot satisfy the adiabatic condition; therefore Cs 
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will not be simply proportional to 12' and will have some 
nonzero value after 12' is reduced to zero. For example, 
if 12' is turned on adiabatically, but turned off diabati- 
cally, Cs will have a final amplitude with magnitude 
I fl'Ci/A I . Estimates of these effects in specific con- 
texts are given by Poyatos et al. [261] and James [61]. 
Since it should be possible to separate internal states 
by relatively large frequencies and further discriminate 
from internal spectator states with polarization, the most 
important task will be to suppress transitions to mo- 
tional spectator levels. The fundamental problem is that 
in order to discriminate between carrier and sideband 
transitions, the ion must be able to tell that it is moving. 
If the time of the operation Top is small compared to the 
inverse of the frequency of the motional mode, this is 
impossible. Therefore, the fundamental limit on Top is 
given approximately by 1/cocom- To maximize computa- 
tional speed, we desire to approach this condition as 
closely as possible. This can be aided by using Rabi 
envelopes which can suppress certain spectral compo- 
nents [281,282]. Since the general problem will involve 
more than one spectator level (more motional states), a 
more general approach might be to apply optimal con- 
trol theory [283-286]. 

4.4.6.3 Tailoring of Laser Fields 

As discussed in the last section, the speed of an ion- 
trap quantum computer will be fundamentally limited 
by excitation of motional sidebands other than the de- 
sired one. To see this in a specific case, the conditional 
dynamics for the first CN gate described in Sec. 3.3 
occurs in step (lb). That step uses a lir pulse on 1 st blue 
sideband of the I T) 1 1) <-> I aux) 1 0) transition. To make 
this step fast, we would like to increase the field inten- 
sity and decrease the interaction time accordingly. How- 
ever, as the pulse becomes shorter, the spectrum of the 
pulse becomes wider. Eventually, the spectral width of 
this pulse is larger than co^ so that we drive the relatively 
strong and unwanted I T) 1 1) -<->■ I aux) I 1 ) carrier transi- 
tion. In general, we also drive I T) I 1 ) -<->■ I aux) In') tran- 
sitions, where n'> \. Most of these unwanted transitions 
can be suppressed by making all terms in the expansion 
in of Eq. (28) equal to zero except the term proportional 
to dB/dz (or dE/dz for electric dipole transitions). In this 
case, we suppress coupling to all unwanted I aux) levels 
except the I aux) 1 2) state whose resonant frequency is 
detuned from the I T) 1 1 ) -<->■ | aux) 1 0) transition by 2(0^ . 

Suppression of even or odd order terms in Eq. (28) 
can be accomplished by appropriate positioning of 
standing waves relative to the ion [40,61,287,288]. To 
illustrate with an example, we first examine the transi- 
tion I si) 1 0) ■<->■ I T) 1 1). For simplicity, we assume single 
photon electric-dipole (rather than stimulated-Raman) 



transitions and consider only motion along the z direc- 
tion. We want to suppress all terms in the expansion of 
the field [Eq. (28)] except the dE/dz term. We need to 
synthesize, using M separate standing wave laser beams, 
a wave which has the form 



E = Eo 



2 C„sin(^„z) 



cos(&)0, 



(134) 



where k„ = k„- z, and where z = corresponds to the 
equilibrium position of the ion in question. The case of 
M =\ has already been suggested in various contexts to 
suppress coupling to the carrier (An = 0) transition (see 
Ref. [40] and Sec. 4.4.4 above). For M> 1, we want to 
satisfy 



C,ytr+C2fe"+... + Cm^£ = 0, 



:3,5,...2M -1.(135) 



In this case, the next sideband spectator level which has 
nonzero coupling is one where \n' — n\ = 2M + 1 . This 
transition is detuned by 2Mco^ from the transition of 
interest and the coupling is suppressed by approximately 
7]^'^ = {kiZof^ compared to i7i,o or /2o,i. 

Second, we consider carrier transitions |sl)|n)-<->- 
IT) In), where m = 0,1. With the stimulated-Raman 
technique on hyperfine transitions, transitions of the 
type |sl)|n)-<->- |T)|n'), n + n\ are highly suppressed 
by the use of copropagating beams. For single photon 
transitions, we need to synthesize, using separate laser 
beams, a wave which has the form 



£ = £o 



2 C„cos(^„z) 



cos.{wt). 



(136) 



For M> 1 , we want 



C,yt,"+C2fe"+... + CM^i& = 0, 



n = 2,4, ... 2(M- 1). 



(137) 



In this case, the next spectator level which has nonzero 
coupling is one where \n' — n\ = 2M. This transition is 
detuned by 2Mco^ from the carrier and the coupling is 
suppressed by approximately 77™ = (kiZof^ compared 
to i7i,i or iToo. For the carrier transitions, we do not have 
the problem with the remaining spectator level as noted 
for the first sideband case above. On the other hand, for 
logic operations which use the carrier transition (Sec. 
2.3 and Ref. [174]), we do not want rj to be too small, 
or else the gates take to long to implement. Therefore 
the suppression of the higher-order sidebands may not 
be as great as in the case of logic using first sidebands. 
For qubits coupled by single-photon transitions, it 
may be difficult to suppress couplings to motional 
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modes in the x and y directions. This will be trae if we 
carry out logic on a string of ions in a linear trap and use 
the axial COM mode. The laser beams must be at an 
angle with respect to the z axis to be able to select 
particular ions. This kind of problem can be eliminated 
by using the x ot y COM mode for logic and making k 
orthogonal to the other directions. Finally, we remark 
that the use of standing waves will give rise to optical 
dipole potentials which can shift the ion motional fre- 
quencies. These shifts should be small, and they can be 
incorporated into the definitions of the motional oscilla- 
tion frequencies. 

4.4.6.4 Spontaneous Emission 

The internal atomic states of trapped ions, which 
store quantum bits of information, must be protected 
from spontaneous emission, at least for the duration of 
the computation. This excludes the possibility of "error 
correction," (See Sec. 3.3) which can tolerate a certain 
level of errors due to spontaneous emission. As dis- 
cussed in Sec. 4.2.1, for qubit levels coupled by single 
photon optical transitions, this may be accomplished by 
employing long-lived energy levels which do not have an 
allowed electric dipole coupling, such as metastable 
electronic levels with a quadrupole or intercombination 
coupling to the ground state. Although the interaction of 
these states with the vacuum (causing spontaneous 
emission) is reduced, their interaction with an external 
field for use in quantum logic operations is also reduced. 
This results in a fundamental limit on the accuracy of 
each operation by roughly the ratio of the spontaneous 
emission rate to the Rabi frequency ^ = F/fi. In the case 
of optical transitions, fi cannot be increased indefi- 
nitely, since at optical intensities beyond about lO'"* W/ 
cm^, the atom is quickly photoionized. This amounts to 
inaccuracies due to spontaneous emission on the order 
of ^= 10"" - 10"' [289]. Even this limit may be too 
optimistic, as the two-level approximation breaks down 
before photoionization occurs, and the coupling to other 
electronic levels must also be considered [188,289]. 
This results in inaccuracies due to spontaneous emission 
on the order of ^= 10"' — 10"", depending on the par- 
ticular ion used. To understand the nature of the prob- 
lem, we examine a simplified system. More detailed 
treatments are given in Refs. [188] and [289]. 

Consider the situation shown in Fig. 9. We assume 
that levels I -l) and | T) comprise the qubit states. We 
want to drive coherent transitions between these two 
levels but, in general, we must consider spontaneous 
Raman scattering from other non resonant spectator 
levels. For this approximate treatment, we will consider 
only the effects of the nearest, most strongly coupled 



state which is designated as state U) in the figure. For 
simplicity, we assume that levels I s) and I T) are cou- 
pled only to the ground state I -l) (and not each other) by 
the applied radiation, and that they decay by the same 
coupling process to the ground state with rates F and 
Ft- To allow for different coupling processes for the I s) 
and I T) states, we let Ft = kF^ and we neglect w^ factors 
in the differences of the lifetimes (this would be valid if 
the 14-)^ IT) transition frequency is much larger than 
A). The worst case is given by considering sideband 
transitions between levels I -l) and I T). From Sec. 2.3.1, 
the resonant Rabi frequency is given by Hi — Tjilr.i — 
rjK^^fis^l. We assume that during the operation, the ion 
has a probability of about 0.5 of being in the | i) and 
I T) state. With these assumptions, the ratio ^ of sponta- 
neous emission (from both upper levels) to Rabi fre- 
quency is 



i-fiA'-^M^J-m{^--J- «*' 



where ^ = Of/irjA)^, and the factors of V2 come from the 
probabilities the ion is in the I 4-) or I T) state. The most 
optimistic answer is given by minimizing ^ with respect 
to K which leads to ^— FJirjA). As an approximate 
"best case" we take FJl-n = 25 MHz, A/l-n = 10"* Hz, 
njl-n = 10 MHz, 7? = 0.1, giving ^=- 2.5 X 10"" and 
requiring k — 0.4^ and Ft — 10""/^. Therefore, weakly 
allowed transitions are desirable if single photon optical 
transitions are used for qubit transitions. 

In the case of two-photon stimulated-Raman transi- 
tions between stable electronic ground states, the ratio of 
spontaneous emission rate to Rabi frequency is approxi- 
mately giR = yJigVAR), where y^ ~ FgVAl is the off- 
resonant spontaneous emission rate, g is the resonant 
single-photon Rabi frequency of each laser beam, and 
Ar is the detuning of the Raman beams from the excited 
state (Sec. 2.3.3). This results in an inaccuracy F/Ar due 
to spontaneous emission, which is independent of opti- 
cal intensity. Since Raman transitions between S elec- 
tronic ground states are effective only when the detun- 
ing Ar is not much greater than the fine structure 
splitting of the atom [290] ; this results in an inaccuracy 
^SR due to spontaneous emission in range from about 
10""* ('Be-') to 10"' ('"Hg^), depending on the particular 
ion used. Spontaneous emission from spectator elec- 
tronic levels should not significantly affect this limit, 
provided that their splitting from the virtual excited 
state significantly exceeds Ar and that the single photon 
resonant Rabi frequencies coupling the ground states to 
the spectator levels are not much bigger than g [188]. 
These appear to be reasonable assumptions for most 
candidate ions. 
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The decohering effects of spontaneous emission can 
be overcome by error correction schemes. Error correc- 
tion is complicated by the fact that when spontaneous 
emission occurs, the atoms may decay to states outside 
the original set of computational basis states. However, 
this situation can, in principle, be detected by optically 
pumping the ions back to the computational basis and 
applying the error correction schemes [291,292]. 

Spontaneous emission decoherence could, in princi- 
ple, be nearly eliminated by driving single-photon tran- 
sitions between ground-state-hyperfine or Zeeman lev- 
els with rf or microwave radiation since spontaneous 
emission from these levels is negligible. Refs. [3] and 
[293] discuss this possibility, where inhomogeneous 
magnetic fields couple the internal and motional states. 
The speed of sideband operations is limited by the size 
of the field inhomogeneity one can obtain. From Eq. 
(28), we want a coupling Hamiltonian Hi= — jXxidBj,/ 
dz)z- Assume dBj,/dz = idBj,/dz)oCosa)t, ^^^ = ^(-M Sx, 
and z = Zo(a + a^) as in Sec. 2.3.1. For resonance at the 
first red sideband (o)= Wq — wj. Hi is given by Eq. (27) 
with Tqfl = — fjuyiidB^I dz)oZol{^fi)- If we take fjuM = AIb, 
Zo = 10 nm, then to achieve OJIt; = 1 MHz, we require 
idBJdz)o == 290 T cm"'; a difficult task. Moreover, it 
would be difficult to address selected ions because of 
the long wavelength of the radiation relative to typical 
ion spacings (see also Sec. 5.3). 

4.4.7 Mode Cross Coupling During Logic 
Operations 

In the preceding sections, we have assumed that when 
transitions are driven between I l^j I ««:) and I T)^ | ni,) 
involving a single mode of motion k, the other 3L-1 
spectator modes of motion are not affected because 
coupling to them is nonresonant. However, when the 
sum or difference frequency of two or more spectator 
modes is near the frequency of the desired mode-A: tran- 
sition (coi, I «t ~ ««: I ), higher order couplings can entan- 
gle the I -l)} I rik) and I T)^ I ni,) states with the spectator 
mode states. 

Equation (33) describes the general interaction 
Hamiltonian between the internal levels of ion j and all 
3L modes of motion. By expanding the exponential in 
Eq. (33) to all orders, we find 



H'ii = An' 



sJU 1 



1=1 b,,df=0 



b,\d,\ 



vibrational quanta, accompanied by the raising or lower- 
ing of the internal electronic levels of ion 7. In general, 
we must account for all terms in Eq. (139) which do not 
vary rapidly in time, or terms in which the resonance 
condition is nearly met: Xi(bi — di)u)i— 8=0)^(11^' — 
til,). Although detailed treatment of this problem is be- 
yond our intent, a couple of comments may be made. 

In general, we must account for all the terms in Eq. 
(139) which cause significant errors in the overall com- 
putation we are trying to carry out. For two or more 
trapped ions, some combination of modes will nearly 
always satisfy the resonance condition. However, this 
may occur only for high orders of bi and di , and if the 
Lamb-Dicke criterion is met, the contributions are van- 
ishingly small. The terms that will cause problems are 
the ones that are close to satisfying the resonance condi- 
tion and are relatively low order in bi and di. If the 
Lamb-Dicke criterion is satisfied, it will always be pos- 
sible to avoid these spurious couplings, but it may be at 
the expense of making the Rabi frequency so small (in 
order to avoid coupling to relatively nearby off -resonant 
terms) that the operations will become too slow. 

To understand this problem in the context of a simple 
example, we assume that a cross-mode coupling of this 
type occurs when two modes, p and q, have frequencies 
which satisfy the condition n^Wp — n^Wg — 0, or ± Wj. 
corresponding to possible extraneous mode coupling on 
the carrier, first blue sideband, or first red sideband of 
the logic operations (assumed to utilize mode k). This 
additional resonance condition yields, to lowest order in 
the Lamb-Dicke parameters, the effective Hamiltonian 

Hy' = AO'sJ 1 + iriiatc'""'' + ale""') 



_!_ (iVpa})''''(iVgagT'' ^Hn, 



Up !«g ! 
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where the resonance conditions are given by S= w — 
o)o = 0, + (Ok, OT — u)t. A specific example is relevant to 
the NIST single 'Be* ion experiments. Here, mode k was 
the X oscillation, and modes p and q are identified with 
the z and y oscillations of the single ion in the trap. In 
this experiment rij, = k ■ xxq, riy = k ■ yyo, ri^ = k ■ zzo, 
and o)j, — o)^ — <Wy . (The frequency relationship 
ft)^ = (u^ — (U3, is a consequence of Maxwell's equations 
for a quadrupole rf trap in the absence of static poten- 
tials applied to the electrodes [211].) We assume that the 
desired transition is the first blue sideband of mode 
x(S= (Oj,). In this case, the resonant part of Eq. (140) 
becomes 



where £2' = /2exp[— ^hSiirjiy]. This equation describes 
the processes of each mode / gaining or losing (bi — d{) 



Hi =^ Ail' 



S+\ irj^al — rj^rjyalay + Oirf) \ + h.c. 



(141) 
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The term proportional to aj is the desired anti-Jaynes- 
Cummings operator, and the term proportional to a^Oy 
can entangle the internal state with the other spectator 
modes (z and y), leading to errors. 

For logic operations on a string of ions in a linear 
trap, we will assume that all other mode frequencies are 
higher. With the use of stimulated-Raman transitions, 
we can make AA:||2 and restrict our attention to spectator 
modes along the z axis. Nevertheless, as L becomes 
large, nearby resonances of the type shown in Eq. (140) 
will become harder to avoid. These coupling terms al- 
ways scale as products of Lamb-Dicke parameters. Thus 
if the spectator mode Lamb-Dicke parameters are small 
enough, or if at least one Lamb-Dicke parameter is ap- 
proximately zero, the higher order unwanted resonances 
may be sufficiently suppressed. Furthermore, if the 
spectator modes are cooled to near the zero-point energy 
((n) « 1), then any couplings in Eq. (140) with powers 
of the annihilation operator a, will be absent most of the 
time. Hence, in large registers, it will probably be im- 
portant to cool all modes to near the zero-point energy. 



5. Variations 

5.1 Few-Ion Accumulators 

Many of the problems anticipated above, such as un- 
wanted coupling to adjacent ions or spectator modes, 
will be aggravated by a large number of ion qubits in a 
single ion trap "register." Moreover, implementation of 
quantum error correction schemes will require highly 
parallel processing [294]. Therefore, a multiplexing 
scheme for ion qubit registers is desirable. One possibil- 
ity is to perform all logic in minimal accumulators 
which hold one or two ions at a time [65] . Ions would be 
shuffled around in a "super-register" and into and out of 
(multiple) accumulators which are well shielded from 
the other ions. The shuffling could be accomplished 
with interconnected linear traps with segmented elec- 
trodes as shown schematically in Fig. 10; construction 
of such traps appears to be possible with the use of 
lithographic techniques [295]. Single-bit rotations on 
the mth ion would be accomplished by moving that ion 
into an accumulator. Logic operations between ions m 
and k would be accomplished by first moving these ions 
into an accumulator. An accumulator could also hold a 
second species of ion (say Mg*) which could be used to 
provide laser cooling to the | n = 0) state of the motional 
mode used for the logic operations, if necessary. There- 
fore, for logic operations, an accumulator would hold 
two computational ions and the auxiliary ion. This 
scheme should make it easier to select ions with laser 
beams because it should be possible to address one ion 
while nulling the laser intensity on the other ion, even 



with very high trap frequencies (see Sec. 4.4.4). The 
very small number of logic ions in an accumulator (1 or 
2) would make extraneous mode coupling much easier 
to avoid. The main disadvantage appears to be that com- 
putational speed is reduced because of the time required 
to shuffle ions in and out of the accumulator and provide 
laser cooling with the auxiliary ion, if required. How- 
ever, energy shifts of the ion's internal structure, due to 
the electric fields required to move the ion, need not be 
severe. For example, to move a 'Be"^ from rest to a 
location 1 cm away (and back to rest) in 1 (xs would 
require a field of less than 50 V/cm. Electric fields of 
this order should give negligible phase shifts in qubits 
based on hyperfine structure (Sec. 4.2.3). The accompa- 
nying phase shift caused by time dilation for a transition 
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Fig. 10. Schematic diagram of an ion trap "super-register" containing 
few-ion accumulators in which logic operations are performed. We 
assume the trap is generically equivalent to the linear trap of Fig. 1. 
The rods of the trap in Fig. 1 are replaced by segmented electrodes in 
which adjacent segments are at the same rf potential but where the 
segments support different static potentials. This allows ions to be 
selectively moved around an in and out of ion accumulators. In the 
accumulator shown a logic operation between ions k and m is indi- 
cated. The configuration of the laser beams is chosen to null the 
intensity on ion m while performing an operation on ion k. The beams 
must be reconfigured to perform an operation on ion m while leaving 
ion k imaffected. The auxiliary ion may be required to perform laser 
cooling each time new logic ions are moved into the accumulator. 
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frequency coo/2ti = 1.25 GHz (hyperfine structure in 
'Be*) would be less than 1 [jurad. 



5.2 Multiplexing With Internal States 

In principle, it should be possible to multiplex quan- 
tum information into all 3L modes of ion motion. This 
will probably be useful in experiments which use a 
small number of ions. However, because of the higher 
potential for decoherence of motional modes, it might 
prove more useful to multiplex quantum information in 
multiple internal states. This would give more work 
space with a smaller number of ions. Quantum logic 
within multiple internal states would most likely follow 
the ideas of NMR quantum computing [165,166]. A 
simple example is given by the eight ground state hyper- 
fine levels in 'Be* which can be labeled like states of 
three coupled spin-V2 particles I €,, €2, €3) where 
e, G {0,1 }. For instance, a Toffoli gate [140] is realized 
by driving a tt pulse on the \f,Mf) = 1 2,2) -^ 1 1,1) 
transition in the ground state hyperfine manifold of 'Be* 
ifwe label the states 1 2,2)= 1 1,1,0) and 1 1,1)= 1 1,1,1), 
in the notation \f,Mf)= I fi, €2, €3). Technically, the 
readout of specific internal quantum states would be 
more complicated since the methods employed in the 
experiments so far distinguish between only two possi- 
bilities, seeing fluorescence or not seeing fluorescence, 
on any given ion. However, by appropriate mappings, 
superpositions of any two internal states could be 
mapped onto two states of "readout" ions and the detec- 
tion accomplished as described above. To multiplex in- 
formation into internal states, we need a way to map 
qubits from any ion into the internal levels of other ions. 
To give a specific example of how this can be accom- 
plished, we consider the case of 'Be* where we might 
wish to multiplex into the \f,Mf)= |2,2), |1,1), |2,1), 
and I 1 ,0) states shown in Fig. 1 1 . This figure is essen- 
tially the same as Fig. 5, but for simplicity of notation, 
we label these hyperfine states as 1 0), 1 1), 1 0'), and 1 1 ') 
respectively. We want to accomplish the transformation 

^,- = «|0)A|0)B+/3|0)A|l)B+r|l)A|0)B+8|l)A|l)B 



«|0)A+/3|l)A+r|0')A+8|l% |0) 



\9Mr] 



(142) 

between ions A and B. This can be realized with the 
transformation 



M — Ca'.B ■ Ca,B ■ Cb,A' ■ Cb,A ■ HA,l,0't 



(143) 



Internal state 
energy 



— — 1' fi-U 



(2,-2) 



(2,-1) 



(2,0) 



(2,2) „ 



where the operator IIa.i.o' denotes a it transition between 
the states |1)a and |0')a on ion A and the operators Cc,t 



Fig. 1 1. Hyperfine levels of the 2s ^S 1/2 ground state of 'Be* in a weak 
magnetic field (not to scale). The energy levels are designated by 
horizontal lines. Above the lines, the levels are represented by atomic 
physics labels (F,Mf) where F is the total angular momentum and Mf 
is the projection of the angular momentum along the magnetic field 
axis. The separation of Zeeman substates in the different F manifolds 
is approximately equal to 0.7 X lO'^Bo Hz where Bo is expressed in 
teslas. The separation of the F = 1 and F = 2 manifolds is approxi- 
mately 1.25 GHz at Bo = 0. For simplicity of notation, in Sec. 5.2, we 
make the identifications \F,Mf)= 1 2,2)= |0), 1 1,1)= 1 1), 1 2,1) = 

lO'), |i,o)^|r). 



are defined in Sec. 3.3. For example, Ca;b denotes a 
controlled-not operation with the |0')a and |r)A states 
acting as the control bit states and the |0)b and | 1)b states 
acting as the target bit states. 

5.3. High-Z Hyperfine Transitions 

The possibility of using highly charged ions as qubits 
poses some interesting possibilities. For simplicity, we 
consider using hyperfine states in high-Z hydrogenic 
ions as qubit levels where Z is the nuclear charge. 
More complicated high-Z ions appear to show similar 
features. If Zis high enough, hyperfine transitions occur 
in the optical region of the spectrum. Therefore, the 
required high field gradients necessary for quantum 
logic would be provided by the laser fields used to drive 
the hyperfine transitions. The potential advantages 
are: (1) Off- resonant spontaneous scattering, which is a 
source of decoherence in other cases (Refs. [188], [289], 
and Sec. 4.4.6.4), is essentially eliminated since the first 
excited 2p level is at much higher energy above the 
ground state. (2) The trap binding of high-Z hydrogenic 
ions can typically be stronger that singly-ionized atoms, 
thereby increasing motional frequencies and reducing 
the time of fundamental operations. The disadvantages 
are (1) High-Z ions are hard to produce. (2) The life- 
times of the upper hyperfine states are short enough that 
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spontaneous emission decoherence from these states 
cannot be neglected. (3) Detection of internal states may 
be difficult. 

The ion production problem is not intractable; for 
example, recent storage ring measurements have been 
able to determine the transition wavelengths and life- 
times of hydrogenic ^093^82+ |^296] and '"'Ho""^ [297]. 
Effort in some laboratories is directed to transferring 
similar high-Z ions to low energy ion traps. The hyper- 
fine frequency of a hydrogenic ion can be estimated to 
be Whft/2TT =2.54X 10** Z^F |^„| Hz, where F and F-1 are 
the values of angular momentum in the hyperfine states 
and g„ is the nuclear g -factor {^x. = iijgn, (x„ = nuclear 
magneton, /= nuclear spin) [298]. The spontaneous 
decay rate from the excited state is estimated to be 
Trad = 1/T,,d- 1.01X10-"^ [MJ{2F,+ \)] («hfs/2T7)^ s-', 
where the term in square brackets [298] depends on 
matrix elements and is on the order of 1. For example, 
the hyperfine transition in ^"'Bi*^* has a wavelength of 
244 nm and the upper level decays with a lifetime of 
351 |xs [296]. For "^'Ho'^'^'', the transition wavelength is 
573 nm and the decay time is 2 to 3 ms [297]. These 
experimental numbers agree roughly with the formulas 
above. These wavelengths are reasonable for quantum 
logic, but the lifetimes are somewhat short. However 
since the lifetime scales as 1/Z', it might be possible to 
use lower Z, longer wavelength transitions. The readout 
of the qubits appears difficult, because coupling to a 
higher electronic level and observing fluorescence is 
impractical due to the short wavelengths involved. In 
principle, however, the qubit information from these 
ions could be transferred to other ions whose states are 
easier to read, perhaps using small accumulators as 
described in the previous section and in Sec. 4.3. 

6. Other Applications 
6.1. Quantum Correlations 

In a classic Paper, Einstein, Podolski, and Rosen 
(EPR)[299] argued that quantum mechanics provides an 
incomplete description of physical reality and specu- 
lated that it might be an approximation to some theory 
which would provide a complete description. Today 
such theories, which supplement quantum mechanics 
with additional, as-yet unobserved parameters, are 
called hidden-variables theories. If such a theory obeys 
some physically reasonable conditions forbidding action 
at a distance, it is called a "local hidden-variables 
theory." 

Simple quantum logic gates performed on small 
numbers of trapped ions can lead to interesting experi- 
ments which may shed light on the viability of local 
hidden-variables theories. For example, as described by 
Cirac and Zoller [1] and Sec. 3.4, using controlled-not 
gates, we can generate the state 



\fr- 



V2 



(|i)i|i).+e'^|T),|t).). 



(144) 



In the context of EPR, this is called a Bell state. Accord- 
ing to Bell's theorem [300], such a state allows us to 
distinguish between quantum mechanics and all local 
hidden- variables theories. If the states of the two en- 
tangled particles are detected outside of each other's 
light cones, then, for particular sets of measurements, 
we may derive inequalities which all local hidden- 
variables theories must obey, but which quantum 
mechanics violates. The experiments performed by 
Clauser and Shimony [301] and Aspect et al. [302] 
provide strong evidence against local hidden-variables 
theories. Their work used polarization measurements on 
entangled pairs of photons. In their experiments, the 
detection of the photons' polarization states occurred 
outside each others' light cones. Thus, the measurement 
on one photon could not possibly have affected the other 
measurement, which closed possible "loopholes" in the 
proof of quantum mechanics over other explanations. 

However, some loopholes still remain open. Since the 
photon detection in the Aspect, et al. experiments was 
not 100 % efficient, the group had to make assumptions 
that the photons they measured were a "fair" sample of 
the whole population of events. Thus, their experiments 
do not rule out the (seemingly implausible) possibility 
of local hidden-variables theories in which the hidden 
variables cause some sub-ensemble of the photon 
pairs to preferentially interact with the measurement 
apparatus. 

In the system of two ions, we may detect the state of 
either ion with nearly 100 % efficiency through the use 
of "electron shelving" (Sec. 2.2.1). On the contrary, it 
will be difficult to perform measurements on two ions 
outside each other's light cone. Such a measurement 
would require separating the ions by a distance larger 
than the speed of light times the measurement time or 
transferring quantum information over large distances 
[151]. (In principle, the ions could be first entangled and 
then placed in different traps which could be separated 
by large distances before measurements were per- 
formed.) Nonetheless, an experiment with two trapped 
ions could be viewed as complementary: the photon 
experiments definitively close loopholes of causality, 
and the ion experiments could close loopholes due 
to detection inefficiency. Such experiments have the 
additional appeal of studying EPR on massive particles 
[303]. EPR states of atoms have recently been created 
in an atomic beam using the methods of cavity QED 
[304]; if detection efficiency can be improved, these 
experiments could also close loopholes due to detection 
inefficiency. Moreover, even though measurements of 
quantum correlations between entangled ions cannot be 
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easily performed outside each other's light cone, one 
can argue strongly that the ions cannot transfer infor- 
mation by any known mechanism. Therefore, if the 
observed correlations violate Bell's inequalities, the 
correlations are established by some new force of nature 
or are, in fact, inherent in the structure of quantum 
mechanics. 

An intriguing possibility for ions is the possibility of 
making "GHZ states" [54,305,306]. For three ions, the 
GHZ state is 



\p- 



1 



V2 



(|i),|i)2|i)3+e-lT),|T).|T)3). 



(145) 



For such a state, a single measurement can distinguish 
between the predictions of quantum mechanics and 
those of any local hidden- variables theory [54,305]. 

Aside from this. Bell states, GFIZ states, and 
Schrodinger-cat states are highly entangled, and are thus 
of inherent interest for the study of uniquely quantum 
behavior. As the experiments improve, it will be interest- 
ing to push the size of entangled states to be as large as 
possible. The question is not whether we can make states 
which have the attributes of Schrodinger cats, but how 
big can we make the cats? Certain theories which 
address the measurement problem will be amenable to 
experimental tests, for example, quantitative limits on 
spontaneous wavefunction collapse theories [307,308] 
can be established. The isolation from the environment 
exhibited by trapped ions, coupled with the control 
possible over their quantum state and high detection 
efficiency make them an interesting laboratory for the 
study of fundamental issues in quantum mechanics. 

6.2. Simulations 

The nonlinearities with respect to motional raising 
and lowering operators inherent in the coupling 
Hamiltonians of Eqs. (14), (37) and (48) can lead to a 
rich variety of dynamics between the ion motion and 
internal levels. Some of these dynamics for a single ion 
have been discussed in Sec. 3.2. An interesting system 
which can be simulated with these couplings is a 
"phonon maser" which provides vibrational ampli- 
fication by stimulated emission of radiation [309]. This 
would be an analog of the micromaser with some 
interesting differences such as the effects of recoil in the 
pumping process. 

The Hamiltonians of Eqs. (14), (37) and (48) for a 
single ion can be consolidated into the general interac- 
tion Hamiltonian 



where 8 = w-Wo, where we have used the rotating-wave 
approximation, and where we have neglected phase 
factors of the field. The parameter e is equal to 1 when 
internal state transitions are involved and 6 = when the 
internal state is unchanged. In Eq. (146), k is the 
wavevector of the field for single photon transitions or 
k — > fci - A:2 when stimulated-Raman transitions are 
used. Similarly, w is the frequency of the applied field 
for single-photon transitions and w — > Wli - (^li when 
stimulated-Raman transitions are used. In an interaction 
picture of the ion's motion, this Hamiltonian becomes 



H[=/iQ 



5;e--(*'-*n exp(i[r7,(fl,.e-" 



+a] e.'''J')+h.c. 



(147) 



where ri^ = k-xxo, a^ is the lowering operator for the 
X motion of frequency (a^ , etc. Now, assume that H is 
small enough, and that, in general, (a^, u)y, and u)^ 
are incommensurate, so that we can excite only one 
spectral component of the possible transitions induced 
by this interaction. For a particular resonance condition 
8 = -i^u)^--(!yU)y-i^u)^, and in the Lamb-Dicke limit, 
we find 



H[^/iQ&'>'si n 
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+h.c. 



(148) 



the two mode case where €= £^ = is considered by 
Drobny and Hladky [310], and in a different excitation 
scheme in Ref. [110]. If the Lamb-Dicke limit is not 
rigorously satisfied, we must consider higher order non- 
linear corrections to this expression; specific examples 
are discussed in Refs. [28] and [110]. These nonlinear 
terms are the origin of the high-order corrections to the 
Rabi frequencies (Eq. (18)). The case of carrier, and first 
red and blue sidebands on internal state transitions (e.g., 
6=1, £^= iy= 0, /z = 0, ±1) are used extensively for 
quantum logic and generation of nonclassical motional 
states and are discussed above. The case 6 = 0, 



t y t 17 \J. V 7 



1 has been used to create coherent 



Hi = ;l/25'e*-"-*"'*i+ft.c. 



(146) 



[21] and Schrodinger cat [47] states of motion and is 
discussed in Sec. 3.1. The case 6 = 0, /j: = ^3, = 0, 
|/z| = 2 has been used to create squeezed states; this is 
discussed in Ref. [21]. A realization of the Hamiltonian 
H{'^S^{aJf+h.c. (6=1, /,= /3=0, /z = -2) has 
been reported by Leibfried, et al. [132]. This is similar 
to the case of two-photon excitation in cavity QED 
analyzed by Buck and Sukumar [311] and Knight [312]. 
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Interactions proportional to S+a^ay+h.c. , S+a^a^+h.c. , 
and a^a^+h.c. might be used to generate the maximally 
entangled state of Eq. (60) without the need to address 
individual ions [194]. An example of an interesting new 
case would perhaps be the realization of three-phonon 
downconversion (e.g., e = 0, /^ = 3, Z^, = - 1, /^ = 0). 
This case is analogous to three-photon downconversion 
in quantum optics (see Refs. [1 10],[313], and references 
therein). Here, it corresponds to driving a two-mode 
resonance using stimulated-Raman transitions where 
&)li - &)l2 = «j - 3&)j:. a suggestion to realize a 
Hamiltonian o^a^a^+kc. (e = 0, 4 = 0, /^ = 2, €y = -1) 
is discussed by Agarwal and Banerji [314]. Clearly, a 
large number of possibilities could, in principle, be 
realized just for a single ion; moreover, the number of 
possibilities increases dramatically if we consider all 
modes of motion for multiple trapped ions. The only 
limitation on how high \fj\ in Eq. (148) can be is that 
n be chosen sufficiently small that couplings to other 
(unwanted) resonances are avoided. This will require 
that decoherence be small enough to see the desired 
dynamical behavior before coherence is lost. 

Various forms of interactions which satisfy the 
requirements of quantum nondemolition (QND) mea- 
surements [20,25,315] of ion motion or quantum 
feedback [316] can be extracted from Eq. (148). QND 
experiments employing dispersive interactions are 
considered by Retamal and Zagury [36]. These schemes 
rely on measuring the ion's internal state however, 
which almost always involves recoil heating, thereby 
destroying the state we wish to preserve. This could be 
circumvented by coupling the ion's motion to a cavity 
field which then serves as the probe [317]. 



6.2.1 Mach-Zehnder Boson Interferometer With 
Entangled States 

Realization of the various Hamiltonians indicated in 
Eq. (148) can lead to simulation of various devices of 
practical interest. As an example, consider a Mach- 
Zehnder interferometer which acts on two modes of 
oscillation of a single trapped ion; to be specific, we will 
consider the x and y modes of motion. The analogy with 
a Mach-Zehnder interferometer for bosons is that 
the two input modes to the boson interferometer are 
replaced by the x and y modes of ion oscillation. The 
(50/50) beamsplitters in the boson interferometer are 
replaced by an operator B± = exp[±iTi{a^ay+aj,ay/4] 
[193,318,319]. This operator can be realized by 
applying the interaction in Eq. (148) with e=0, and 
^x = —Sy = 1 for a time given by flr\„r\yt = ttM. A 
differential phase shift between the two arms of the 



interferometer can be simulated by shifting the relative 
phases of the fields in Eq. (148) between successive of 
applications B±. In a particle (e.g., boson) interfero- 
meter, one typically measures the number of particles in 
either one or both output modes. In a single ion experi- 
ment, we have only one convenient observable, the inter- 
nal state of the ion (either 14-) or |T). Nevertheless, we 
can characterize the action of the phonon interferometer 
by repeating the experiment many times and measuring 
the density matrix of the output state [132]. 

It will be most interesting to characterize the action of 
the interferometer for various nonclassical input states. 
One example is the two-mode Fock state |nj:)j:|M},)}, 
[195]. This state could be prepared by applying the 
techniques described in Sec. 3.2 sequentially to the ion's 
X and y modes. This state is interesting because it has 
been shown that one could approach the Heisenberg 
uncertainty limit in a Mach-Zehnder interferometer by 
measuring the distribution of bosons in the output 
modes [195,196,197]. An alternative technique for 
studying the action of a beamsplitter on the two-mode 
Fock states has been suggested by Gou and Knight [23] 
when C0j,= cOy. Here, a beamsplitter could be simulated 
by first preparing |mj:)j: !«},)}, along two orthogonal axes 
and then probing along two other axes (x' and y') which 
are rotated (in the xy plane) with respect to the first. This 
technique could also be used to analyze, for example, the 
(|0);,.|2)3,. + |2);,.|0)y)/V2 state from an initially pre- 
pared |1)j:|1)j, state [23]. 

Another interesting state to consider for the phonon 
interferometer is the (\N)^\0)y+(\0)^\N)y)/2"^ state 
(which is the desired state after the first beam splitter). 
This state has been shown to yield exactly the 
Heisenberg uncertainty limit for an interferometer for 
any value of A^ [194], if after the second beamsplitter, 
we measure the number of particles N{x) in the x output 
port. The result of this measurement is assigned the 
value (-l)"*''l This state could be prepared from the 
|4')|0)j:|0)j, state by the following two steps: 

(1) Apply a 17/2 pulse on the A^th blue sideband of mode 
X [e = 1, 4 = -N, fj, = in Eq. (148)]; this creates the 

state (|i)|0),+ |T)|Ar),)|0),/2"l 

(2) Apply a tt pulse on the A^th blue sideband of mode 
y [e= 1, /, = fy=-N in Eq. (148)]; this creates the 

state |T)(|Ar),|0),+ |0),|Ar),)/2"^ 

After the second beamsplitter, we have a state which can 
be written as 



^fi„al=|T) 2 Cn^\n^)x\N-n^)y 



(149) 



314 



Volume 103, Number 3, May-June 1998 

Journal of Research of the National Institute of Standards and Technology 



In principle, we would like to measure it^, record the 
value N(x), and assign the value (-1)"*^'. Effectively, 
this assignment can be accomplished if we can find an 
interaction M which provides the mapping 

JV 



+ 1^) E C„,e'«"'>,),|A^-n,), 



(150) 



After this mapping, we need only measure the internal 
state; if the ion is found in the |T) state we assign the 
value +1 to the measurement; if the ion is found in 
the |4')state, we assign the value -1. The mapping 
M can be achieved by applying radiation with k \\x 
at the carrier frequency (€=1, €^ = 0) and insuring 
/3„,„f = 2TTm ± n;c IT where m is an integer. From Eq. (18), 
we have 



/3„,„ t~ Ote 2 



\-n^v}[ 1 + 



lil 



-ni 



(151) 



Therefore, if we make /3exp(- 17^/2) f = 2ottt and 
r}^(l + r}^/4 = (2m~^) we achieve the desired mapping as 
long as the contribution to the phase from the term 
proportional to n} in this equation is small compared to 
TT. Therefore we require m»N^I% or, equivalently. 

One final example of a two-mode phonon interfero- 
meter which directly yields Heisenberg 1/A^ phase sensi- 
tivity is a "beamsplitter" which creates the state 
(|i)|0),|A^),+|T)|A^),|0),)/V2. This state can be cre- 
ated by starting with an initial (|i)+|t))|A/^/2);,|A/^/2)3,/ 
V2 dual Pock state (Sec. 3.2). Next, we apply NI2 
TT-pulses alternating between the two interaction 
Hamiltonians Hi = nrj^rjyiS+a^ay+h.c.) and 

H2 = Orj^rjyiS^xay+h.c. ) which can be realized in the 
Lamb-Dicke regime. In this way, the state of motion of 
the ion is stepped through 



it)if.i).if-i),.|i)if-i)4.i), % 

|i)||,2)4-2)..|T)|f-2)4.2), % 



(152) 



|i)|A^),|0),+ |T)|0),|A^), 



where, in this example, we assume N/2 is even. The 
interactions Hi and H2 follow from Eq. (148) with e=l. 



f^ = -1, f^ = 1, /^ = and e=l, 4 = 1, f^ = -1, €, = 0, 
respectively. The ^th pulse has Rabi frequency 
Or],r]y[(N/2+k)(N/2-k+l)]"^ in the Lamb-Dicke re- 
gime. After a relative phase is accumulated in the two 
"paths" of the interferometer (simulated by adjusting 
the phase of the laser pulses as discussed above), we 
reverse the above steps and apply a final "tt/2" pulse on 
the carrier. Upon measuring the probability of occu- 
pation in state |4-) or |T), the interference fringes exhibit 
1/A^ phase sensitivity. 

If the Lamb-Dicke criterion is not satisfied, the two 
components of the wavefunction superposition may 
experience different Rabi frequencies during each pulse, 
leading to undesired evolution. The exact Rabi frequen- 
cies of the two components of the wavefunction upon 
application of the ^th pulse (interaction Hi or H2) follow 
fromEq. (18): 



/3a = 



/3r = 



nr,,r,y e-"^'"'^"^* U,2.*-i(t].') U/2-* iVy) 



\/(N/2+k)(N/2-k+l) 



(153) 



/JTj.Tj, e-'"'"-'^"/' a,2.UVy) Ul2-k iV^) 

\/(N/2+k)(N/2-k+l) 



where Oa is the Rabi frequency between states 
\N/2+k-l%\N/2-k+l)y and \NI2+k)JNI2-k)y and /3b 
is the Rabi frequency between states \NI2-k+\)x\NI 
2+k-\)y and \NI2-k\\NI2+k)y. The only differences 
between /3a and /3b are the arguments of the associated 
Laguerre polynomials. Thus, as long as 7}:,= rjy, the 
system will evolve as in Eq. (152), even when the Lamb- 
Dicke criterion is not satisfied. 



6.2.2 Squeezed-Spin States 

As an example which demonstrates the advantages of 
using entangled states in spectroscopy, we discuss the 
following simple experiment which can be carried out 
on a single ion. Let us suppose we are interested in 
measuring, with maximum signal-to-noise ratio, the 
Zeeman frequency between states of a / = 1 manifold in 
an atom. To be specific, consider that we are interested 
in measuring the Zeeman frequency of the AMp = ± 1 
transitions in the ^SiaiF = 1) ground-state hyperfine 
multiplet of a 'Be"^ ion (Figs. 5 and 1 1). This manifold 
is composed of the three levels \F,Mf) G. {|l,0,), 
1 1, ± 1 )}. We will assume that the applied field is small 
enough that the frequency of the 1 1,-1) ^ 1 1,0) transi- 
tion is equal to the frequency of the 1 1 ,0) — > 1 1 , 1 ) transi- 
tion. A straightforward way to measure the Zeeman 
frequency would be to prepare the atom in the |l,l) 
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(or 1 1,-1) state, drive the Zeeman transition (using the 
Ramsey method), and then measure J^ . We will assume 
this detection can be accomplished with 100 % detec- 
tion efficiency, so the noise in the measurement is 
limited by the quantum statistics in the measurement 
process [100]. After many measurements, requiring a 
total averaging time t, a certain precision in the Zeeman 
frequency would be obtained. However an alternative 
measurement strategy, described below, would require 
an averaging time t/2 to reach the same measurement 
precision. 

The basic idea is outlined in Ref . [9] ; for / = 1 , it is 
equivalent to the more general technique described by 
[194]. We prepare the atom in the state (/((O) = 1 1,0) and 
then apply the Ramsey fields. Subsequently, we measure 
the probability of the ion to remain in the |l,0) state. 
(This measurement scheme has been used in the exper- 
iments of Abdullah et al. [320]; the noise in these exper- 
iments, however, was not limited by projection noise.) 
This measurement is equivalent to measuring the opera- 
tor I- J} , where Jis the identity operator, and therefore 
equivalent to measuring the operator J^ , as discussed by 
Wineland, et al. [9], or the operator YliCr^i, as discussed 
by Bollinger, et al. [194]. After application of the first 
Ramsey pulse to the |l,0) state, the ion is in the state 
(1 1,1> + 1 1,-1))/ V 2; this is equivalent to the maximally 
entangled state of Eq. (60) which could be formed by 
two spin-1/2 particles. As described in Sec. 3.3, we 
would expect the time required to reach a certain 
measurement precision to be reduced by a factor of 
2 (L = 2) over the case of uncorrected particles, repre- 
sented by starting the ion in the state (/((O) = |l,l) or 
1 1 ,- 1 ) state and finally measuring J^ . 

We could carry out this experiment in the following 
way. We first optically pump the ion to the |2,2) hyper- 
fine state. We then apply two successive tt pulses which 
carry out the transformations |2,2) — > |2,1) — > 1 1,0) with 
stimulated-Raman transitions. We then perform Ramsey 
spectroscopy on the |l,0) state. We can measure the 
probability of subsequently finding the atom in the 1 1 ,0) 
state by first reversing the order of the two tt pulses and 
then measuring the probability of finding the ion in the 
1 2,2) state. 

6.3 Mass Spectroscopy and NMR at the Single 
Quantum Level 

Quantum logic operations may be useful in precision 
measurements other than spectroscopy (Sec. 3.4). For 
example, Ref. [48] discusses a method to measure 
cyclotron resonance frequencies of single ions at the 
quantum level; this technique essentially employs quan- 
tum logic to distinguish between motional quantum 
states. This capability would provide mass spectroscopy 



at such low energies that anharmonic shifts (electric 
potential and relativistic) would be very small. (At 
present, however, high-resolution mass spectrometry 
experiments are limited by other effects such as mag- 
netic field drifts [77,78,79].) In mass spectroscopy, the 
basic problem is to measure a "test" ion's cyclotron 
frequency and compare it to the cyclotron frequency of 
a "measurement" ion (in the same magnetic field) which 
serves as the transfer standard. 

One idea [48] is to store two ions in separate Penning 
traps which are stacked along their symmetry axis and 
share a common endcap and a common, axial magnetic 
field. The axial frequencies of the ions are adjusted to be 
the same, in which case the axial harmonic oscillators 
are coupled through the charge in the common endcap. 
The technique might work as follows: The axial motion 
of the "measurement" ion (for example, 'Be*) is first 
laser cooled to the n = level. When the ions' axial 
modes are resonantly coupled, this cooling can be trans- 
ferred to the test ion's axial motion; if the coupling is left 
on for a certain amount of time, the energies in the ions' 
axial oscillations are exchanged. Subsequently, the 'Be"^ 
ion's axial oscillation is recooled so that both ions are 
cooled to the zero-point energy. This axial cooling is 
then transferred to the test ion's cyclotron mode by 
parametric coupling, after which the axial motion of 
the measurement ion must be recooled to the zero-point 
energy. An external field is then applied to weakly 
excite the test ion's cyclotron motion. When the ampli- 
tude of this field is adjusted appropriately and the 
resonance condition is met, the n = to n = 1 transition 
in the test ion is driven with high probability. The steps 
above are then reversed so that if the test ion's cyclotron 
motion was excited to the n=\ level, the 'Be* ion's axial 
motion is now in the n=\ level. When the resonance 
condition for the test ion's cyclotron frequency is not 
met, the 'Be* ion's motion remains in the n = level. 
Discrimination between the n = and n=\ axial level is 
then performed using quantum logic operations applied 
to the 'Be* ion as discussed in Sec. 3. For example, after 
the steps above, an axial red sideband tt pulse will excite 
a I-1-) — > |T) transition in the 'Be* ion conditioned on 
whether or not the 'Be* ion's axial mode was in the 
M = or n = 1 state. The test ion's cyclotron frequency 
can be referenced to a spin flip frequency in 'Be* (either 
electron or nuclear) which then acts as a transfer 
standard. Finally, by performing a cyclotron resonance 
measurement on a second test ion in the same fashion 
we can find the ratio of the two test ions' cyclotron 
frequencies and therefore derive their mass ratio. 

Ref. [48] also suggested that these ideas could be 
applied to measure magnetic moments of test ions. In 
this section, we describe a variation on the method dis- 
cussed there. The basic idea is to perform NMR on an 
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unknown magnetic moment (the test ion) at the single 
spin level. Using quantum logic operations, the spin flip 
is detected in a measurement ion to which the test ion is 
coupled. We will describe the ideas in the context of a 
specific example; but the techniques are easily general- 
ized to measuring the magnetic moments of other ions. 

Suppose we want to measure the ratio of proton and 
antiproton magnetic moment. Such a measurement is of 
high current interest and can provide a second test of 
CPT on baryons in addition to the precise mass com- 
parisons already performed [79]. Precision measure- 
ments of the magnetic moment ratios coupled with the 
mass ratios should test much of the same physics in the 
context of CPT as precision spectroscopy of antihydro- 
gen [321]. A measurement scheme for the proton/ 
antiproton moment ratio has been suggested previously 
by Quint and Gabrielse [322]. This method would 
employ the same basic ideas as in the electron g-2 
experiments of Dehmelt and coworkers [108], where the 
proton's (or antiproton's) magnetic moment energy is 
transferred parametrically into its cyclotron energy 
(which is then detected) by use of applied inhomoge- 
neous (oscillating) magnetic fields. The apparent draw- 
back to this scheme is that it is slow because of the 
weakness of the parametric coupling and the difficulty 
of detecting small changes in the cyclotron energy. The 
method suggested here incorporates the same kind of 
parametric coupling but is potentially more efficient 
because it is sensitive to transfer at the single quantum 
level through the use of quantum logic techniques. 

The basic idea of the method we propose here is to 
first store a proton and a 'Be"^ ion in trapping potentials 
which are in close proximity and in a common mag- 
netic field (Fig. 12). The proton spin flip frequency is 
then compared to the 'Be* spin flip frequency (elec- 
tron or nuclear), which acts as a transfer standard, effec- 
tively calibrating the (common) magnetic field. We then 
perform the same kind of measurement on a simulta- 
neously trapped antiproton and 'Be"^ ion. By combining 
the two measurements, we determine the proton/ 
antiproton magnetic moment ratio. 

First consider the measurement of the proton to 'Be"^ 
spin flip frequency ratio. We simultaneously store a 
single proton and 'Be"^ ion in a double-well potential as 
indicated in Fig. 12. We assume that the ions are con- 
fined in a direction perpendicular to the z axis by a 
linear rf trap combined with a superimposed static 
magnetic field B = BoZ. (One particular geometry for a 
Paul trap in a strong magnetic fields is described by 
Bate, et al. [323].) The quadrupole electrodes are 
segmented to provide a static double-well potential 
along the z direction. We assume the electrode segments 
are capacitively coupled together so the rf electrode 
potentials are independent of z . A vertical wire (shown 
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Fig. 12. Schematic diagram of a trap for simultaneous storage of two 
ions with different charge to mass ratios, q\lm\ and q2/m2. We assume 
the trap is generically equivalent to the linear trap of Fig. 1 . The rods 
of the trap in Fig. 1 are replaced by segmented electrodes in which 
adjacent segments are at the same rf potential but where the segments 
support different static potentials in order to make a double well 
potential along z. Therefore in this figure, the view is from above, 
where we see two (of four) segmented electrodes which have replaced 
two of the rods in the trap of Fig. 1 . On the left is a cross section of 
a wire through which an rf oscillating current is sent; this wire gener- 
ates an rf magnetic field at the site of the ions as explained in the text. 



in cross section in the figure) is at one end of the trap; 
oscillating currents in this wire generate an oscillating 
magnetic field which can be used to drive the spin flip 
transitions of the proton and nuclear spin flip transition 
in the 'Be* ion (the electron spin flip transition in 'Be* 
could be driven with injected microwave radiation). 
Oscillating currents in this wire will also provide a 
parametric coupling between the proton spin and axial 
motion (below). 

The experiment could proceed as follows: We first 
tune the static potentials so that the proton and 'Be* 
axial frequencies are the same. By performing Doppler 
laser cooling on the 'Be* ion's axial and radial modes, 
the proton's axial motion is cooled. The proton's radial 
modes will also likely be cooled through the Coulomb 
coupling. If not, we must apply an inhomogeneous 
rf field which parametrically couples the proton's radial 
and axial frequencies. We now uncouple the proton 
and 'Be* ion by adjusting their axial frequencies to 
be different. We then cool the 'Be* axial motion to the 
ground state (Sec. 3.1) and follow this by switching 
the axial modes back into resonance. The Coulomb 
interaction between the ions gives rise to a coupling 
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?i^2/(4'n'eo|Z2-Zi|), where Zi and Z^ are the respective 
axial positions of the ions and q\ and ^2 are, here, the 
charges of the proton and 'Be* ion. We have neglected 
the effects of induced charges in the trap electrodes, this 
will not be a large effect if the distance between ions is 
on the order of or less than the distance of each ion to 
the trap electrodes. For small amplitudes of oscillation 
we can write Z=Z,o+z,, !G{1,2} where Z,o are the 
equilibrium positions of each ion. In the limit of weak 
coupling and to lowest order in Zj, the interaction 
between ions is given by 



H„ 



q\q2Z\Z2 

271 €od^ 



(154) 



where d= |Zio-Z2o|. This gives rise to a frequency split- 
ting between axial normal modes of [48] 



8&)i_2 = ir/iexch = 



qiqi 



2'T!€od^ w^\m\m2 



(155) 



where m, are the ion masses and 4xch is the time to 
exchange the energy of the two ions' axial energies after 
they are coupled together. At the time 4xch, the ions are 
again decoupled, leaving the proton in its axial ground 
state. The 'Be"^ ion is subsequently recooled to its axial 
ground state. 

Assume the proton is initially in the lower energy 
state IT). We now apply an oscillating magnetic field at 
a frequency &>„ which is near the proton spin flip 
frequency Wo = g^ii^Bolfi where g^ is the proton gyro- 
magnetic ratio expressed in terms of the Bohr magneton 
(gp ~ 1.521X10"^). The amplitude and duration of this 
field is adjusted to make a it pulse if the resonance 
condition &>„ = &)o is satisfied; that is, the proton spin 
undergoes the transition |T) ^ |i). We now apply an 
inhomogeneous magnetic field at a frequency w^ by 
driving a current I cos co J through the vertical wire, 
shown on the left in Fig. 12. If we neglect the shielding 
effects of the trap electrodes, this provides a field at the 
site of ion 1 equal to B(t) = x ixqI cos co J/ (2ti p) and a 
field gradient dBJdz = -I cos co J/ (2ti p^) where pis the 
distance from the wire to the proton. This oscillating 
field gives rise to a coupling which has the form of 
Eq. (28). If the resonance condition coa =coo-co^ is satis- 
fied, this coupling takes the form of Eq. (27) with 
Oi = gp/aB/ao/zo/(4'T7^p^) where zo is the zero-point 
motion for the proton (as in Sec. 2) . We assume that the 
duration ^coupling of this inhomogeneous field is adjusted 
to give a complete transfer |~l)|0) -^ |t)|l)(/3ircoupiing = 
17/2). We now couple the axial modes for a time i^xch so 
that if the proton was in the |m = 1) state after the last 
step, the 'Be* ion is now in the |« = 1 ) state. Finally, the 
|n = 1) motional state is detected on the 'Be* ion using 



the methods as discussed above and in Sec. 3. If co„ is 
nonresonant, the modes remain in the |n = 0) state. 

If the proton is initially in the \-l) state and co„ is 
resonant, then the field at frequency cOa =coo-co^ has no 
effect, and the ion remains in the |n = 0) motional state 
and gives a false "no signal." However, when &>„ is 
reapplied, the |m = 1) state is produced giving a signal. 
After reinitializing (that is, preparing both ions in their 
axial ground states), signals are always produced if &>„ 
remains resonant. If &>„ is nonresonant and the proton 
initial state is |T), no signals are ever produced. However 
if &)„ is nonresonant and the proton initial state is \-l), we 
produce a false signal. After reinitializing and repeating 
the experiment, we produce no signals. This is true on 
subsequent tries as long as &>„ remains nonresonant. 
Therefore to reduce the effects of false signals, we 
should repeat each try several times for each value of &>„ 
and discard the first measurement. Interlaced with 
measurements of the proton spin flip frequency, we 
measure the 'Be* spin flip frequency by driving with 
resonant rf radiation and using the techniques outlined 
in Sec. 3 for detection. This allows us to monitor and 
correct for magnetic field drifts on a fairly short time 
scale. 

The antiproton/'Be* comparison is accomplished 
similarly except for the important difference that axial 
potentials must be trapping for one species and nontrap- 
ping for the other at a given location. Apparently, the 
biggest liability in the scheme presented here is the 
same as that of the proposal of Quint and Gabrielse 
[322], namely, the weakness of the parametric coupling 
between proton magnetic moment and axial motion. 
For (1)^/ 217=1 MHz, we find zo(proton) = 7 1 nm. With 
/= 1 A, and p = 1 mm, we find /3i/2'jt ~ 0.15 Hz. For 
these same conditions and d = 0.5 mm, we have 
?,xeh ~ 27 ms. At a field of 5 T, coo/2'n ~ 106 MHz. At 
very high resolutions we want resonance linewidths of 
less than 1 Hz which, in turn, requires long resonance 
times for the proton spin flip. Therefore the long times 
required for ^coupling need not be a serious liability. 

We have not considered details of proton, antiproton, 
or 'Be* transfer into the trap, however this might be 
accomplished by adapting a scheme similar to that 
described in Sec. 4.1. The trapping arrangement we 
show in Fig. 12 is essentially the same as the coupled 
trap idea of Ref. [48]. By removing the common endcap 
in the coupled-trap scheme, we arrive at the situation 
depicted in Fig. 12. In either case, we can show that the 
coupling between ions is approximately given by Eq. 
(154) where d is the overall distance between ions 1 and 
2. The potential advantage of the scheme described 
here is that d can probably be made smaller than 
in a coupled trap, thereby reducing 4xch- We have 
assumed the use of 'Be* ions, but many other ions 
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would work. A potential advantage of 'Be* ions is that 
the axial potential wells required to make Wi = Wi are 
not as different as for other choices of ion 2. In a similar 
spirit, we could measure the proton (or antiproton) 
magnetic moment in terms of other atomic parameters. 
For example, the proton spin flip frequency vs 'Be* spin 
flip frequency could be combined with a separate 
measurement of the 'Be* spin flip frequency compared 
to its cyclotron frequency to yield a measurement of the 
proton spin flip frequency to 'Be* cyclotron frequency 
in the same magnetic field. Coupled with an accurate 
value of the electron to 'Be* mass ratio, these measure- 
ments yield an accurate value of gp. An important 
systematic effect to consider in these measurements is 
magnetic field inhomogeneity. Field homogeneity could 
be checked to high accuracy by moving the 'Be* ion to 
various locations in the trap, thereby mapping the 
magnetic field. Equally important in the proton/antipro- 
ton comparison is to insure that the proton and anti- 
proton are in the same location. It appears this could be 
accomplished by insuring the 'Be* is always at the same 
location (using optical means) and adjusting the trap 
potentials to always yield the same values of w^ and 
8a)i,2. By this method, it appears that accuracy signifi- 
cantly better that 1 part in 10' could be achieved. 

6.4 Quantum State Manipulation of Mesoscopic 
Mechanical Resonators 

Much of what we have discussed concerns manipula- 
tion of the mechanical oscillation of atoms or atomic 
ions. It is perhaps interesting to speculate on the possi- 
bility of applying similar techniques to the manipulation 
of macroscopic or mesoscopic mechanical resonators in 
the quantum regime [119,120]. If studies can be per- 
formed at the quantum level, new sensors at the single 
phonon level could be built [324]. One approach is to 
make the mechanical resonators small enough and 
the temperature low enough that Aw>k^T [325]; this 
approach may obviate the idea discussed here. The idea 
we examine here is an extension of the idea of coupling 
the oscillatory motion of two ions together; here we 
consider coupling the motion of an ion (or the COM 
mode of a collection of ions) to a single mode of 
mesoscopic mechanical resonator. Although current 
technology appears to prohibit performing such experi- 
ments at the quantum level, the development of 
mesoscopic resonators is rapid and such experiments 
may be possible in the future. Similar considerations 
regarding coupling of ions to piezoelectric resonators 
were discussed in Ref. [48]. 

To be specific, we consider the situation sketched in 
Fig. 13. The mechanical resonator is assumed to be 
a silicon beam resonator, fixed at both ends; we will 




Fig. 13. Schematic diagram of an ion coupled to mechanical 
resonator. The bridge-type cantilever is essentially the same as that 
reported by Cleland and Roukes [325], but is assumed to support a 
metallized sphere supporting charge q\. An ion ^2 is confined by the 
potentials on the metallized sphere and additional electrodes shown 
schematically as A, B, and C. The Coulomb coupling between the 
charges provides the coupling between the two harmonic oscillators. 

assume the conditions realized in the experiments of 
Cleland and Roukes [325]. The beam resonator has a 
length D , thickness in the z direction equal to Az and a 
thickness in the x direction (out of the plane of Fig. 1 3) 
equal to Ax. If we take the conditions of Fig. 3 of 
Cleland and Roukes [325] as a guide, we have 
D = 1.1 (xm, Az = 0.8 (xm, and Ax = 0.33 |jum. We 
assume that the beam is metalized near its center and 
can support a charge q ^. For simplicity, we will assume 
the metalization is confined to a spherical shell of radius 
R |, but the exact geometry is not so important. Neglect- 
ing the dielectric effects of the Si beam (7=4 K), the 
capacitance of the metallized sphere is approximately 
equal to C^ = 4TT€oRi. The sphere can therefore support 
a charge qi = ChVh where Vb is the potential on the 
sphere. An atomic ion of charge ^2 and mass OT2 is 
trapped by a combination of electrodes at a distance d 
from the cantilever. One of the electrodes is the 
charged, metallized sphere; the other electrodes are 
indicated schematically as A, B, and C. To estimate the 
mass m2 of the mechanical resonator we assume it has 
a mass equal to half of the beam's mass and that this 
mass is concentrated in the metallized sphere. For the 
conditions assumed here, and p(Si) = 2.33 g • cm~^, we 
find nil ~ 2.4X 10''^ g. We arrange the trapping poten- 
tials so that the z -oscillation frequency of the ion is 
equal to the beam oscillation frequency; in this case, we 
realize two coupled oscillators as described in the 
previous section. If we take d = 5 |jum, Vb= 1000 V, 
(D/2tt- = 70MHz, ^2 = ? (proton), m2 = 9 u ('Be*), we 
find from Eq. (155) that /exch = 6.4 (xs. To reliably work 
in the quantum regime of the mechanical oscillator, /exch 
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must be smaller than t', the time for the mechanical 
oscillator to make a transition from its \n = 0) ground 
state to the |m = 1) state. If we assume the conditions of 
Cleland and Roukes [325] where T = 4 K, and 
2 = 2X10^ Eq. (64) gives f * = 0.0382 (xs. This is 
clearly too short to sympathetically cool the beam 
resonator's mode to the ground state, although some 
cooling could be achieved as outlined by in Ref. [48]. In 
the future, it might be possible to make the ratio f Vfexch 
larger than 1, perhaps through higher 2's and lower 
ambient temperatures [or perhaps if more exotic ion 
species with Z2 » 1 become available (Sec. 5.2)], it 
may be possible to manipulate the quantum motion of 
mesoscopic mechanical resonators by these techniques. 

7. Summary/Conclusions 

We have attempted to identify some of the important 
practical effects which must be taken into account in 
order to create arbitrary, entangled quantum states of 
trapped ions. We have taken a "passive stabilization" 
approach in that we try to anticipate and suggest ways to 
guard against the physical effects causing decoherence. 
Ultimately, complicated manipulations, such as lengthy 
quantum computations, are expected to benefit from 
and/or require some form of active error correction. 
Indeed, some of the near-term future experiments will 
probably demonstrate some of these schemes. In 
Sec. IV, we have listed some of the potential sources of 
decoherence in trapped ion experiments; here we 
speculate on what appear to be the most important of 
these. 

Motional decoherence is discussed in Sec. 4.1. 
In the NIST single 'Be"" ion experiments [17,21,45,47, 
131,132,21 1], heating appears to be the most important 
source of decoherence, primarily because, at the present 
time, its source is still unidentified. Various possibilities 
were discussed and, although the heating may be caused 
by some fundamental effect, we speculate that it is 
probably caused by some, as yet undetected, source of 
added electric field noise. (Because the current experi- 
ments employ electronic filtering at the motional 
frequencies, direct observation of fluctuating potentials 
on the electrodes has been precluded.) Future experi- 
ments will be able to resolve this. Moreover, once this 
source of noise is understood, the ion becomes an 
extremely sensitive detector of fluctuating potentials 
appearing on the electrode surfaces. In any case, it will 
probably be desirable to eventually operate the ion 
trap at cryogenic temperatures in order to, for example, 
reduce the effects of ion loss due to background 
gas collisions. A cryogenic environment will have the 
added benefit of reducing sources of electronic noise 
and associated heating. 



Various sources of internal-state decoherence are 
discussed in Sec. 4.2. In the current experiments, 
decoherence is dominated by fluctuations in magnetic 
fields acting on qubit transitions which are strongly field 
dependent. It appears that this problem can be highly 
suppressed by the use of magnetic shielding and, eventu- 
ally, use of qubit transitions which are insensitive 
to magnetic field to first order. We expect that internal 
state decoherence will be negligible compared to 
motional decoherence. 

Decoherence induced during application of the logic 
pulses (Sec. 4.4) may be the most troublesome. Many of 
the sources of this type of decoherence are primarily 
technical, for example, caused by intensity fluctuations 
in the laser pulses which induce transitions, or 
stray light impinging on ions not directly addressed. 
The more fundamental causes of decoherence are 
(1) coupling to internal or motional states other than the 
intended ones (spectator level problem), (2) coupling to 
unintended motional modes (cross-mode coupling), and 
(3) fluctuations in the Rabi rates due to excitation of the 
3L-\ unused motional modes (Debye-Waller factors). 
The first two of these effects appear to be a question of 
speed. Coupling to unwanted spectator modes and 
cross-mode coupling can always be avoided by making 
the operations slow enough that the extraneous 
couplings are suppressed by spectral selection. This has 
the negative effect of allowing more time for decoher- 
ence and increasing the required time for a given 
computation. Fluctuations caused by fluctuating Debye- 
Waller factors are, in principle, reduced as the number 
of ions increases because of the averaging effects of 
many modes [Eq. (128)]. However, it is also likely that 
excitation of these modes is harder to avoid as the 
number of ions (and unwanted spectator modes) 
increases. It will therefore be desirable to laser-cool all 
modes to the zero-point state. 

As the number of trapped ions increases, it will be- 
come increasingly difficult to avoid these three types of 
effects. Therefore, some sort of multiplexing scheme 
will be necessary when large numbers of qubits are 
involved. In Sec. 5.1, we have presented one possible 
solution where the ions are broken up into smaller 
numbers of independent groups or registers. The ions 
are then connected by moving ions around between 
registers. It may also be advantageous to multiplex 
quantum information within multiple internal states of 
ions; this is briefly described in Sec. 5.2. 

Stimulated-Raman transitions between long-lived 
qubit states (such as ground-state hyperfine levels) ap- 
pear to offer significant advantages over single-photon 
optical transitions. Single -photon transitions require 
high laser frequency stability, whereas stimulated- 
Raman transitions require only high relative frequency 
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stability between the Raman beams, which is techni- 
cally easier to accomplish. Stimulated-Raman transi- 
tions also provide the ability to select the magnitude and 
direction of the effective k vector (k^n) by choosing 
different directions for the beams. This has the advan- 
tage of, for example, allowing k^ff to be parallel to the 
axis of ions in a linear trap, thereby suppressing 
coupling to radial modes, while still allowing the Lamb- 
Dicke parameter to be controlled by adjusting the angle 
of the beams. At the same time, spatial selection of ions 
along the axis of the trap can be good since each Raman 
beam can be at a relatively steep angle relative to the 
trap axis. With single photon transitions, selection of 
ions and modes by appropriately choosing the direction 
of A: can only be obtained for the radial modes in a linear 
trap. Stimulated-Raman transitions have the potential 
disadvantage of inducing significant ac-Stark shifts 
(Sec. 4.4.3). However, for anticipated operating condi- 
tions, the effects of Stark shifts can be suppressed 
relative to the effects of laser amplitude fluctuations, 
which affect both single-photon and stimulated-Raman 
transitions. 

Decoherence in the ion trap system can probably be 
substantially reduced over what has been obtained in 
experiments so far. How far this reduction can be carried 
is an issue which must be resolved experimentally. 
Decohering effects may eventually be controlled to such 
a level that fault-tolerant error correction schemes might 
be employed to achieve computations of arbitrary 
length. This may only require a single operation fidelity 
of 10'^ to 10'*^ [185,326]. If this condition is met, speed 
will become an important issue because of the poten- 
tially large amount of overhead (increase in number of 
required qubits and operations in fault tolerant 
schemes). As discussed in Sec. 4.4.6, the Rabi rate for 
any operation is limited to approximately the motional 
mode frequency. In principle, mode frequencies can be 
substantially increased beyond what is currently 
achieved (= 10 MHz). An obvious direction to pursue is 
to make smaller traps with higher trapping potentials; 
however, this aggravates the problem of addressability, 
and will increase the coupling of the ions to the elec- 
trodes thereby increasing decoherence. The optimum 
conditions must, again, be resolved experimentally. 



8. Appendix A. Entangled States and 
Atomic Clocks 

We compare the use of entangled vs nonentangled 
states in an atomic clock under conditions and con- 
straints different than those considered in Sec. 3.4. We 
assume that the resources available are a given number 
of atoms L and a total observation time which is longer 



than any other time scale in the problem. We assume 
that decoherence of the internal states is negligible 
during the Ramsey free precession time Tr. To make an 
atomic clock, we want to steer, or "lock," a reference 
oscillator ( "local oscillator") to the atoms' resonance. 
Typically, we can find a local oscillator whose rms 
frequency fluctuations Acolo(t) over short averaging 
times T are smaller than the intrinsic fluctuations given 
by projection noise. However, for long times, the 
frequency fluctuations of the local oscillator are worse 
than those given by projection noise. (If this were not the 
case, the local oscillator would be a better clock, obviat- 
ing the need for the atoms.) By measuring the atomic 
populations after each application of the Ramsey radia- 
tion [3,9,100], we can generate an error signal which 
steers the local oscillator to the center of the atomic 
resonance with a servo time constant Ts. The servo can 
make a correction after a few independent measure- 
ments on the atoms; for simplicity, we assume Ts <^ Tr 
(Eq. (A4) below). To make the error signal as large as 
possible, we want the atomic linewidth Aco^ as small as 
possible. The atomic linewidth can be expressed as 
AcOs, = it/(L^^^~^^Tr), where €=1/2 for nonentangled 
states and e=l for entangled states [3,9,100]. There- 
fore, we want Tr as large as possible. However, if we 
make Tr too large, the local oscillator fluctuations 
A&)lo(7r) will be larger than Aco a, thereby giving no 
useful signal. This is the constraint which tells us 
whether entangled or nonentangled states are more 
useful. 

To analyze this problem in detail, we make the follow- 
ing assumptions: (1) We assume the (free-running 
or unlocked) local oscillator has fluctuations over 
averaging time r equal to 



Aa^oir) = Ct". 



(Al) 



This is assumed to hold over the range of values of r 
which include Tr and Ts. We will consider only values of 
n > -1/2, since, for n < -1/2, the local oscillator will 
better than the atoms for r sufficiently long. (2) We 
assume the intrinsic atomic clock stability, limited by 
projection noise, is given by (Sec. 3.4) 



A(Wn,eas(T) = " 



(A2) 



(3) For a given servo time constant Ts, we assume 

Awlo(Ts) = KiAwmc^siTs) (A3) 

Normally, in the case where Aco^o < Aw^eas for short 
times T, we would think of adjusting Ts so that /Ti ^ 1 ; 
that is, the local oscillator is locked to the atomic reso- 
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resonance at or before the time t when the local oscil- 
lator fluctuations become worse than the projection 
noise. However, when Ki> 1, the locked local oscillator 
stability can eventually reach the stability given by 
projection noise (approaching it as 1/t) given adequate 
servo gain [327-329]. Therefore we will allow K, to be 
larger than 1 . 

Along with these assumptions, we impose two con- 
straints: (1) The servo time constant Ts must be longer 
than the sampling time Tr. (We will assume Tr is much 
larger than the dead time, that is, the time for optical 
pumping, detection, etc.) This is expressed by the condi- 
tion 



Ts = K2TR 



K2>1. 



(A4) 



For entangled states, each measurement gives one of two 
possible values (Sec. 3.4). Therefore, K2 must be larger 
for entangled states than nonentangled states. (2) For 
T = Tr, the atomic linewidth must be greater than the 
local oscillator fluctuations or else the error signal used 
to correct the local oscillator frequency is ambiguous. 
This is expressed by the condition 



Aco^iTR) = T7/(L<^^-"rR = K.Acoi^o (Tr), K,>1. 



(A5) 



Eq. (A5) gives the sampling time Tr as a function of K^ 



Tr = 



CK^L^ 



(A6) 



Plugging this into Eq. (A2), we find that the stability of 
the locked local oscillator at long times t is 



CK, 



1 (/ie+1/2) 

2(/i+l T ~ n+1 



1 



77= ■ (A7) 



This is the main result of Appendix A. From this expres- 
sion, we see that Aw^eas is smaller for entangled states 
when n > 0; however, the gain is not as significant as 
when we assumed Tr to be fixed. We also see that 
Awmeas IS Smaller for nonentangled states when n < 0. 
These results are due to the constraint expressed in 
Eq. (A5). From the above, we also find 



K,: 



K, 



■V 



(A8) 



For nonentangled states on large numbers of atoms, Ki 
and Tr are both much larger than for entangled states, 
and a very long time may be required to achieve the 
intrinsic atomic clock stability given by Eqs. (A2) and 
(A7). therefore, as a practical issue, we may wish to 
constrain ^"1 to be equal to 1 . In that case we find, from 
Eqs. (A3) and (A5) two values of Rr. We must take the 
smaller of these, which results from the value derived 
from Eq. (A3) (with Ki= 1 ). Plugging this value of Tr 
into Eq. (A2), we find the stability of the locked oscilla- 
tor to be 



CKT"^ 



/(2^\ 



1 



V~T- 



(A9) 



For large L, entangled states will give a smaller value of 
Awmeas for all valucs of m > -1/2; however this result is 
not fundamental and simply comes from the constraint 
that Ki is equal to 1 . 



9. Appendix B. Master Equation for tlie 
Density Matrix of a Radiatively 
Damped Harmonic Oscillator 

Equation (62) is the master equation for the density 
matrix p of a single harmonic oscillator, including 
radiative damping terms. This equation is equivalent to 
the following system of coupled first-order differential 
equations for the diagonal and off-diagonal number- 
state matrix elements of p : 

{n\p\n) = p„„ = 7(M+l)(n+l)p„+i„+i 

-y{2nn+n+n)p„„+yrinp„.i „_, 



{m\p\n) = p^„ = y(n+ 1 )V(m+ l)(n+ l)p„+i „+i 
+ 7nVmMp„_i„_i 
-^y[2n(m+n+l)+(m+n)]pmn 



(Bl) 



Some special cases of these equations were given in Eq. 
(63), where they were used to estimate the time for an 
ion to make a transition from the state 10). Cohen- 
Tannoudji [330] treated a similar system, in which the 
incoherent excitation was absent (n = 0), but a 
monochromatic perturbation was present. The time 
derivative of the thermal average of the number of 
vibrational quanta is 
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2 nPnn = y(n+l) 2 n(n+l)p„+i 



-yn 2 np„, 

n = 



-y{2n+l) 2 n^ Pnn 

n = 

H-yn 2 n^ Pn-\n-\ 

n= 1 

The sums can be simplified to 

CO CO 

2 n(n+l)p„+i„+i = 2 (n+l-l)(n+l)p„+i„+i 

CO 

n = 

-2 (m+1)p„+i„+i 



n=0 n=0 



and 



(B2) 



(B3) 



2 nVn-in-i= 2 [(n-l)^+2n-l]p„_i„_i 

n= 1 n= 1 

= E [(n-l)'+2(«-l)+l]p„_,„_i 

n= 1 

CO 03 CO 

= 2 "^ P™ +2 2 «P™ +2 P"" 



n=0 n=0 



= 2 "^ P™ +2 2 «P"" +1- 
Collecting terms, we have 

2 «p™= [r(M+i)-r(2M+i)+'yM] 2 «V"^ 

+ [-y{n+l) -yn+2yn]2j npnn+ yn 



(B4) 



= -72 «P™+ T": 

n = 



-^ = -y{n)+yn. 



(B6) 



(B5) 



In the steady state, {n) = n, independent of y (see also, 
Ref. [121], p. 97). 
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Note added in proof 

Some recent experiments are relevant to the subjects 
of this paper. In addition to the work at NIST [72,191] 
and Garching [219], other groups have realized 
Doppler-cooled "arrays" of ions which might be used as 
quantum registers [331-335]. The motional modes of 
two trapped ions have recently been laser-cooled to the 
ground state [336] and differentially addressed to make 
approximate Bell states [337]. Recent experiments on 
neutral atoms confined in optical lattices have employed 
sideband laser cooling using stimulated-Raman transi- 
tions [338,339]; these experiments have achieved cool- 
ing to (n)=0 with high efficiency [338,340]. Quantum 
state engineering and manipulation, and schemes for 
quantum logic will follow [340-342]. 
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